Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,247,844:1 +G 15.2% coding (72/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,8431.G15.2% 378.8 / 70.5 112coding (71/663 nt)yqjAgeneral envelope maintenance protein; DedA family inner membrane protein; putative multidrug efflux transporter
Reads supporting (aligned to +/- strand):  ref base . (53/42);  new base G (10/7);  total (63/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

Only 100 of 117 total aligned reads displayed.

GCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTG  >  NC_000913/3247745‑3247940
                                                                                                   |                                                                                                 
gCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCgatg                                                                                                  >  1:1249270/1‑100 (MQ=255)
 cGTTAACCGAGTCACCAGGAATAATGAATGGAACTTTTGACCCGATTGCGACAAGCCCTGTGGGCGGAGGATTCTGAAACCCTGGCCAAACCTTCgatg                                                                                                  <  2:520889/99‑1 (MQ=255)
      aCCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGc                                                                                            <  1:1414715/100‑1 (MQ=255)
       ccGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCa                                                                                           >  1:1975884/1‑100 (MQ=255)
       ccGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCATCGATGGATTGGcc                                                                                           <  2:141231/100‑2 (MQ=255)
       ccGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACACTGGCCAATCCATCGATGATTGGCa                                                                                           >  2:1747166/1‑100 (MQ=255)
         gATTCACCAGGAATAATTAATGGAACTTGTGTCCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTCGCCAATCATCGATGGATTGGccag                                                                                         <  2:553843/100‑4 (MQ=255)
         gATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCAtg                                                                                         >  2:1954353/1‑100 (MQ=255)
         gATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCATCGATGGATTGGccag                                                                                         <  1:1418534/100‑4 (MQ=255)
         gATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCATCGATGGATTGGccag                                                                                         <  1:1574649/100‑4 (MQ=255)
          aTTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCAtgt                                                                                        <  1:1820955/100‑1 (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:632308/1‑100 (MQ=255)
             cACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCAtgttgt                                                                                     >  1:752709/1‑100 (MQ=255)
             cACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCAtgttag                                                                                     >  1:451629/1‑98 (MQ=255)
              agcAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGAGGGCGCAGGCTTTTGAAACCCTGAACAGTCCATCGACGATAGGCATGTTGTa                                                                                    >  2:105527/3‑100 (MQ=255)
              aCCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTa                                                                                    >  1:907115/1‑100 (MQ=255)
              aCCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTa                                                                                    <  2:2115541/100‑1 (MQ=255)
              aCCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTa                                                                                    >  1:247321/1‑100 (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                                                                                <  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                                                                      >  2:2008941/1‑100 (MQ=255)
                            atgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCttgtttg                                                                      <  2:2229827/100‑1 (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)
                                  cTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAAttt                                                                >  1:418092/1‑100 (MQ=255)
                                  cTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAAttt                                                                >  1:744719/1‑100 (MQ=255)
                                    tttGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTg                                                              <  2:375993/100‑1 (MQ=255)
                                        aCCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTcc                                                          >  1:904965/1‑100 (MQ=255)
                                        aCCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTcc                                                          <  2:2388497/100‑1 (MQ=255)
                                         cccAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCt                                                         <  1:731788/100‑1 (MQ=255)
                                           cAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTg                                                       >  1:1040238/1‑100 (MQ=255)
                                           cAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTg                                                       >  2:274131/1‑100 (MQ=255)
                                              ttGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGaa                                                     <  2:841862/99‑1 (MQ=255)
                                                  tgcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACgg                                                >  1:1476729/1‑100 (MQ=255)
                                                  tgcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACgg                                                >  1:6517/1‑100 (MQ=255)
                                                  tgcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACgg                                                <  1:677545/100‑1 (MQ=255)
                                                  tgcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACg                                                 >  2:1246327/1‑99 (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)
                                                           tgtgGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTcc                                       >  1:1278520/1‑100 (MQ=255)
                                                           tgtgGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTcc                                       <  1:1665072/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)
                                                                 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: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)
                                                                                cccTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCg                   >  1:914402/1‑99 (MQ=255)
                                                                                 ccTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGAc                 >  2:2322399/1‑100 (MQ=255)
                                                                                 acatgCCAATCATCGATGGATTGGCATGTTGTATTTTGTCTTGTTTGTTATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTTCCg                       >  2:1300183/5‑94 (MQ=37)
                                                                                  catgCCAATCATCGATGGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACa                >  1:1493250/4‑100 (MQ=255)
                                                                                  catgCCAATCATCGATGGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGCGCGACa                >  1:1633347/4‑100 (MQ=255)
                                                                                   atgCCAATCATCGATGGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAg               >  1:619144/3‑100 (MQ=255)
                                                                                   atgCCAATCATCGATGGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCAGGCGGCCTTTTTACCGGGCGACAg               >  2:1480370/3‑100 (MQ=255)
                                                                                    tgCCAATCATCGATGGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGt              >  1:1795443/2‑100 (MQ=255)
                                                                                         aTCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACt         <  2:2047262/100‑1 (MQ=255)
                                                                                           ccATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTgg       >  1:1241349/1‑100 (MQ=255)
                                                                                           ccATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTgg       <  1:1495382/100‑1 (MQ=255)
                                                                                            cATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGt      >  2:1712849/1‑100 (MQ=255)
                                                                                            cATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGt      >  1:2309852/1‑100 (MQ=255)
                                                                                             aTCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTgg       <  2:1138719/98‑1 (MQ=255)
                                                                                                gatgatTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTg  <  2:1098443/100‑1 (MQ=255)
                                                                                                gatgatTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTg  <  1:1346158/100‑1 (MQ=255)
                                                                                                gatgatTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTg  >  1:1098448/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
GCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTG  >  NC_000913/3247745‑3247940

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: