Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,479,012 A→T 15.9% intergenic (+123/+22) 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,0120AT15.9% 279.4 / 41.9 107intergenic (+123/+22)rraB/yjgMprotein inhibitor of RNase E/GNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (50/40);  new base T (6/11);  total (56/51)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

Only 100 of 107 total aligned reads displayed.

AGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGC  >  NC_000913/4478916‑4479110
                                                                                                  |                                                                                                  
 gTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCTGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTg                                                                                                <  1:2294655/100‑1 (MQ=255)
 gTTCCTGCAACCTGGTGAATTTGTAGCCAGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTg                                                                                                >  1:44100/1‑100 (MQ=255)
  ttCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGa                                                                                               <  1:1530213/100‑1 (MQ=255)
    ccTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAac                                                                                             <  1:1729999/100‑1 (MQ=255)
    ccTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAac                                                                                             >  1:229094/1‑100 (MQ=255)
     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)
              ggTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCAt                                                                                   <  2:299367/100‑1 (MQ=38)
              ggTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAt                                                                                   >  2:154675/1‑100 (MQ=255)
              ggTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAt                                                                                   >  1:780655/1‑100 (MQ=255)
              ggTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAt                                                                                   >  2:513350/1‑100 (MQ=255)
              ggTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAt                                                                                   >  1:1897854/1‑100 (MQ=255)
                tGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttt                                                                                 >  1:2745619/1‑100 (MQ=255)
                tGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttt                                                                                 <  2:137350/100‑1 (MQ=255)
                tGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttt                                                                                 >  2:2675508/1‑100 (MQ=255)
                tGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TGGTTAATGTTCGAATGCCCGCGAGTTGAACACGCggg                                                                                      <  2:1731136/95‑1 (MQ=255)
                 gAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtttt                                                                                >  1:1900993/1‑100 (MQ=255)
                 gAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtttt                                                                                <  1:617131/100‑1 (MQ=255)
                  aaTTTGTATCCCTGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATTTTCGAATGCCCGCGAGTTTAACACGCGGGCAtttt                                                                                <  1:117314/99‑1 (MQ=39)
                   aTTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtttttt                                                                              >  1:178836/1‑100 (MQ=255)
                   aTTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtttttt                                                                              >  1:1449741/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)
                        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                                                             >  2:1569331/1‑100 (MQ=255)
                                    ggCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACAt                                                             >  1:34263/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)
                                         tCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCa                                                        <  1:831777/100‑1 (MQ=255)
                                          cATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCAc                                                       >  1:1689641/1‑100 (MQ=255)
                                          cATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCAc                                                       <  2:1302346/100‑1 (MQ=255)
                                             gTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTc                                                    <  1:1996176/100‑1 (MQ=255)
                                             gTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTc                                                    <  2:1287718/100‑1 (MQ=255)
                                               cacaTCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTcac                                                  >  2:2181253/1‑100 (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)
                                                                     tacg‑‑tTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGcc                            <  1:598913/100‑1 (MQ=255)
                                                                     tacg‑‑tTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGc                             >  2:1331617/1‑99 (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:742179/1‑100 (MQ=255)
                                                                                   gTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTc                >  1:935910/1‑100 (MQ=255)
                                                                                   aaaaaaaTGCCCGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTc                >  1:1783826/6‑100 (MQ=37)
                                                                                    ttCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACCTGGCCCGTGCAGCAAAGAGCATAGTCg               >  2:2135658/1‑100 (MQ=255)
                                                                                    ttCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCg               >  1:1867302/1‑100 (MQ=255)
                                                                                     tCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGa              >  1:1924758/1‑100 (MQ=255)
                                                                                       gAATGCCCGCGAGTTGAACACGTGGGCATTTCTTTAGTGTTCGCGCATCATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGata            <  2:1859760/100‑1 (MQ=255)
                                                                                       gAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCTCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGata            <  2:1522791/100‑1 (MQ=255)
                                                                                       gAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGata            <  1:1358872/100‑1 (MQ=255)
                                                                                       gAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGata            <  2:238271/100‑1 (MQ=255)
                                                                                       gAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGata            <  2:2745601/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       >  2:2729745/1‑100 (MQ=255)
                                                                                            cccGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTc       >  2:1111020/1‑100 (MQ=255)
                                                                                            cccGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTc       >  1:2628064/1‑100 (MQ=255)
                                                                                            cccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTc       >  2:2619651/1‑100 (MQ=255)
                                                                                            cccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTc       >  1:1110025/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)
                                                                                                  |                                                                                                  
AGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGC  >  NC_000913/4478916‑4479110

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: