Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,368,221 T→C 5.4% L132P (CTG→CCG) ‡ yhcG → DUF1016 family protein in the PD‑(D/E)XK nuclease superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,368,2210TC 5.4% 461.9 / 10.3 144L132P (CTG→CCG) ‡yhcGDUF1016 family protein in the PD‑(D/E)XK nuclease superfamily
Reads supporting (aligned to +/- strand):  ref base T (69/64);  new base C (4/4);  total (73/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 146 total aligned reads displayed.

TCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGTTT  >  NC_000913/3368073‑3368369
                                                                                                                                                     |                                                                                                                                                    
tctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCt                                                                                                                                                      >  4:1102973/1‑149 (MQ=255)
  tGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGt                                                                                                                                                    >  3:358903/1‑149 (MQ=255)
   gCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTc                                                                                                                                                   <  4:84459/149‑1 (MQ=255)
    cAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCt                                                                                                                                                  <  1:858922/149‑1 (MQ=255)
    cAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCt                                                                                                                                                  >  1:945869/1‑149 (MQ=255)
    cAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑ACTGTCt                                                                                                                                                  <  3:1499286/149‑1 (MQ=255)
     aaaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTg                                                                                                                                                 >  1:1392465/1‑149 (MQ=255)
        aaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACAGC‑GTCTGCCa                                                                                                                                              <  3:727626/149‑1 (MQ=255)
        aaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCa                                                                                                                                              <  2:687349/149‑1 (MQ=255)
        aaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑AGCGTCTGCCa                                                                                                                                              <  2:78990/149‑1 (MQ=255)
         aaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCTTTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑AGCGTCTGCCaa                                                                                                                                             <  3:1095220/149‑1 (MQ=255)
           cTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑AGCGTCTGCCAAAt                                                                                                                                           <  1:1261457/149‑1 (MQ=255)
            ttACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATc                                                                                                                                          >  3:300920/1‑149 (MQ=255)
               cGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCttt                                                                                                                                       <  2:326091/149‑1 (MQ=255)
                gCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTc                                                                                                                                      <  4:301262/149‑1 (MQ=255)
                gCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTc                                                                                                                                      <  2:1458728/149‑1 (MQ=255)
                 ccAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTcc                                                                                                                                     >  3:624787/1‑149 (MQ=255)
                  cAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTccc                                                                                                                                    <  4:531714/149‑1 (MQ=255)
                          gCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACAGC‑GTCTGCCaaa                                                                                                                                            <  1:604663/133‑1 (MQ=255)
                                 tACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCAc                                                                                                                     <  4:323100/149‑1 (MQ=255)
                                    ctctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGTAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCAc                                                                                                                     >  1:1364694/1‑146 (MQ=255)
                                    ctctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCAc                                                                                                                     <  2:1364694/146‑1 (MQ=255)
                                     tctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTAc                                                                                                                 >  3:958482/1‑149 (MQ=255)
                                     tctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTAc                                                                                                                 <  1:1581340/149‑1 (MQ=255)
                                     tctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTAc                                                                                                                 >  3:910932/1‑149 (MQ=255)
                                            tCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGtt                                                                                                          >  4:224549/1‑149 (MQ=255)
                                               aCATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTAc                                                                                                       >  3:564135/1‑149 (MQ=255)
                                                    ttGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTcc                                                                                                  >  1:1467955/1‑149 (MQ=255)
                                                    ttGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGCTTACTTTcc                                                                                                  >  2:152232/1‑149 (MQ=255)
                                                     tGAAATTCACCAGACAATGGCTAGCGAATTAACACCATTGGGAATTCCGCAGGCATCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCGAAATCTTGCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCg                                                                                                 <  2:646810/149‑1 (MQ=255)
                                                          ttCGCCAGACATTGTCTGGTGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTaaa                                                                                            >  1:32502/1‑149 (MQ=255)
                                                            cACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTaaaaa                                                                                          <  4:352093/149‑1 (MQ=255)
                                                             gccAGACATTGTCTGGTGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAAc                                                                                         >  2:1591047/2‑149 (MQ=255)
                                                                aGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCt                                                                                      <  4:472295/149‑1 (MQ=255)
                                                                 gACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTg                                                                                     >  1:1435305/1‑149 (MQ=255)
                                                                    aaTGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAgg                                                                                                                          <  2:1449080/109‑1 (MQ=255)
                                                                    aaTGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAgg                                                                                                                          >  1:1449080/1‑109 (MQ=255)
                                                                       gTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGccc                                                                               <  3:700070/149‑1 (MQ=255)
                                                                            gCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGc                                                                          >  2:145984/1‑149 (MQ=255)
                                                                             cGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCt                                                                         <  3:628462/149‑1 (MQ=255)
                                                                                aTTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCtttt                                                                      >  2:1447225/1‑149 (MQ=255)
                                                                                aTTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCtttt                                                                      >  2:1050353/1‑149 (MQ=255)
                                                                                 ttAACACCATTGGGAATTCCGCAGATACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTa                                                                     <  1:1057834/149‑1 (MQ=255)
                                                                                 ttAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTTCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTa                                                                     <  1:594910/149‑1 (MQ=255)
                                                                                 ttAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTa                                                                     <  4:39515/149‑1 (MQ=255)
                                                                                 ttAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTa                                                                     <  3:1007021/149‑1 (MQ=255)
                                                                                   aaCACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATg                                                                   >  4:1356875/1‑149 (MQ=255)
                                                                                    acacCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGa                                                                  >  1:1427617/1‑149 (MQ=255)
                                                                                     cacCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaa                                                                 >  2:583945/1‑149 (MQ=255)
                                                                                        cgatGGGAATTCTGCAGACGGTGTCTGCGGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaaaaa                                                              >  1:1186909/4‑149 (MQ=38)
                                                                                        cATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaaaaa                                                              >  3:824417/1‑149 (MQ=255)
                                                                                          ttGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAaga                                                            >  4:1132108/1‑149 (MQ=255)
                                                                                          ttGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAaga                                                            <  4:147807/149‑1 (MQ=255)
                                                                                          ttGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAaga                                                            >  1:263430/1‑149 (MQ=255)
                                                                                          ttGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAaga                                                            >  4:167740/1‑149 (MQ=255)
                                                                                            gggAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTCTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAgaga                                                          <  4:958482/149‑1 (MQ=255)
                                                                                                aTTCTGCAGACGGTGTCTGCGGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCt                                                      >  4:356423/1‑149 (MQ=255)
                                                                                                  tCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTcc                                                    >  3:106973/1‑149 (MQ=255)
                                                                                                       aGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAAAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCtgt                                               <  4:151272/149‑1 (MQ=255)
                                                                                                         acacCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTgg                                             <  2:1206692/149‑1 (MQ=255)
                                                                                                                  gCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGt                                    >  2:1547035/1‑149 (MQ=255)
                                                                                                                  gCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGt                                    <  3:752126/149‑1 (MQ=255)
                                                                                                                  gCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGt                                    >  3:540956/1‑149 (MQ=255)
                                                                                                                    aGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGtt                                                                                                          >  2:999141/1‑77 (MQ=255)
                                                                                                                    aGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGtt                                                                                                          <  1:999141/77‑1 (MQ=255)
                                                                                                                    aGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTc                                  <  3:1536074/149‑1 (MQ=255)
                                                                                                                        ttCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCa                              <  2:1359311/149‑1 (MQ=255)
                                                                                                                           ccATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACttt                                                                                                    >  4:893887/1‑76 (MQ=255)
                                                                                                                           ccATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCt                           >  3:1398589/1‑149 (MQ=255)
                                                                                                                            cATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTa                          >  1:799487/1‑149 (MQ=255)
                                                                                                                             aTCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGc                                                                                       >  4:1513539/1‑87 (MQ=255)
                                                                                                                             aTCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTg                                     <  3:1322513/137‑1 (MQ=255)
                                                                                                                             aTCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTg                                     >  4:1322513/1‑137 (MQ=255)
                                                                                                                             aTCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTag                         >  4:967038/1‑149 (MQ=255)
                                                                                                                               ggggaaaGATTTGGCCGAC‑AGCGTCTGCCAAATCTTTGCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaaa                                                                <  4:758694/104‑1 (MQ=37)
                                                                                                                               ggggaaaGATTTGGCAGAC‑AGCGTCTGCCAAATCTTTGCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaaa                                                                >  3:758694/5‑108 (MQ=37)
                                                                                                                                ggCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGc                      <  2:1392465/149‑1 (MQ=255)
                                                                                                                                ggCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGc                      <  3:974072/149‑1 (MQ=255)
                                                                                                                                 gCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCg                     >  3:1232529/1‑149 (MQ=255)
                                                                                                                                  gaaaGATTTGGCAGAC‑AGCGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCgg                    >  2:1503377/2‑149 (MQ=255)
                                                                                                                                  cAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCgg                    >  2:764466/1‑149 (MQ=255)
                                                                                                                                   aaaGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTc                                                                                                   <  1:1261994/69‑1 (MQ=255)
                                                                                                                                   aaaGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGc                   >  3:889046/1‑149 (MQ=255)
                                                                                                                                    aaGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAgaga                                                          >  3:1493647/1‑109 (MQ=255)
                                                                                                                                      gATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCaaa                <  3:224549/149‑1 (MQ=255)
                                                                                                                                        tttGGCAGACAGC‑GTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTAt                                                                    <  2:1188258/95‑1 (MQ=38)
                                                                                                                                        tttGGCAGAC‑AGCGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTAt                                                                    >  1:1188258/1‑95 (MQ=255)
                                                                                                                                         ttGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTg             <  2:1420351/149‑1 (MQ=255)
                                                                                                                                          tGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGc            <  2:1427617/149‑1 (MQ=255)
                                                                                                                                          tGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGc            >  3:906394/1‑149 (MQ=255)
                                                                                                                                           ggCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAgaga                                                          >  3:646356/1‑102 (MQ=255)
                                                                                                                                           ggCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAgaga                                                          <  4:646356/102‑1 (MQ=255)
                                                                                                                                           ggCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCg           <  4:624787/149‑1 (MQ=255)
                                                                                                                                           ggCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCg           <  3:1356875/149‑1 (MQ=255)
                                                                                                                                              aGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGAcc        >  3:1173784/1‑149 (MQ=255)
                                                                                                                                              aGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGAcc        >  1:971945/1‑149 (MQ=255)
                                                                                                                                              aGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGACCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCtgt                                               >  4:765524/1‑110 (MQ=255)
                                                                                                                                              aGAC‑GCTGTCTGCCAAATCGTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTTACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGAcc        >  2:304584/1‑149 (MQ=255)
                                                                                                                                                aC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCa      <  4:106973/149‑1 (MQ=255)
                                                                                                                                                     tGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGttt  <  4:1398589/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGTTT  >  NC_000913/3368073‑3368369

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑