Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,972,673 A→C 6.7% E51A (GAG→GCG) ‡ rffG → dTDP‑glucose 4,6‑dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,972,6730AC 6.7% 326.4 / 8.1 122E51A (GAG→GCG) ‡rffGdTDP‑glucose 4,6‑dehydratase
Reads supporting (aligned to +/- strand):  ref base A (45/67);  new base C (7/2);  total (53/69)
Fisher's exact test for biased strand distribution p-value = 3.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 132 total aligned reads displayed.

AGAAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGAAAC  >  NC_000913/3972525‑3972820
                                                                                                                                                    |                                                                                                                                                   
agaAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcgag                                                                                                                                                    >  2:342153/1‑150 (MQ=255)
 gaAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcgagc                                                                                                                                                   >  2:437139/1‑150 (MQ=255)
 gaAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcgagc                                                                                                                                                   >  2:1726916/1‑150 (MQ=255)
   aaaaTTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAgcgc                                                                                                                                                 >  2:1525469/1‑150 (MQ=255)
   aaaaTTCTGATAACAGGTGGTGCCGGGTTTATTGGATCGGCGCTGGTGCGTTATATCATAAACGATACGAGCGAGGCGCTGGTATTGTTCGAGAAGCTGACCTACGCCGGATCCCGGATCTCCCTGGCACCGTTCGCGCAAaacgagcgg                                                                                                                                                 >  2:131735/1‑149 (MQ=255)
       ttCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGCTCGCTTTg                                                                                                                                             <  2:1479100/150‑1 (MQ=255)
         cTGATAACAGGTGGTGCCGGGGTTATTGGCTCGGCGCAGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGAGAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGcc                                                                                                                                           <  1:945775/150‑1 (MQ=255)
         cTGATAACAGGCGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTAAATCATCAACGAAACGAGCGACGCGGTGGTAGCGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGcc                                                                                                                                           <  1:1814496/150‑1 (MQ=255)
           gctcaCAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCtt                                                                                                                                         <  2:1168760/146‑1 (MQ=255)
            aTAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCttt                                                                                                                                        >  2:364591/1‑150 (MQ=255)
             tAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTg                                                                                                                                       <  2:1768543/150‑1 (MQ=255)
             tAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTg                                                                                                                                       <  1:1765273/150‑1 (MQ=255)
             tAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTg                                                                                                                                       <  1:1609537/150‑1 (MQ=255)
             aaaCAGGTGGTGCCGGGATTATTGGCTCGGCGCTGGTGCGTTATATCATCATCGAAACGAGCGACGCGGTGGTAGGGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTg                                                                                                                                       <  1:226801/149‑1 (MQ=255)
               aCAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAATCCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTgag                                                                                                                                     <  1:1041686/150‑1 (MQ=255)
                  ggtggtGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGaaa                                                                                                                                  <  2:620212/150‑1 (MQ=255)
                   gtggtgCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAg                                                                                                                                 <  1:1003766/150‑1 (MQ=255)
                     ggtgCCGGGTTTATTGGCGCGGCGCTGGTGCGTTATATCATCAACGACACGAGCGACGCGGTGGTAGTGGGCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGtt                                                                                                                               <  1:964806/150‑1 (MQ=255)
                      gtgCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGATACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGGGCGCTTTGCCTGTGAGAAAGttt                                                                                                                              >  1:913589/1‑149 (MQ=255)
                      gtgCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTg                                                                                                                              <  2:679370/150‑1 (MQ=255)
                       tgCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGa                                                                                                                             <  1:626575/150‑1 (MQ=255)
                        gCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACAGGTCGCGCAAAGCGAGCGCTTTGCCATTGAGAAAGTTGat                                                                                                                            >  1:1053790/1‑150 (MQ=255)
                         ccGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTTGCACCGGTCGCGCCAAGCGATCGCTTTTACATTGAGAAAGTTtata                                                                                                                           >  1:1690730/1‑150 (MQ=255)
                         ccGGGCTTCTTGGCTCGCCGCTGGTGCGTTATATCATCACCGAAACGCGCGACGCGGGGGCAGTGGTCGAAAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAATGTTGata                                                                                                                           <  2:1226204/150‑1 (MQ=255)
                               ttATTGGCTCGGCGCTGGTGCGTTATATCATCAACTAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCg                                                                                                                     <  1:897926/150‑1 (MQ=255)
                               ttATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCACAGCGAGCGCTTTGCCTTTGAGTAAGTTGATATCTGCg                                                                                                                     >  1:311969/1‑150 (MQ=255)
                                tATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTGTGCCTTTGAGAAAGGTGATATCTGCGa                                                                                                                    >  1:1654629/1‑150 (MQ=255)
                                  ttGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATc                                                                                                                  <  1:1449246/150‑1 (MQ=255)
                                   tGGCGCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCg                                                                                                                 <  2:963480/150‑1 (MQ=255)
                                    ggCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCgg                                                                                                                >  2:1008466/1‑150 (MQ=255)
                                              tGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGGCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTgg                                                                                                      <  1:834066/150‑1 (MQ=255)
                                               ggTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGGCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTGGATATCTGCGATCGGGCAGAACTGGc                                                                                                     <  1:1382486/150‑1 (MQ=255)
                                                gTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCa                                                                                                    <  2:389470/150‑1 (MQ=255)
                                                 tGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGCACGCATTGCCCTTGAGAAATTTGAAAACAGCGTTCGGGCAGAAcatgccc                                                                                                   >  2:656267/1‑144 (MQ=255)
                                                    gTTATATCATCAACGAACCGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGGCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCAcgcg                                                                                                <  2:1670328/150‑1 (MQ=255)
                                                    gTTATATCATCAACGAAACGAGCGACGCGGGGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCAcgcg                                                                                                <  1:1849424/150‑1 (MQ=255)
                                                           catcaACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCAc                                                                                         <  1:850149/150‑1 (MQ=255)
                                                            ataaacgaaacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCATTTCCTTTGAGAAATTTTATATATGCGAACGGGACGCAAAGGCACGCGGATTCAca                                                                                        >  1:622204/4‑149 (MQ=255)
                                                                acgaaactaGCGACGCGGTGGTAGTGGTCTATAAGCTGACCTACGCCGGATACCTGATGTCGCTGGCACCTGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGTCAGAACTGGCACGCGTATTCACTGAGc                                                                                    <  1:540190/150‑1 (MQ=255)
                                                                   aaacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATc                                                                                 <  2:1719402/150‑1 (MQ=255)
                                                                    aacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCa                                                                                >  1:1179775/1‑150 (MQ=255)
                                                                      cgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGc                                                                              >  2:5782/1‑150 (MQ=255)
                                                                      cgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATATGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGc                                                                              >  1:671907/1‑150 (MQ=255)
                                                                            ccgcgGTGGTAGTGGTCGATACGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACt                                                                        <  2:256028/149‑1 (MQ=255)
                                                                               cgGTGGTAGTGGTCGATAAGCTTACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACtgtg                                                                     <  2:349650/150‑1 (MQ=255)
                                                                               cgGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCGTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACtgtg                                                                     >  2:678202/1‑150 (MQ=255)
                                                                                ggtggtAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTgtgt                                                                    >  1:1333392/1‑150 (MQ=255)
                                                                                ggtggtAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTgtgt                                                                    <  1:1316054/150‑1 (MQ=255)
                                                                                 gtggtAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTc                                                                   <  1:1409871/150‑1 (MQ=255)
                                                                                   ggtAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCAt                                                                 <  2:1028172/150‑1 (MQ=255)
                                                                                   ggtAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCAt                                                                 >  2:1110329/1‑150 (MQ=255)
                                                                                    gtAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCGTTGGGTAAGTTGGTATCTGCGATAGGGCAGAACTGGGCAGCGTTTTAAATGAGCATCAGCCAGACTGTGTCATg                                                                >  2:618179/1‑150 (MQ=255)
                                                                                     tAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCCAACCGACCGCGTTGCCTTTTAGAAAGTTGAAAAATGCGATCGGGCAGAACTGGACCCCGTAATAACATGCCAACAGCCCAACTGTGTCaatc                                                               >  2:1619708/1‑147 (MQ=255)
                                                                                      aGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCa                                                              <  1:1638462/150‑1 (MQ=255)
                                                                                        tGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCGGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATc                                                            <  1:537428/150‑1 (MQ=255)
                                                                                         ggTCGATAAGCTGACGTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCATTGCGTTTGAGCAAGTTGATATCTCCGATCGGGCAGAAATGGCACGCGTATTCAATGAGAATAAGCCAGACGGTGTCCTGCATCt                                                           >  2:1250189/1‑150 (MQ=255)
                                                                                          gTCGATAAGCTGACCTACGCCGGAAAGCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGAAATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTg                                                          <  1:122176/150‑1 (MQ=255)
                                                                                          gTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGCGCGCTTTGCCTTGGAGAAAGTTGATATCTGCGATCGGGCAGACCTGGCACGCGTATTCACTGAGCATCCGCCAGACTGTGTCATGCATCTg                                                          >  1:1055432/1‑150 (MQ=255)
                                                                                           gcgctAAGCTGCCCTACGCCGGAAACCTGATGTCGCTGGCACCGGGCGCGCAAAGCGAGCGCTTTGCCTGTGAGTAAGTTGTTCTCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCAGCTgg                                                         <  2:334109/146‑1 (MQ=255)
                                                                                            cGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATATGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGGGTCATGCATCTGGc                                                        >  1:1645528/1‑150 (MQ=255)
                                                                                            cGATAAGATGCCCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAACTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCAACTGGc                                                        >  1:1007522/1‑150 (MQ=255)
                                                                                             gATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCGTTGAGAAAGTTGATATCTGCGATCGGGCAGAACAGGCACGCGTATAAACTGAGCATAAGCCAGACTGTGTCATGCAACTGGcc                                                       >  2:194790/1‑149 (MQ=255)
                                                                                                aaGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGcc                                                    <  1:1008466/150‑1 (MQ=255)
                                                                                                 aGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGAAATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCg                                                   <  2:1689470/150‑1 (MQ=255)
                                                                                                  gCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCATCACTTGCACGCGTTTTACCTTACCATCAGCCAGACTGTGTCATGCTTCTGGCacgcga                                                  >  1:1172257/1‑145 (MQ=255)
                                                                                                  gCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGa                                                  <  1:342153/150‑1 (MQ=255)
                                                                                                     tacatACGCAGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATATGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAg                                               <  2:315376/146‑1 (MQ=255)
                                                                                                     gACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAg                                               <  1:1525469/150‑1 (MQ=255)
                                                                                                     gACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAg                                               <  1:1644667/150‑1 (MQ=255)
                                                                                                      aCCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGc                                              <  2:488973/150‑1 (MQ=255)
                                                                                                         tACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCATAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCAt                                           >  1:809592/1‑150 (MQ=255)
                                                                                                         tACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCCATCGGGCAGAACTGGCCCGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCAt                                           >  1:996700/1‑150 (MQ=255)
                                                                                                           cGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGt                                         <  1:326374/150‑1 (MQ=255)
                                                                                                           cGCCGGAAACCTGATGACGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGt                                         <  1:1730862/150‑1 (MQ=255)
                                                                                                            gCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGtt                                        <  1:626463/150‑1 (MQ=255)
                                                                                                            gCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGtt                                        <  1:364591/150‑1 (MQ=255)
                                                                                                              cGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGa                                      >  2:1028393/1‑150 (MQ=255)
                                                                                                                  aaCCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTTGCCGt                                  >  2:20318/1‑150 (MQ=255)
                                                                                                                      tGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTTGCACGCGTATTACCTGAGCATCAGCCAGCCGGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTa                              >  2:701411/1‑150 (MQ=255)
                                                                                                                       gATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGCCCGTTCTAt                             >  2:95491/1‑150 (MQ=255)
                                                                                                                       gATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTAt                             <  1:281959/150‑1 (MQ=255)
                                                                                                                         tGTCGCTGGCACCGGGCGGGCAAAGCGAGCGCTTTTTCTTTGAGAATGTTGATATCTGCGTTCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATATGGCAGCCGAAAGCCATGTTGACCGTTCTATTg                           <  1:1772135/150‑1 (MQ=255)
                                                                                                                          gTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCGTTGAGAAAGTTGATATCTGCGGTCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGCCCGGTCTATTGa                          >  2:812902/1‑150 (MQ=255)
                                                                                                                            cGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACg                        >  2:658074/1‑150 (MQ=255)
                                                                                                                            cGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACg                        >  1:72937/1‑150 (MQ=255)
                                                                                                                              cTGCCACCGTCCGCGCCAAGCGAGCCCCTCGCCTGTGACCAAGAAGAACTCTGCGATCGGGCAGAACTGGCAGGCGTATTCACAGAGCATCACCCAGACTGTGTCATGGATCTGGCAGCCGACAGCCATGTTGACCGTTCTATTGACGGc                      <  2:307353/150‑1 (MQ=255)
                                                                                                                                  cACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCCCGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCgg                  >  2:521544/1‑150 (MQ=255)
                                                                                                                                    ccGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGCGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCa                >  1:1369262/1‑150 (MQ=255)
                                                                                                                                       caaaggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCg             >  2:1861977/6‑150 (MQ=255)
                                                                                                                                        tCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACGGGCACGCGGATTCAGTGAGCATCTGCCTGACTGTGTCATGCTTCAAGCAGACGAAGGACATGTTGCACGTGCTATTGACGGCCCGGCAGCGt            >  1:358778/1‑150 (MQ=255)
                                                                                                                                         cgcgCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGtt           >  2:241728/1‑150 (MQ=255)
                                                                                                                                         aaggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCCCTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGtt           >  1:1559450/4‑150 (MQ=255)
                                                                                                                                          gcgcAAAGCGAGCGCTTTGCCTTTGAGAAAGATGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTAATGCATCTGCCAGCCGAAAGCCATGTTGACCGTTCTATTTACGGCCCGGCAgaggtt          >  1:361679/1‑145 (MQ=255)
                                                                                                                                           cgcAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTa         >  2:600057/1‑150 (MQ=255)
                                                                                                                                           cgcAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGACATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTa         <  2:1638512/150‑1 (MQ=255)
                                                                                                                                              aaaGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTg      <  2:48338/150‑1 (MQ=255)
                                                                                                                                                 gcgagcgCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGaaa   <  1:1008144/150‑1 (MQ=255)
                                                                                                                                                  cgagcgCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGAAAc  <  1:501608/150‑1 (MQ=255)
                                                                                                                                                  cgagcgCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGAAAc  <  1:362107/150‑1 (MQ=255)
                                                                                                                                                  cgagcgCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGAAAc  <  2:1179775/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGAAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGAAAC  >  NC_000913/3972525‑3972820

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: