Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,479,020 A→T 17.3% intergenic (+131/+14) rraB → / ← yjgM protein inhibitor of RNase E/GNAT family putative N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,479,0200AT17.3% 252.4 / 42.6 103intergenic (+131/+14)rraB/yjgMprotein inhibitor of RNase E/GNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (48/37);  new base T (7/11);  total (55/48)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

Only 100 of 111 total aligned reads displayed.

CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCT  >  NC_000913/4478921‑4479119
                                                                                                     |                                                                                                   
cggCAACCTGGTGAATTTGTAGCCCGGATAAAGCATTCATGTGACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAaca                                                                                                     >  1:1859768/3‑100 (MQ=255)
cTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAaca                                                                                                     <  1:718902/100‑1 (MQ=255)
 tGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGAATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGGTGAacac                                                                                                    >  1:841948/1‑100 (MQ=255)
    aaCCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACAcgcg                                                                                                 <  2:2645795/100‑1 (MQ=255)
    aaCCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACAcgcg                                                                                                 <  2:1968143/100‑1 (MQ=255)
     aCCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCgg                                                                                                >  2:2035228/1‑100 (MQ=255)
      ccTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCggg                                                                                               <  1:1543606/100‑1 (MQ=39)
      ccTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATACCCGCGAGTTCAACAcacg                                                                                                 >  1:988975/1‑98 (MQ=255)
       cTGGTGAATTTGTAGCCCGTATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTCTTAATGTTCGAATGCCCGCGTGTTCAACTCGCggg                                                                                               <  2:1670630/99‑1 (MQ=38)
       cTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGc                                                                                              <  1:1434682/100‑1 (MQ=39)
       cTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGc                                                                                              <  1:261301/100‑1 (MQ=39)
       cTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGc                                                                                              <  1:992098/100‑1 (MQ=39)
       cTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGc                                                                                              >  1:715053/1‑100 (MQ=255)
         ggTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCAt                                                                                            <  2:299367/100‑1 (MQ=38)
         ggTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAt                                                                                            >  2:154675/1‑100 (MQ=255)
         ggTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAt                                                                                            >  2:513350/1‑100 (MQ=255)
         ggTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAt                                                                                            >  1:780655/1‑100 (MQ=255)
         ggTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAt                                                                                            >  1:1897854/1‑100 (MQ=255)
           tGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttt                                                                                          <  2:137350/100‑1 (MQ=255)
           tGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttt                                                                                          >  2:2675508/1‑100 (MQ=255)
           tGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttt                                                                                          >  1:2745619/1‑100 (MQ=255)
           tGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TGGTTAATGTTCGAATGCCCGCGAGTTGAACACGCggg                                                                                               <  2:1731136/95‑1 (MQ=255)
            gAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtttt                                                                                         <  1:617131/100‑1 (MQ=255)
            gAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtttt                                                                                         >  1:1900993/1‑100 (MQ=255)
             aaTTTGTATCCCTGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATTTTCGAATGCCCGCGAGTTTAACACGCGGGCAtttt                                                                                         <  1:117314/99‑1 (MQ=39)
              aTTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtttttt                                                                                       >  1:178836/1‑100 (MQ=255)
              aTTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGGACG‑‑TTGTTAATGTTCAAATGCCCGCGTGTTCAACTCGCGGGCAttcgaa                                                                                       <  2:1174720/100‑5 (MQ=37)
               tttGTAGCCCGAATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCAttcgaac                                                                                      <  1:924824/100‑6 (MQ=37)
                ttGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCAttcgaaca                                                                                     <  1:1207869/100‑7 (MQ=37)
                ttGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTa                                                                                     >  1:1021766/1‑100 (MQ=255)
                 tGTAGCCCGGATAAGGCATTCATGTCACATCCGGCACAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTag                                                                                    <  2:584757/100‑1 (MQ=255)
                  gTAGCCCGGATAAGTCATTGATGTCACATCGGGCATAAACAAAATGTGCG‑‑TTGTTAATTTTGGAATGCCCGCGAGTTGAACACTCGGGCAtttgttt                                                                                      <  2:2734926/97‑1 (MQ=255)
                  gTAACCCGAATAGGGCAGTCCCGTCACATCCAGCATAGACAAAACGTACG‑‑TTGTTAATTTTCGAGTGCCCGCGTGTTCTCCTCGCGGGAAtt                                                                                           <  2:1647846/92‑1 (MQ=255)
                   tAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTagag                                                                                  <  1:1565850/100‑1 (MQ=255)
                   tAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTagag                                                                                  <  2:1055657/100‑1 (MQ=255)
                    aGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGt                                                                                 >  1:1522798/1‑100 (MQ=255)
                     acccGGATAAGGACTTCATCTTCCATCCGACATAAACAAAACGTACG‑‑TTGGTGATGGGCGAATGCCCGGGAGGTGAACACTCGGGCGGGTTTTTAgtgtg                                                                                >  2:1882348/2‑99 (MQ=255)
                         ggATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGGGTTcgcg                                                                            >  1:872731/1‑100 (MQ=255)
                              aGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACa                                                                       >  1:1914206/1‑100 (MQ=255)
                              aGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACa                                                                       >  2:2062385/1‑100 (MQ=255)
                              aGGCATTCATGTCACATCCGACATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACa                                                                       <  2:2112743/100‑1 (MQ=255)
                               ggCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACAt                                                                      >  1:34263/1‑100 (MQ=255)
                               ggCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACAt                                                                      >  2:1569331/1‑100 (MQ=255)
                                gCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATc                                                                     >  1:574003/1‑100 (MQ=255)
                                gCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATc                                                                     >  1:2548379/1‑100 (MQ=255)
                                    tCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTTATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCa                                                                 <  2:1231591/100‑1 (MQ=255)
                                     cATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCAc                                                                <  2:1302346/100‑1 (MQ=255)
                                     cATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCAc                                                                >  1:1689641/1‑100 (MQ=255)
                                        gTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTc                                                             <  2:1287718/100‑1 (MQ=255)
                                        gTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTc                                                             <  1:1996176/100‑1 (MQ=255)
                                           acaTCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTcaca                                                          >  1:1477901/1‑100 (MQ=255)
                                            caTCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACaa                                                         >  1:1344258/1‑100 (MQ=255)
                                              tCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATc                                                       <  1:732874/100‑1 (MQ=255)
                                               ccGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCg                                                      <  1:976845/100‑1 (MQ=255)
                                                 ggCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGAc                                                    >  1:448117/1‑100 (MQ=255)
                                                   cATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACAt                                                  >  1:701046/1‑100 (MQ=255)
                                                   aacattaacaaCGTACGTTTTGTTTATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGAc                                                    >  2:1264098/10‑100 (MQ=37)
                                                       aaCAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGcc                                              <  2:2367180/100‑1 (MQ=255)
                                                        aCAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGccc                                             >  1:1339603/1‑100 (MQ=255)
                                                        aCAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGccc                                             >  1:1973567/1‑100 (MQ=255)
                                                        aCAAAACGTACG‑‑TTGTAAATGTTCGAATGCCCGCGAGATGAACACGCGGGCATTTTTTTAGAGTTCGCGCAAAATCCGCACTTCACAATCGACATGGccc                                             >  2:319679/1‑100 (MQ=255)
                                                         cAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCg                                            >  2:1868222/1‑100 (MQ=255)
                                                                      ttgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAgc                                 <  1:2143417/100‑1 (MQ=255)
                                                                       tgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAgcg                                <  1:370048/100‑1 (MQ=255)
                                                                           aaTGTTCGACTGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGCTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATa                            >  2:1244193/1‑100 (MQ=255)
                                                                           aaTGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATa                            <  1:1502098/100‑1 (MQ=255)
                                                                            ctaaaaaaaTGCCCGCGTGTTCAACTCGAGCGCATCTTTTTAGAGTTCGCGCACCAAACGCACTGCACACTCGACATGGCCCGTGACGCCCACCGCAtca                           >  2:2003610/8‑98 (MQ=12)
                                                                              gTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTc                         >  1:935910/1‑100 (MQ=255)
                                                                              gTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTc                         >  1:742179/1‑100 (MQ=255)
                                                                              aaaaaaaTGCCCGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTc                         >  1:1783826/6‑100 (MQ=37)
                                                                               ttCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACCTGGCCCGTGCAGCAAAGAGCATAGTCg                        >  2:2135658/1‑100 (MQ=255)
                                                                                tCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGa                       >  1:1924758/1‑100 (MQ=255)
                                                                                  gAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCTCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGata                     <  2:1522791/100‑1 (MQ=255)
                                                                                  gAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGata                     <  1:1358872/100‑1 (MQ=255)
                                                                                  gAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGata                     <  2:238271/100‑1 (MQ=255)
                                                                                  aaaTGCCCGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGata                     >  1:1250698/2‑100 (MQ=38)
                                                                                    aTGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGatat                    <  2:324406/99‑1 (MQ=255)
                                                                                      gCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCt                 <  2:348017/100‑1 (MQ=255)
                                                                                      gCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCt                 <  2:854756/100‑1 (MQ=255)
                                                                                       cccGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTc                >  1:2628064/1‑100 (MQ=255)
                                                                                       cccGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTc                >  2:2729745/1‑100 (MQ=255)
                                                                                       cccGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTc                >  2:1111020/1‑100 (MQ=255)
                                                                                       cccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTc                >  1:1110025/1‑100 (MQ=255)
                                                                                       cccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTc                >  2:2619651/1‑100 (MQ=255)
                                                                                        ccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCa               >  1:759023/1‑100 (MQ=255)
                                                                                         cgcgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCaa              >  2:1411433/1‑100 (MQ=255)
                                                                                         cgcgAGTTGAACACGCGGGCATTTTTTTAGAGCTCGCGCAACCTCCGCACTTCACAATCGACAGGGCCCGTGCAGCCAACCGCCTACTCGATACGCTCaa              >  2:1208883/1‑100 (MQ=255)
                                                                                           cgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAg            <  2:1021762/100‑1 (MQ=255)
                                                                                           cgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAg            <  1:2062398/100‑1 (MQ=255)
                                                                                            gAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGACCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGc           <  2:2570002/100‑1 (MQ=255)
                                                                                             tgTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGcc          >  1:1158235/2‑100 (MQ=255)
                                                                                                tGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCaa       >  2:2193362/1‑100 (MQ=255)
                                                                                                tGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCa        <  1:200886/99‑1 (MQ=255)
                                                                                                 gAACACGCGGGCATTTTTTTAGAGATGGCGCAACATCCGCCCTTCACAATCGACATGGCTCGTGCAGCCAAGCGCATAGTCGATATGTTCaaa             >  2:1684088/1‑93 (MQ=255)
                                                                                                  aaCACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAAt     >  1:615648/1‑100 (MQ=255)
                                                                                                    cacGCTGGCATTTTTTTAGACTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGc   <  2:2095790/100‑1 (MQ=255)
                                                                                                    cacGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGc   <  1:2381592/100‑1 (MQ=255)
                                                                                                    cacGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGc   >  2:590026/1‑100 (MQ=255)
                                                                                                    cacGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGc   >  2:878848/1‑100 (MQ=255)
                                                                                                     acGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCt  >  2:1323230/1‑100 (MQ=255)
                                                                                                     |                                                                                                   
CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCT  >  NC_000913/4478921‑4479119

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