Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,479,016 G→C 17.8% intergenic (+127/+18) 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,0160GC17.8% 275.1 / 49.1 107intergenic (+127/+18)rraB/yjgMprotein inhibitor of RNase E/GNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (48/39);  new base C (8/11);  total (56/51)
Fisher's exact test for biased strand distribution p-value = 3.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 110 total aligned reads displayed.

GTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAA  >  NC_000913/4478917‑4479114
                                                                                                     |                                                                                                  
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:229094/1‑100 (MQ=255)
   ccTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAac                                                                                                 <  1:1729999/100‑1 (MQ=255)
    cTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAaca                                                                                                <  1:718902/100‑1 (MQ=255)
        aaCCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACAcgcg                                                                                            <  2:2645795/100‑1 (MQ=255)
        aaCCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACAcgcg                                                                                            <  2:1968143/100‑1 (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:992098/100‑1 (MQ=39)
           cTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGc                                                                                         <  1:261301/100‑1 (MQ=39)
           cTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGc                                                                                         >  1:715053/1‑100 (MQ=255)
             ggTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCAt                                                                                       <  2:299367/100‑1 (MQ=38)
             ggTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAt                                                                                       >  1:1897854/1‑100 (MQ=255)
             ggTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAt                                                                                       >  2:513350/1‑100 (MQ=255)
             ggTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAt                                                                                       >  1:780655/1‑100 (MQ=255)
               tGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttt                                                                                     >  1:2745619/1‑100 (MQ=255)
               tGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttt                                                                                     <  2:137350/100‑1 (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)
                   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: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)
                                                                  cgtacg‑‑tTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAg                                  <  2:109810/100‑1 (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: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)
                                                                                   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                <  2:2745601/100‑1 (MQ=255)
                                                                                      gAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGata                <  2:238271/100‑1 (MQ=255)
                                                                                      gAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGata                <  1:1358872/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)
                                                                                                 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)
                                                                                                     |                                                                                                  
GTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACG‑‑TTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAA  >  NC_000913/4478917‑4479114

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: