Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,247,837 Δ1 bp 15.5% coding (65/663 nt) yqjA → general envelope maintenance protein; DedA family inner membrane protein; putative multidrug efflux transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,247,8360C.15.5% 313.9 / 59.2 110coding (64/663 nt)yqjAgeneral envelope maintenance protein; DedA family inner membrane protein; putative multidrug efflux transporter
Reads supporting (aligned to +/- strand):  ref base C (51/40);  new base . (10/7);  total (62/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

Only 100 of 116 total aligned reads displayed.

CAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGG  >  NC_000913/3247738‑3247935
                                                                                                  |                                                                                                    
cAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATcc                                                                                                     <  2:1238653/100‑1 (MQ=255)
cAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATcc                                                                                                     >  1:353692/1‑100 (MQ=255)
cAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATcc                                                                                                     <  2:292644/100‑1 (MQ=255)
 aGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCa                                                                                                    >  2:2157245/1‑100 (MQ=255)
       gCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCgatg                                                                                             >  1:1249270/1‑100 (MQ=255)
        cGTTAACCGAGTCACCAGGAATAATGAATGGAACTTTTGACCCGATTGCGACAAGCCCTGTGGGCGGAGGATTCTGAAACCCTGGCCAAACCTTCgatg                                                                                             <  2:520889/99‑1 (MQ=255)
           tAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTg                                                                                         >  1:1417955/1‑100 (MQ=255)
             aCCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGc                                                                                       <  1:1414715/100‑1 (MQ=255)
              ccGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCa                                                                                      >  1:1975884/1‑100 (MQ=255)
              ccGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCATCGATGGATTGGcc                                                                                      <  2:141231/100‑2 (MQ=255)
              ccGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCATCGATGGATTGGcc                                                                                      <  1:889891/100‑2 (MQ=255)
              ccGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACACTGGCCAATCCATCGATGATTGGCa                                                                                      >  2:1747166/1‑100 (MQ=255)
                gATTCACCAGGAATAATTAATGGAACTTGTGTCCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTCGCCAATCATCGATGGATTGGccag                                                                                    <  2:553843/100‑4 (MQ=255)
                gATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCATCGATGGATTGGccag                                                                                    <  1:1418534/100‑4 (MQ=255)
                gATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCATCGATGGATTGGccag                                                                                    <  1:1574649/100‑4 (MQ=255)
                 aTTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCAtgt                                                                                   >  1:493742/1‑100 (MQ=255)
                 aTTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCATCGATGGATTGGccagg                                                                                   <  2:5522/100‑5 (MQ=255)
                  ttCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCa                                                                                      <  1:1031881/96‑1 (MQ=255)
                  ttCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCATCGATGGATTGGccaggg                                                                                  <  2:831962/100‑6 (MQ=38)
                  ctCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCAtgtt                                                                                  <  2:451627/99‑1 (MQ=255)
                    cACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCAtgttgt                                                                                >  1:752709/1‑100 (MQ=255)
                    cACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCAtgttgt                                                                                >  1:632308/1‑100 (MQ=255)
                     aCCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTa                                                                               >  1:907115/1‑100 (MQ=255)
                     aCCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTa                                                                               <  2:2115541/100‑1 (MQ=255)
                     aCCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCACTGTGGGCGCAGGATTTTGAAACCCTGGCCACGCCATCGATGATTGGCATGATGTa                                                                               >  2:105485/1‑100 (MQ=255)
                       cAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTAtt                                                                             >  1:1008683/1‑100 (MQ=255)
                        aGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTAttt                                                                            <  1:716456/100‑1 (MQ=255)
                        aGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTAttt                                                                            >  1:496779/1‑100 (MQ=255)
                         ggAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTAtttt                                                                           >  2:711305/1‑100 (MQ=255)
                         ggAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTAtttt                                                                           <  1:177596/100‑1 (MQ=255)
                         ggAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTAtttt                                                                           >  1:720818/1‑100 (MQ=255)
                           aataatGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGt                                                                         <  1:1770011/100‑1 (MQ=255)
                            ataatGAATGGAACTTGGAACCCAAGTGCTGCAAGCAGTGCGGGCGCAGGATTTGGAAACCCTGGCCAATACATCGATGATTGGCAGGTGGAATTTCGt                                                                         <  2:835316/99‑1 (MQ=255)
                             taatGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCt                                                                       <  2:2078985/100‑1 (MQ=255)
                             taatGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAAGCCATCGATGATTGGCATGTTGTATTTTGTCt                                                                       <  2:1350854/100‑1 (MQ=255)
                              aatgaatgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCtt                                                                      >  1:296323/1‑100 (MQ=255)
                               atgaatgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCACTCCATCGATGATTGGCATGTTGTATTTTGTCttg                                                                     >  2:1873258/1‑100 (MQ=255)
                               atgaatgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCttg                                                                     <  1:2274519/100‑1 (MQ=255)
                                 gaatgGAACGTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCttgtt                                                                   <  1:1387765/100‑1 (MQ=255)
                                   atgGAACTTTTTACCCAATTGCTGCAAGCTCTGGCGGCGCAGGATTTTGAAAGCTTGGCCAAACCATCGATAATTGGCATGTTGTATTTTGTCttgtttg                                                                 <  2:850540/100‑1 (MQ=255)
                                   atgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCttgtttg                                                                 <  1:581254/100‑1 (MQ=255)
                                   atgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCttgtttg                                                                 >  2:2008941/1‑100 (MQ=255)
                                    tgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCttgtttgt                                                                >  2:1221688/1‑100 (MQ=255)
                                    tgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCttgtttgt                                                                <  1:1688174/100‑1 (MQ=255)
                                      gAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCATGATTTTGAAACCATGGCCAATACGTCGATGATTGGCATGTTGTAGTTTGTCTAGTTTGTa                                                               >  2:1230624/1‑99 (MQ=255)
                                      gAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTaa                                                              >  1:1536468/1‑100 (MQ=255)
                                      gAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTaa                                                              >  1:705431/1‑100 (MQ=255)
                                         cTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAAttt                                                           >  1:418092/1‑100 (MQ=255)
                                           tttGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTg                                                         <  2:375993/100‑1 (MQ=255)
                                               aCCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTcc                                                     <  2:2388497/100‑1 (MQ=255)
                                               aCCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTcc                                                     >  1:904965/1‑100 (MQ=255)
                                                cccAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCt                                                    <  1:731788/100‑1 (MQ=255)
                                                  cAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTg                                                  >  2:274131/1‑100 (MQ=255)
                                                  cAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTg                                                  >  1:1040238/1‑100 (MQ=255)
                                                     ttGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGaa                                                <  2:841862/99‑1 (MQ=255)
                                                         tgcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACgg                                           >  1:1476729/1‑100 (MQ=255)
                                                         tgcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACgg                                           <  1:677545/100‑1 (MQ=255)
                                                         tgcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACgg                                           <  1:1794760/100‑1 (MQ=255)
                                                         tgcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACg                                            >  2:1246327/1‑99 (MQ=255)
                                                          gcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGgc                                          <  1:324280/100‑1 (MQ=255)
                                                          gcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGgc                                          >  1:451211/1‑100 (MQ=255)
                                                            aaGCCCTGTGGGGGCAGGATTTTGAAACCCTGGCCATTCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGgc                                          <  2:1049277/98‑1 (MQ=255)
                                                            aaGCCCTGTGGGCTCAGGATTTTGAAACCCTGTCCAATCCATCGATGATTGGCATGTTGTATTTTGACTTGTTTGTAATTGTGTTCCTTGAAAACGgctt                                        <  2:118203/100‑1 (MQ=255)
                                                            aaGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGgctt                                        <  1:1876595/100‑1 (MQ=255)
                                                            aaGCCCTGCGGGCGCAGGATATTGAAACCCTGGCCAACCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGgc                                          >  2:539534/1‑98 (MQ=255)
                                                             aGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGgcttg                                       <  2:1778642/100‑1 (MQ=255)
                                                                     gggCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCcggc                               <  1:627490/100‑1 (MQ=255)
                                                                       gcgcAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCcggcgg                             >  1:1807353/1‑100 (MQ=255)
                                                                       gcgcAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCcggcgg                             >  1:396046/1‑100 (MQ=255)
                                                                        ctcAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCggcggc                            <  2:1441946/98‑1 (MQ=255)
                                                                        cgcAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAACACGGCTTGCTTCCggcggc                            >  1:1390260/1‑100 (MQ=255)
                                                                        cgcAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCggcggc                            >  2:220446/1‑100 (MQ=255)
                                                                        cgcAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTGGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCcggcgg                             >  2:2024538/1‑99 (MQ=255)
                                                                         gcAGGATTTTGAAACCCTTGCCAATCAATCGATGATTGGCATGTTGTATTTTGCCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCcggcgg                             <  1:1957662/98‑1 (MQ=255)
                                                                         gcAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGcc                           <  2:728978/100‑1 (MQ=255)
                                                                            ggATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCttt                        <  2:2299776/100‑1 (MQ=255)
                                                                            ggATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCttt                        >  1:2047280/1‑100 (MQ=255)
                                                                            ggATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCttt                        <  1:1000536/100‑1 (MQ=255)
                                                                              aTTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCttttt                      >  2:1395347/1‑100 (MQ=255)
                                                                               gtttGAAACCCTGGCCAATCCATCGATGATTGGCATGTGGTATTTTGTCTTGTTTGTCATCTTGTTCGTTGATTACGGCTTGCTTCCGGCGGCCttttt                      >  2:1672895/2‑99 (MQ=255)
                                                                                tttGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTAc                    >  1:241194/1‑100 (MQ=255)
                                                                                tttGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTAc                    >  1:2264623/1‑100 (MQ=255)
                                                                                 ttGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTAcc                   >  1:699586/1‑100 (MQ=255)
                                                                                 ttGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTAcc                   >  2:2302968/1‑100 (MQ=255)
                                                                                     aaCCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGc               >  1:1634476/1‑100 (MQ=255)
                                                                                      naacatgCCAATCATCGATGGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCg              >  2:1725070/7‑100 (MQ=38)
                                                                                      caacatgCCAATCATCGATGGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCg              >  1:1582691/7‑100 (MQ=38)
                                                                                       cccTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCg              >  1:914402/1‑99 (MQ=255)
                                                                                        ccTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGAc            >  2:2322399/1‑100 (MQ=255)
                                                                                        acatgCCAATCATCGATGGATTGGCATGTTGTATTTTGTCTTGTTTGTTATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTTCCg                  >  2:1300183/5‑94 (MQ=37)
                                                                                        acatgCCAATCATCGATGGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGAc            >  1:1244569/5‑100 (MQ=255)
                                                                                         catgCCAATCATCGATGGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGCGCGACa           >  1:1633347/4‑100 (MQ=255)
                                                                                          atgCCAATCATCGGTGGATTGGCATGTTGTATTTTGTCTTGTTTGTAGTTTTGTTCCTTGGAAACGGTTTGATTCAGGAGGCCTTTTTACCGGGCGaa            >  2:1903500/3‑97 (MQ=255)
                                                                                          atgCCAATCATCGATGGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAg          >  1:619144/3‑100 (MQ=255)
                                                                                          atgCCAATCATCGATGGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCAGGCGGCCTTTTTACCGGGCGACAg          >  2:1480370/3‑100 (MQ=255)
                                                                                           tgCCAATCATCGATGGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGt         >  1:1795443/2‑100 (MQ=255)
                                                                                            gCCAATCCATCGATGATTGGCGTGTTGTATTGTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGtt        <  2:1061255/100‑1 (MQ=255)
                                                                                                aTCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACt    <  2:2047262/100‑1 (MQ=255)
                                                                                                  ccATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTgg  >  1:1241349/1‑100 (MQ=255)
                                                                                                  ccATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTgg  <  1:1495382/100‑1 (MQ=255)
                                                                                                  |                                                                                                    
CAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGG  >  NC_000913/3247738‑3247935

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: