Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,201,035 T→A 5.7% F253I (TTT→ATT)  yehI → DUF4132 domain‑containing protein YehI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,201,0350TA 5.7% 223.7 / 3.2 70F253I (TTT→ATT) yehIDUF4132 domain‑containing protein YehI
Reads supporting (aligned to +/- strand):  ref base T (38/28);  new base A (2/2);  total (40/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

CTGATCCTCCCGGAAAAACCAGAAATAGCCAATGAGTTGGTAGGCCTTGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGCCCATAAAGAAGATTGCGGCAGTCTGCTGGTACAA  >  NC_000913/2200885‑2201184
                                                                                                                                                      |                                                                                                                                                     
cTGATCCTCCGGAAAAACCAGAAATAGCCAATGAGTTGGTAGGCCTTGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGatatat                                                                                                                                                       >  2:1416173/1‑150 (MQ=255)
  gATCCTCCCGGAAAAACCAGAAATAGCCAATGAGTTGGTAGGCCTTGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATAtt                                                                                                                                                      <  2:976905/150‑1 (MQ=255)
      ctccCGGAAAAACCAGAAATAGCCAATGAGTTGGTAGGCCTTGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGc                                                                                                                                                  <  2:546951/150‑1 (MQ=255)
         ccGGAAAAACCAGAAATAGCCAATGAGTTGGTAGGCCTTGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGAt                                                                                                                                               <  2:910038/150‑1 (MQ=255)
          cGGAAAAACCAGAAATAGCCAATGAGTTGGTAGGCCTTGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATc                                                                                                                                              >  1:1467577/1‑150 (MQ=255)
           ggAAAAACCAGAAATAGCCAATGAGTTGGTAGGCCTTGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCg                                                                                                                                             <  1:1004351/150‑1 (MQ=255)
              aaaaCCAGAAATAGCCAATGAGTTGGTAGGCCTTGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCgaga                                                                                                                                          <  1:463453/150‑1 (MQ=255)
               aaaCCAGAAATAGCCAATGAGTTGGTAGGCCTTCAATGCCCGCGAACGCATTTTCATTCTAAGGAGGGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGAGCGAGaa                                                                                                                                         <  2:1149041/150‑1 (MQ=255)
                aaCCAGAAATAGCCAATGAGTTGGTAGGCCTTGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAg                                                                                                                                        <  2:1349335/150‑1 (MQ=255)
                  ccAGAAATAGCCAATGAGTTGGTAGGCCTTGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGc                                                                                                                                                  >  2:509079/1‑138 (MQ=255)
                  ccAGAAATAGCCAATGAGTTGGTAGGCCTTGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGc                                                                                                                                                  <  1:509079/138‑1 (MQ=255)
                  ccAGAAATAGCCAATGAGTTGGTAGGCCTTGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGcc                                                                                                                                      >  2:759650/1‑150 (MQ=255)
                           gCCAATGAGTTGGTAGGCCTTGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATg                                                                                                                             <  1:123263/150‑1 (MQ=255)
                             cAATGAGTTGGTAGGCCTTGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTc                                                                                                                           >  2:998591/1‑150 (MQ=255)
                             cAATGAGTTGGTAGGCCTTGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTc                                                                                                                           >  2:202748/1‑150 (MQ=255)
                              aaTGAGTTGGTAGGCCTTGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCg                                                                                                                          >  1:1492388/1‑150 (MQ=255)
                               aTGAGTTGGTAGGCCTTGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGc                                                                                                                         <  2:1244058/150‑1 (MQ=255)
                                    ttGGTAGGCCTTGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGaa                                                                                                                    <  2:422469/150‑1 (MQ=255)
                                      ggTAGGCCTTGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGaaaa                                                                                                                  >  2:304126/1‑150 (MQ=255)
                                      ggTAGGCCTTGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGaaaa                                                                                                                  >  2:518638/1‑150 (MQ=255)
                                      ggTAGGCCTTGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATAATTAGCGATCGAGAAGCCAGCTACATGTCGCATGaaaa                                                                                                                  >  1:170523/1‑150 (MQ=255)
                                         aGGCCTTGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCa                                                                                                               >  1:79772/1‑150 (MQ=255)
                                          ggCCTTGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCAc                                                                                                              >  2:853722/1‑150 (MQ=255)
                                           gCCTTGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACt                                                                                                             >  1:922312/1‑150 (MQ=255)
                                               ttaaTGCCCTGGAACTCATTTTCATTCTAAGGAGTGGGTAATAGTTGTTGCTTATGACCCCACAGCGGTGAGAAAACTCAAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCgg                                                                                                         <  2:513378/148‑1 (MQ=255)
                                               tGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGCACCCACAGCCGTCAGAAAACTCGAACACTACTGGAGCCAGGATATATATAGCGATCGAGAAGCCAGCTACAAGTCGCATGAAAACCACTTCgg                                                                                                         >  2:845553/1‑150 (MQ=255)
                                                gAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGc                                                                                                        <  1:1381711/150‑1 (MQ=255)
                                                   tGCCCGAGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTAc                                                                                                     >  2:186920/1‑150 (MQ=255)
                                                     cccGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTAcgc                                                                                                   >  2:773020/1‑150 (MQ=255)
                                                                 ttttCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGcgccgc                                                                                       >  2:215068/1‑150 (MQ=255)
                                                                        tCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGTATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTg                                                                                <  1:459316/150‑1 (MQ=255)
                                                                        tCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTg                                                                                <  2:79772/150‑1 (MQ=255)
                                                                        tCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTg                                                                                <  2:196086/150‑1 (MQ=255)
                                                                         cTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTgc                                                                               >  2:124275/1‑150 (MQ=255)
                                                                         cTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTgc                                                                               >  2:1461884/1‑150 (MQ=255)
                                                                             ggAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGa                                                                           <  2:1492388/150‑1 (MQ=255)
                                                                                  ggTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAg                                                                      >  2:788367/1‑150 (MQ=255)
                                                                                   gTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAgg                                                                     <  1:1662406/150‑1 (MQ=255)
                                                                                   gTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAgg                                                                     >  1:1567273/1‑150 (MQ=255)
                                                                                     tAAAAGTTGTTGCTAATGACACCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGAc                                                                   >  2:403073/1‑150 (MQ=255)
                                                                                       aaaGTTGTTGCTAATGACCCCCCAGCGGTGAGAACACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTg                                                                 <  2:1188196/150‑1 (MQ=255)
                                                                                             gttgCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGcc                                                           <  2:574913/150‑1 (MQ=255)
                                                                                             gttgCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGcc                                                           >  1:1651690/1‑150 (MQ=255)
                                                                                               tgCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCAt                                                         >  2:239658/1‑150 (MQ=255)
                                                                                               tgCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCAt                                                         >  1:909657/1‑150 (MQ=255)
                                                                                                gCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCAtt                                                        >  2:1069648/1‑150 (MQ=255)
                                                                                                   aaTGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCcg                                                     >  2:288904/1‑150 (MQ=255)
                                                                                                   aaTGACCCCACAGCGGTGAGAAAACTAGAACACTACTGGAGCCAGCATATATTTACCGATCGAGAAGCCAGATACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTccc                                                     >  2:572722/1‑149 (MQ=255)
                                                                                                      gACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCgcgc                                                  >  1:1490022/1‑150 (MQ=255)
                                                                                                       aCCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGATAACCACTTCGGCTACGCGGCCTGCGCCGCGCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGcc                                                 <  1:124275/150‑1 (MQ=255)
                                                                                                       aCCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGcc                                                 >  1:386610/1‑150 (MQ=255)
                                                                                                        ccccACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCt                                                >  1:1168318/1‑150 (MQ=255)
                                                                                                           cacaGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTcgc                                             <  1:773020/150‑1 (MQ=255)
                                                                                                           cacaGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATAATTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTcgc                                             >  2:225701/1‑150 (MQ=255)
                                                                                                                 ggTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTa                                       <  1:1416173/150‑1 (MQ=255)
                                                                                                                  gTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTAt                                      <  2:922312/150‑1 (MQ=255)
                                                                                                                      gaAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGccc                                  >  2:40846/1‑150 (MQ=255)
                                                                                                                      gaAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGccc                                  <  1:202748/150‑1 (MQ=255)
                                                                                                                      gaAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGccc                                  >  2:34734/1‑150 (MQ=255)
                                                                                                                      gaAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGccc                                  >  1:452234/1‑150 (MQ=255)
                                                                                                                       aaaaCTGGAACACTACTTGAGCCAGGATATAATTATCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGCCCa                                 <  2:170523/150‑1 (MQ=255)
                                                                                                                           cTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGCCCATaaa                             >  1:1277824/1‑150 (MQ=255)
                                                                                                                            tCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGCCCATAaag                            <  2:1595875/150‑1 (MQ=255)
                                                                                                                              gAACACTACTGGAGCCAGGATATAATTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGCCCATAaagaa                          <  1:225701/150‑1 (MQ=255)
                                                                                                                                   ctactGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGCCCATAAAGAAGATTg                     <  1:304126/150‑1 (MQ=255)
                                                                                                                                        ggAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTAGGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGCCCATAAAGAAGATTGCGGCa                >  2:1418117/1‑150 (MQ=255)
                                                                                                                                        ggAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGCCCATAAAGAAGATTGCGGCa                >  1:876688/1‑150 (MQ=255)
                                                                                                                                        ggAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGCCCATAAAGAAGATTGCGGCa                >  1:685584/1‑150 (MQ=255)
                                                                                                                                         gAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGCCCATAAAGAAGATTGCGGCAg               <  2:1467577/150‑1 (MQ=255)
                                                                                                                                         gAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGCCCATAAAGAAGATTGCGGCAg               <  2:876688/150‑1 (MQ=255)
                                                                                                                                           gCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGCCCATAAAGAAGATTGCGGCAGTc             >  1:840639/1‑150 (MQ=255)
                                                                                                                                            ccAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGCCCATAAAGAAGATTGCGGCAGTct            >  1:1226713/1‑150 (MQ=255)
                                                                                                                                                 atatatTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGCCCATAAAGAAGATTGCGGCAGTCTGCTgg       <  2:1173309/150‑1 (MQ=255)
                                                                                                                                                 atatatTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGCCCATAAAGAAGATTGCGGCAGTCTGCTgg       <  1:239658/150‑1 (MQ=255)
                                                                                                                                                  tatatTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCTATGTATGCCCATAAAGAATATTGCGGCAGTCTGCTGGt      <  2:909657/150‑1 (MQ=255)
                                                                                                                                                  tatatTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGCCCATAAAGAAGATTGCGGCAGTCTGCTGGt      >  1:1391559/1‑150 (MQ=255)
                                                                                                                                                  aatatTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGCCCATAAAGAAGATTGCGGCAGTCTGCTGGt      >  1:444569/2‑150 (MQ=255)
                                                                                                                                                    tatTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGCCCATAAAGAAGATTGCGGCAGTCTGCTGGTAc    >  2:1544776/1‑150 (MQ=255)
                                                                                                                                                     atTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGCCCATAAAGAAGATTGCGGCAGTCTGCTGGTACa   >  2:1629815/1‑150 (MQ=255)
                                                                                                                                                      tttAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGCCCATAAAGAAGATTGCGGCAGTCTGCTGGTACaa  <  2:1490022/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CTGATCCTCCCGGAAAAACCAGAAATAGCCAATGAGTTGGTAGGCCTTGAATGCCCGCGAACTCATTTTCATTCTAAGGAGTGGTTAAAAGTTGTTGCTAATGACCCCACAGCGGTGAGAAAACTCGAACACTACTGGAGCCAGGATATATTTAGCGATCGAGAAGCCAGCTACATGTCGCATGAAAACCACTTCGGCTACGCGGCCTGCGCCGCCCTTTTGCGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGCCCATAAAGAAGATTGCGGCAGTCTGCTGGTACAA  >  NC_000913/2200885‑2201184

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