Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,368,220 C→G 5.6% L132V (CTG→GTG) ‡ 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,2200CG 5.6% 512.1 / 16.6 143L132V (CTG→GTG) ‡yhcGDUF1016 family protein in the PD‑(D/E)XK nuclease superfamily
Reads supporting (aligned to +/- strand):  ref base C (70/65);  new base G (4/4);  total (74/69)
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.71e-01

Only 100 of 146 total aligned reads displayed.

CTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGTT  >  NC_000913/3368072‑3368368
                                                                                                                                                     |                                                                                                                                                    
ctctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑Gc                                                                                                                                                      <  4:52657/149‑1 (MQ=255)
 tctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCt                                                                                                                                                     >  4:1102973/1‑149 (MQ=255)
   tGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGt                                                                                                                                                   >  4:1536074/1‑149 (MQ=255)
   tGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGt                                                                                                                                                   >  3:358903/1‑149 (MQ=255)
    gCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTc                                                                                                                                                  <  4:84459/149‑1 (MQ=255)
     cAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCt                                                                                                                                                 >  1:945869/1‑149 (MQ=255)
     cAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCt                                                                                                                                                 <  1:858922/149‑1 (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                                                                                                                                     <  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)
                                     ctctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTa                                                                                                                 >  4:427086/1‑149 (MQ=255)
                                      tctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTAc                                                                                                                >  3:910932/1‑149 (MQ=255)
                                      tctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTAc                                                                                                                <  1:1581340/149‑1 (MQ=255)
                                      tctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTAc                                                                                                                >  3:958482/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)
                                                              gccAGACATTGTCTGGTGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAAc                                                                                        >  2:1591047/2‑149 (MQ=255)
                                                                  gACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTg                                                                                    >  1:1435305/1‑149 (MQ=255)
                                                                     aaTGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAgg                                                                                                                         >  1:1449080/1‑109 (MQ=255)
                                                                     aaTGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAgg                                                                                                                         <  2:1449080/109‑1 (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                                                             <  4:340632/149‑1 (MQ=255)
                                                                                         cATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaaaaa                                                             >  3:824417/1‑149 (MQ=255)
                                                                                           ttGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAaga                                                           <  4:147807/149‑1 (MQ=255)
                                                                                           ttGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAaga                                                           >  4:1132108/1‑149 (MQ=255)
                                                                                           ttGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAaga                                                           >  1:263430/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)
                                                                                                          acacCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTgg                                            <  2:1206692/149‑1 (MQ=255)
                                                                                                                   gCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGt                                   >  3:540956/1‑149 (MQ=255)
                                                                                                                   gCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGt                                   >  2:1547035/1‑149 (MQ=255)
                                                                                                                   gCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGt                                   <  3:752126/149‑1 (MQ=255)
                                                                                                                     aGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGtt                                                                                                         <  1:999141/77‑1 (MQ=255)
                                                                                                                     aGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGtt                                                                                                         >  2:999141/1‑77 (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                     <  3:974072/149‑1 (MQ=255)
                                                                                                                                 ggCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGc                     <  2:1392465/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)
                                                                                                                                     aaGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAgaga                                                         <  4:1493647/109‑1 (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           >  3:906394/1‑149 (MQ=255)
                                                                                                                                           tGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGc           <  2:1427617/149‑1 (MQ=255)
                                                                                                                                           tGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGc           <  4:267790/149‑1 (MQ=255)
                                                                                                                                            ggCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAgaga                                                         <  4:646356/102‑1 (MQ=255)
                                                                                                                                            ggCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAgaga                                                         >  3:646356/1‑102 (MQ=255)
                                                                                                                                            ggCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCg          <  4:624787/149‑1 (MQ=255)
                                                                                                                                            ggCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCg          <  3:1356875/149‑1 (MQ=255)
                                                                                                                                               aGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGAcc       >  1:971945/1‑149 (MQ=255)
                                                                                                                                               aGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGAcc       >  3:1173784/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)
                                                                                                                                                     |                                                                                                                                                    
CTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAC‑GCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGTT  >  NC_000913/3368072‑3368368

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