Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,177,172 T→A 15.9% intergenic (+44/‑186) tufB → / → secE translation elongation factor EF‑Tu 2/preprotein translocase membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,177,1720TA15.9% 284.2 / 41.7 114intergenic (+44/‑186)tufB/secEtranslation elongation factor EF‑Tu 2/preprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (53/42);  new base A (7/11);  total (61/53)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

Only 100 of 117 total aligned reads displayed.

TCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATA  >  NC_000913/4177073‑4177269
                                                                                                   |                                                                                                 
 ccGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAttt                                                                                                  <  2:1828266/100‑1 (MQ=255)
  cGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAtt                                                                                                   <  2:1649347/98‑1 (MQ=255)
  cGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTg                                                                                                 >  2:2047464/1‑100 (MQ=255)
  cGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTg                                                                                                 <  1:2094084/100‑1 (MQ=255)
  cGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTg                                                                                                 >  1:1598729/1‑100 (MQ=255)
    tGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGAt                                                                                               <  2:157857/100‑1 (MQ=255)
    tGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGAt                                                                                               >  2:1269398/1‑100 (MQ=255)
    tGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGAt                                                                                               <  1:1000653/100‑1 (MQ=255)
    tGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGAt                                                                                               <  1:1886818/100‑1 (MQ=255)
         gcggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCt                                                                                          <  1:1279769/100‑1 (MQ=255)
         gcggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCt                                                                                          <  1:1237303/100‑1 (MQ=255)
         gcggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCt                                                                                          <  2:913238/100‑1 (MQ=255)
          cggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtt                                                                                         <  2:434025/100‑1 (MQ=255)
           ggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCttt                                                                                        <  2:643680/100‑1 (MQ=255)
           ggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttt                                                                                        <  2:931656/100‑1 (MQ=255)
            gcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCtttt                                                                                       <  1:2036536/100‑1 (MQ=255)
            gcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtttt                                                                                       <  2:244076/100‑1 (MQ=255)
            gcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtt                                                                                         <  2:127267/98‑1 (MQ=255)
              cGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTg                                                                                     <  1:656918/100‑1 (MQ=255)
              cGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttttag                                                                                     <  2:2124382/100‑3 (MQ=255)
               gTACCGTTGGGGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGTCGATAACATTTGACTCAATGCGCACTAAAAGGGCATCAAATGATGCCCttttag                                                                                     <  2:226196/99‑3 (MQ=38)
                tACCGTTGGCGGGGGCGGTGTGGCAAAAGTTCTGAGCTAATTGCGGATAAGATTTGACGCCAGGGGCACTATCATGTCATCAAATGATGCCCtatt                                                                                       <  2:475268/96‑1 (MQ=255)
                tACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCa                                                                                   >  2:2261253/1‑100 (MQ=255)
                tACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCa                                                                                   <  1:1907811/100‑1 (MQ=255)
                tACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCa                                                                                   >  2:2022767/1‑100 (MQ=255)
                tACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttttagtg                                                                                   <  1:2265906/100‑5 (MQ=38)
                tACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACAAAAAGGGCATCAAATGATGCCCttttagtg                                                                                   <  2:437295/100‑5 (MQ=38)
                atctgTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCa                                                                                   <  2:2131647/96‑1 (MQ=255)
                 aCCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCAc                                                                                  <  1:1210895/100‑1 (MQ=255)
                  ccGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACg                                                                                 <  2:647803/100‑1 (MQ=255)
                  ccGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACg                                                                                 >  1:2063477/1‑100 (MQ=255)
                    gTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCaagat                                                                               >  1:2131671/1‑96 (MQ=255)
                    gTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACg                                                                                 >  1:1990734/1‑98 (MQ=255)
                       ggCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTAACGCAATGCGCACTAAAAGGGCAACATTTGATGCCCTTTTTGCACGCttt                                                                             >  2:234939/1‑99 (MQ=255)
                                  tGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAAcc                                                                 <  1:214083/100‑1 (MQ=255)
                                    tAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTg                                                               >  2:1817917/1‑100 (MQ=255)
                                    tAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTg                                                               <  2:1831223/100‑1 (MQ=255)
                                     aGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTgg                                                              >  1:954345/1‑100 (MQ=255)
                                      gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGc                                                             <  2:1993120/100‑1 (MQ=255)
                                      gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGc                                                             >  2:1195214/1‑100 (MQ=255)
                                       cAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCt                                                            >  1:1797692/1‑100 (MQ=255)
                                       cAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCt                                                            <  1:2396273/100‑1 (MQ=255)
                                         aaaGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtca                                                          <  2:890554/100‑1 (MQ=255)
                                           aGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtcatc                                                        <  1:2022789/100‑1 (MQ=255)
                                           aGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtcatc                                                        <  1:1156947/100‑1 (MQ=255)
                                           aGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtcagc                                                        <  2:336390/100‑3 (MQ=255)
                                            gTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtcatca                                                       >  1:403517/1‑100 (MQ=255)
                                             ttCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAg                                                      <  2:810869/100‑1 (MQ=255)
                                                tGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGa                                                   >  1:2037636/1‑100 (MQ=255)
                                                tGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGa                                                   >  2:817788/1‑100 (MQ=255)
                                                 gAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGAt                                                  >  2:2137602/1‑100 (MQ=255)
                                                  aGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGAtt                                                 >  2:1758279/1‑100 (MQ=255)
                                                   gCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGAttt                                                >  2:817590/1‑100 (MQ=255)
                                                     tAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTc                                              >  1:1350052/1‑100 (MQ=255)
                                                      aaTTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCt                                             >  1:1893088/1‑100 (MQ=255)
                                                         tGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTg                                          <  2:1776461/100‑1 (MQ=255)
                                                           ccGATAACATTTGACGCAATGCGCACTAAAAGGGCATAATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTc                                        >  2:1908820/1‑100 (MQ=255)
                                                            cGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCa                                       >  1:1603390/1‑100 (MQ=255)
                                                            cGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCa                                       >  2:109480/1‑100 (MQ=255)
                                                             gATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCAt                                      >  1:1920302/1‑100 (MQ=255)
                                                                aaCATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAAt                                   <  1:2022238/100‑1 (MQ=255)
                                                                 aCATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATc                                  >  2:1276495/1‑100 (MQ=255)
                                                                 aCATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATc                                  >  1:2400674/1‑100 (MQ=255)
                                                                  cATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCa                                 >  2:2147068/1‑100 (MQ=255)
                                                                   aTTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCAt                                >  2:2216307/1‑100 (MQ=255)
                                                                    tttGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCAtt                               >  1:492452/1‑100 (MQ=255)
                                                                      tGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGc                             >  2:737660/1‑100 (MQ=255)
                                                                      tGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTGTGCATGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGc                             <  2:822839/100‑1 (MQ=255)
                                                                       gACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTTTCATAATCATTGc                             >  1:657699/1‑99 (MQ=255)
                                                                       gACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCt                            >  1:686913/1‑100 (MQ=255)
                                                                       gACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCt                            >  2:1785941/1‑100 (MQ=255)
                                                                        aCGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTg                           <  1:1344889/100‑1 (MQ=255)
                                                                         cGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTga                          <  1:2098838/100‑1 (MQ=255)
                                                                         cGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTga                          >  1:1050801/1‑100 (MQ=255)
                                                                         cGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATAATTGCTga                          >  1:1861369/1‑100 (MQ=255)
                                                                           cAATGCGCCCTACAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTTCTgaga                        >  2:521397/1‑100 (MQ=255)
                                                                           cAATGCGCACTAAAAGGGCAACATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTgaga                        >  2:1464951/1‑100 (MQ=255)
                                                                            aaTGCGCACTAAAAGGGAATCAAATGATGCCCTTTTTGCACGATTTGGTACCAGAACCTAGCTCATCCGTGATTTTCTTTGTCATAATAATTGCTgag                         >  2:822239/1‑98 (MQ=255)
                                                                              tGCGCACTAAAAGGGCATCATTTGATGCCCTCTTTTCACGCTTTCGTACCAGAACCTGGTTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAg                     <  2:1976279/100‑1 (MQ=255)
                                                                               gcgcACTAAAGTGGCACCATGTGATGCCCTTTTTGCACACTTTCGTACCAGAACGTGGCTCATCCTTGGTGTTCTATCGCACACTCATTGCt                            >  2:1468535/1‑92 (MQ=255)
                                                                                cgtgcaaaaaGGGCATCAAATGATGCCCTATTAGTGCGCTTTCGTACCACAACCTGGCTCATCAGTGATGTTCTTTGTCATAATAATTGCTGAGACAGGc                   >  2:962991/7‑100 (MQ=21)
                                                                                cgcACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAg                     >  1:621021/1‑98 (MQ=255)
                                                                                cgcACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGc                   >  1:1557427/1‑100 (MQ=255)
                                                                                cgcACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGc                   <  2:403514/100‑1 (MQ=255)
                                                                                 gcACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGct                  >  1:1522756/1‑100 (MQ=255)
                                                                                  tgcaaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGctc                 >  2:1892022/5‑100 (MQ=38)
                                                                                    caaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGCTTTTCTTTGTCATAATCATTGCTGAGCCAGGc                   >  2:1868114/3‑96 (MQ=38)
                                                                                    caaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTg               >  1:643397/3‑100 (MQ=255)
                                                                                     tAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGt              >  2:624888/1‑100 (MQ=255)
                                                                                        aaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGa           >  2:1882809/1‑100 (MQ=255)
                                                                                          gggCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAgg         <  2:1522859/100‑1 (MQ=255)
                                                                                          gggCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAgg         <  1:2070416/100‑1 (MQ=255)
                                                                                          gggCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAgg         >  2:628907/1‑100 (MQ=255)
                                                                                           ggCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAggg        <  2:1855136/100‑1 (MQ=255)
                                                                                            gCATCATTTGATGCCCTTTTTGCACGCATTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGc       >  2:369998/1‑100 (MQ=255)
                                                                                            gCATCATTTGATGCCCTTTTTGCACACTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTCTCATAATCATTGCTCACACCCGCTCTGTTAAgcg        >  2:904528/1‑97 (MQ=255)
                                                                                             catcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCg      <  1:2327024/100‑1 (MQ=255)
                                                                                             catcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCg      >  1:1760790/1‑100 (MQ=255)
                                                                                                catTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGCTTTTCTTTGTCATAATCATTGCTGATACAGGCTCTGTTGAGGGCGtat   <  2:560140/100‑1 (MQ=255)
                                                                                                 atTTGATGCCCTTTTTGCACGCTTTCATACAAGAACCTGGCTCATCAGTAATTTTCTTTGTCATAATCATTGCTGAGACAGGCACGGTAGAGGGCGtata  >  2:1286182/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
TCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATA  >  NC_000913/4177073‑4177269

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