Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,368,189 A→G 5.3% A121A (GCA→GCG 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,1890AG 5.3% 497.8 / 12.6 150A121A (GCA→GCGyhcGDUF1016 family protein in the PD‑(D/E)XK nuclease superfamily
Reads supporting (aligned to +/- strand):  ref base A (85/56);  new base G (3/5);  total (89/61)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.73e-01

Only 100 of 153 total aligned reads displayed.

AAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTC  >  NC_000913/3368040‑3368337
                                                                                                                                                     |                                                                                                                                                    
aaGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCa                                                                                                                                                      <  2:289586/150‑1 (MQ=255)
  ggATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGa                                                                                                                                                    >  1:2749931/1‑150 (MQ=255)
  ggATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAAGTAACACCATTGGGACTTCCGCAGTCACCGTCTGCAGa                                                                                                                                                    >  1:2749917/1‑150 (MQ=255)
   gATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGaa                                                                                                                                                   <  1:2649989/150‑1 (MQ=255)
   gATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGaa                                                                                                                                                   <  2:2057925/150‑1 (MQ=255)
        aTGTCTATGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCTGCAGACGGTGTCTGCGGAATTccc                                                                                                                                              <  2:2408178/150‑1 (MQ=255)
        aTGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAACCCGTCCGCCCAATGAGGCTTTTTTCCCTCTCTTTTCCACATGTGGAAATTCACCAGACAAGGTCTGGCGAATTAACACCATTGGGAATTCTGCAGACGGTGTCTGCGGAATTccc                                                                                                                                              <  2:2052989/150‑1 (MQ=255)
        aTGTATAAGGTATAAGCGTGGGTTCTCGGCAAACAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCTGCAGACGGTGTCTGCGGAATTccc                                                                                                                                              <  2:1622080/150‑1 (MQ=255)
           tCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTGTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATc                                                                                                                                           >  2:849923/1‑150 (MQ=255)
            cTACGGTATAAGCGTGGGTTCTCTGCAAAAATCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACAGGTTGAAATTCACCAGACAATGTCTGGCGTATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCg                                                                                                                                          <  2:2731393/150‑1 (MQ=255)
            cTAAGGTATCTGCGTGGGTTCTCTGCAAAAAACTTCCGCCAAATGAGGCTTTTTTTCCTCTTTTTTCAACATGTTGAAATTCTCCAGACAATGTCTGGCGAATTTACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCg                                                                                                                                          <  2:2731477/150‑1 (MQ=255)
            cTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCg                                                                                                                                          <  2:7068/150‑1 (MQ=255)
             tAAGGTATAAGCGTGGTTTCTCTGAAAAAAACTTACGCCAACTGAGTCTTTTTTACCTCTTTTTTAAACATTTTGAAATTCACCAGAAAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCgg                                                                                                                                         <  1:1763604/150‑1 (MQ=255)
             tAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCgg                                                                                                                                         >  2:331892/1‑150 (MQ=255)
                ggTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCaa                                                                                                                                      <  2:266123/150‑1 (MQ=255)
                 gTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCaaa                                                                                                                                     >  2:1238247/1‑150 (MQ=255)
                 gTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCaaa                                                                                                                                     >  1:2333942/1‑150 (MQ=255)
                 gTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCaaa                                                                                                                                     >  1:2295647/1‑150 (MQ=255)
                 gTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATAAACACCATTGGGAATACCGCAGACACCCTCTGCAGAATTCCCATCGGCaaa                                                                                                                                     >  1:1677578/1‑150 (MQ=255)
                       gCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTg                                                                                                                               >  2:1048096/1‑150 (MQ=255)
                        cGTGGGTTCTCTGCAAAAAACTTCCGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTgg                                                                                                                              <  1:2261398/150‑1 (MQ=255)
                        cGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGATTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTgg                                                                                                                              <  2:2265589/150‑1 (MQ=255)
                            ggTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGa                                                                                                                          >  1:1966881/1‑150 (MQ=255)
                              ttCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACg                                                                                                                        >  2:674986/1‑150 (MQ=255)
                                ctctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAcagc                                                                                                                      <  1:210556/150‑4 (MQ=255)
                                 tctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTg                                                                                                                     >  1:1163719/1‑150 (MQ=255)
                                 tctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTg                                                                                                                     >  2:1033837/1‑150 (MQ=255)
                                  ctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCACTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTgc                                                                                                                    <  1:2075095/150‑2 (MQ=255)
                                   tGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTc                                                                                                                   >  1:415856/1‑150 (MQ=255)
                                   tGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTc                                                                                                                   >  1:149925/1‑150 (MQ=255)
                                    gCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTACCATCGGCAAAGATTTGGCAGACGCTGTCt                                                                                                                  >  1:478636/1‑150 (MQ=255)
                                     cAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAAATCCCATCGTCAAAGATTTGGCAGACGCTGTCTg                                                                                                                 >  2:2869879/1‑150 (MQ=255)
                                      aaaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGACTTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGc                                                                                                                >  1:2326284/1‑150 (MQ=255)
                                      aaaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGc                                                                                                                >  1:89786/1‑150 (MQ=255)
                                       aaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCCGACCCCGTCTGCAGACTTCCCATCGGCAAAGATTTGGCAGACGCGGTCTGcc                                                                                                               >  1:2013470/1‑150 (MQ=255)
                                       aaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGcc                                                                                                               <  1:2677589/150‑1 (MQ=255)
                                       aaaaaCTTACGCCAAAGGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGcc                                                                                                               <  1:2603061/150‑1 (MQ=255)
                                         aaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACAGCGTCTGCCaa                                                                                                             <  2:1196810/150‑1 (MQ=255)
                                          aaCTAACGCCAAATGAGGCTTTTTTACCTATTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGTAATTCCGCAGACAGCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCaaa                                                                                                            <  1:1340568/150‑1 (MQ=255)
                                           aCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAAt                                                                                                           >  2:1544923/1‑150 (MQ=255)
                                           aCTTACGCCAAATGAGGCTTTTTTACCTATTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAAt                                                                                                           <  1:1060792/150‑1 (MQ=255)
                                             ttACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCt                                                                                                         >  2:317274/1‑150 (MQ=255)
                                             ttACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCt                                                                                                         >  1:2023923/1‑150 (MQ=255)
                                                 gCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTcc                                                                                                     >  1:65928/1‑150 (MQ=255)
                                                 gCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCGGTCTGCCAAAACTTTcc                                                                                                     >  1:1365212/1‑150 (MQ=255)
                                                  ccAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTccc                                                                                                    <  2:2694009/150‑1 (MQ=255)
                                                  ccAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTccc                                                                                                    >  1:2909871/1‑150 (MQ=255)
                                                    aaaTGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAAATAACACCATTGGGAATTCCGCAGACACCGGCTGCAGCATTCCCATCGGCAAAGATTTGGCAGACCCGGTCTGCCAACTCTTTCCCCt                                                                                                  >  1:1292272/1‑150 (MQ=255)
                                                           gCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGTATTTCCGCAGACACCGTCTGAAGAATTCCCATCGGAAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAg                                                                                           <  1:972274/150‑1 (MQ=255)
                                                             ttttttACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGt                                                                                         >  2:1167686/1‑150 (MQ=255)
                                                             ttttttACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAAGTCCGCAGACACCGTCTGCCGAATTCCCATCGGCAAAGATGTGGAAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGt                                                                                         >  2:1909070/1‑150 (MQ=255)
                                                                tttACCTCTTTTTTCAACATGTTGAAATTCACCAGACATTGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCa                                                                                      <  1:2547962/150‑1 (MQ=255)
                                                                tttACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCa                                                                                      <  2:1645918/150‑1 (MQ=255)
                                                                tttACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCa                                                                                      <  2:1053148/150‑1 (MQ=255)
                                                                     ctctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCCGACGCTGTCTGCCAAAACTTTCCCCTTGCCCCGGGCCCCCTAc                                                                                 >  1:1451159/1‑150 (MQ=255)
                                                                     ctcaTTGTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATATGTCCCCTTGCCCAGGACCCCCTAc                                                                                 >  2:1776751/1‑150 (MQ=255)
                                                                      tctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACg                                                                                >  2:935648/1‑150 (MQ=255)
                                                                      tctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACg                                                                                >  1:416654/1‑150 (MQ=255)
                                                                        ttttttCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCTGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCACATCTTTGCCAGTGCCCAGGGCCGCCTGCggc                                                                              >  1:1834556/1‑148 (MQ=255)
                                                                           tttCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGACTTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGt                                                                           >  2:1147287/1‑150 (MQ=255)
                                                                           tttCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGt                                                                           >  1:2165900/1‑150 (MQ=255)
                                                                                aCATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACt                                                                      >  2:387335/1‑150 (MQ=255)
                                                                                    gTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTcc                                                                  <  1:164949/150‑1 (MQ=255)
                                                                                    gTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTcc                                                                  <  1:164952/150‑1 (MQ=255)
                                                                                     ttGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCCGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCg                                                                 >  2:874513/1‑150 (MQ=255)
                                                                                           ttCGCCAGACATTGTCTGGTGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTaaaa                                                           >  1:728078/1‑150 (MQ=255)
                                                                                           ttCGCCAGACATTGTCTGGTGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTaaaa                                                           >  1:728220/1‑150 (MQ=255)
                                                                                           ttCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTaaaa                                                           >  2:1499099/1‑150 (MQ=255)
                                                                                              aCCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACg                                                        >  2:247519/1‑150 (MQ=255)
                                                                                              aCCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACg                                                        <  2:1677578/150‑1 (MQ=255)
                                                                                                 aGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTg                                                     <  2:2909871/150‑1 (MQ=255)
                                                                                                   aCAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGAc                                                   <  1:674986/150‑1 (MQ=255)
                                                                                                     aaTGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGc                                                 <  1:1279439/150‑1 (MQ=255)
                                                                                                     aaTGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGc                                                 >  2:1828394/1‑150 (MQ=255)
                                                                                                       tGTCGGGCGAATTAACACCCTTGGGATTTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGccc                                               <  1:2413812/150‑1 (MQ=255)
                                                                                                        gTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCAACGTCCGTTTACTTTACGTTAAAAACGCTGACGCCCg                                              >  1:2367576/1‑150 (MQ=255)
                                                                                                           tGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGGATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAg                                           >  2:1951306/1‑150 (MQ=255)
                                                                                                            ggCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGGTTACGTTCCGTTAAAAAAGCTGACGCCCGCAgg                                          >  1:2125678/1‑149 (MQ=255)
                                                                                                             gCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCt                                         <  1:1033837/150‑1 (MQ=255)
                                                                                                                 cttAACCCCCTGGGGAAGTCCGCAGACACCGTCTGCAGAATTCCCATCGGCACCGATTGGGCAGCCGCTGTCTGCCAAATCTTTCCCCGGGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTa                                     <  1:2855211/149‑1 (MQ=255)
                                                                                                                 aTTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTa                                     <  1:1350278/150‑1 (MQ=255)
                                                                                                                  ttAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTAt                                    >  2:723000/1‑150 (MQ=255)
                                                                                                                  ttAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTAt                                    >  2:1925966/1‑150 (MQ=255)
                                                                                                                  ttAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCAGTCTGCCAAATCTGTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTGTAt                                    >  2:2815475/1‑150 (MQ=255)
                                                                                                                    aaCACCCGTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGa                                  <  2:2013470/150‑1 (MQ=255)
                                                                                                                          aTTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAaga                            <  1:1923516/150‑1 (MQ=255)
                                                                                                                           ttGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTGTATGAACAagag                           >  2:2849895/1‑150 (MQ=255)
                                                                                                                           atGGGAATTCTGCAGACGGTGTCTGCGGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAagag                           >  2:2017308/2‑150 (MQ=39)
                                                                                                                                  ttCTGCAGACGGTGTCTGCGGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTcc                    >  2:2152101/1‑150 (MQ=255)
                                                                                                                                  ttCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTcc                    <  1:2528345/150‑1 (MQ=255)
                                                                                                                                      gCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCtg                <  2:2295647/150‑1 (MQ=255)
                                                                                                                                      gCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCtg                <  2:378176/150‑1 (MQ=255)
                                                                                                                                       cAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCtgt               <  1:1168618/150‑1 (MQ=255)
                                                                                                                                       cAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCtgt               <  2:175177/150‑1 (MQ=255)
                                                                                                                                         gACACCGTCTGCAGAATGCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTgg             <  2:1353219/150‑1 (MQ=255)
                                                                                                                                          acacCGTCTGCAGAATTCCCATCGGCAAAGATTTTGCAGACGCTGTCTGCCAAAACATTCCCCTTGCCCAGGTCCACCGACGTCCCGTTACTTTACGTTAAAAACGCTGAAGCCCGCAGCGTTTATGAAAAAGAGACGCTCCGCGGTGGc            >  1:439256/1‑150 (MQ=255)
                                                                                                                                              cGTCTGCAGAATTCCCATCGGCAAAGATTTTGCAGAAGCTGTCTGCCAAATCTTTACCCATGCCCAGGTCCACCAACGTCCGTTTTATTTCCGTTAAAAAAGCTGACGCCCCCAGCGTTTATGTAAAAGAGACGCTCCGCGGTGGCGGGt        >  1:1415726/1‑150 (MQ=255)
                                                                                                                                               gTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTc       >  2:1833139/1‑150 (MQ=255)
                                                                                                                                                  tGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGt    >  2:748178/1‑150 (MQ=255)
                                                                                                                                                    cAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTc  >  2:2114836/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTC  >  NC_000913/3368040‑3368337

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: