Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,663,553 A→C 10.1% intergenic (+26/+337) mdtF → / ← gadW anaerobic multidrug efflux transporter, ArcA‑regulated/transcriptional activator of gadA and gadBC; repressor of gadX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,663,5530AC10.1% 311.0 / 22.0 117intergenic (+26/+337)mdtF/gadWanaerobic multidrug efflux transporter, ArcA‑regulated/transcriptional activator of gadA and gadBC; repressor of gadX
Reads supporting (aligned to +/- strand):  ref base A (56/48);  new base C (3/9);  total (60/57)
Fisher's exact test for biased strand distribution p-value = 7.20e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

Only 100 of 118 total aligned reads displayed.

CAAAACGCGGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGACCCACA  >  NC_000913/3663405‑3663702
                                                                                                                                                    |                                                                                                                                                     
cAAAACGCGGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACAt                                                                                                                                                      >  2:638500/1‑150 (MQ=255)
   aaCGCGGTAGGTACCGGCGTGATGGGCTGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGGGGAACAGCTCTTTGCCCGCTTTAAAAATGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGaa                                                                                                                                                   <  1:1134179/150‑1 (MQ=255)
    aCGCGGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAg                                                                                                                                                  <  2:1610822/150‑1 (MQ=255)
     cgcgGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAgg                                                                                                                                                 <  1:266254/150‑1 (MQ=255)
         gtaggtaCCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCa                                                                                                                                             >  1:775695/1‑150 (MQ=255)
          taggtaCCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAg                                                                                                                                            <  2:1316485/150‑1 (MQ=255)
            ggtaCCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCTGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGcc                                                                                                                                          <  1:1458373/150‑1 (MQ=255)
              taCCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGccct                                                                                                                                        <  2:829334/150‑3 (MQ=255)
               aCCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTa                                                                                                                                       >  2:118266/1‑150 (MQ=255)
                ccGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTAc                                                                                                                                      >  1:933786/1‑150 (MQ=255)
                ccGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTAc                                                                                                                                      <  2:835611/150‑1 (MQ=255)
                ccGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTAc                                                                                                                                      <  1:1613909/150‑1 (MQ=255)
                ccGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGCTGAACCTTCATGTTCAACCTTAc                                                                                                                                      <  2:148098/150‑1 (MQ=255)
                ccGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGCTGAACCTTCATGTTCAACCTTAc                                                                                                                                      <  1:1663048/150‑1 (MQ=255)
                  ggCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTActc                                                                                                                                    <  1:892449/150‑1 (MQ=255)
                  ggCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGCTGAACCTTCATGTTCAACCTTActc                                                                                                                                    <  1:1568082/150‑1 (MQ=255)
                  ggCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGCTGAACCTTCATGTTCAACCTTActc                                                                                                                                    <  1:322503/150‑1 (MQ=255)
                   gCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGCTGAACCTTCATGTTCAACCTTActct                                                                                                                                   <  2:1593471/150‑1 (MQ=39)
                   gCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGCTGAACCTTCATGTTCAACCTTActct                                                                                                                                   <  1:860814/150‑1 (MQ=39)
                   gCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGCTGAACCTTCATGTTCAACCTTActct                                                                                                                                   <  1:1302946/150‑1 (MQ=39)
                    ggggctGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCtt                                                                                                                                  <  1:350500/145‑1 (MQ=255)
                     gTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCttt                                                                                                                                 >  2:1139838/1‑150 (MQ=255)
                      tGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTc                                                                                                                                <  1:570824/150‑1 (MQ=255)
                         tGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTg                                                                                                                             >  2:1364488/1‑150 (MQ=255)
                           ggCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCt                                                                                                                           >  1:1493685/1‑150 (MQ=255)
                                   gTTTGCCGCAACAGTGCTGGCAATTTACTTCGATCCTGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTc                                                                                                                   <  1:1590920/150‑1 (MQ=255)
                                      tGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTcacc                                                                                                                >  2:50463/1‑150 (MQ=255)
                                      tGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGGCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTcacc                                                                                                                <  2:350608/150‑1 (MQ=255)
                                       gCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTcacca                                                                                                               <  1:700398/150‑1 (MQ=255)
                                          gCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCCCCACTc                                                                                                            >  1:1487911/1‑150 (MQ=255)
                                          gCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTc                                                                                                            >  1:1494714/1‑150 (MQ=255)
                                              cAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTcaaca                                                                                                        <  2:1518880/150‑1 (MQ=255)
                                                 tGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAAcc                                                                                                     <  2:700289/150‑1 (MQ=255)
                                                   cTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAg                                                                                                   >  2:215341/1‑150 (MQ=255)
                                                    tGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGc                                                                                                  <  2:1644411/150‑1 (MQ=255)
                                                     ggCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCt                                                                                                 <  2:19639/150‑1 (MQ=255)
                                                     ggCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCt                                                                                                 >  2:1605138/1‑150 (MQ=255)
                                                       cAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTaa                                                                                               >  1:871051/1‑150 (MQ=255)
                                                            tACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTc                                                                                          >  2:1373285/1‑150 (MQ=255)
                                                            tACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACAATTCAACACTCAACAACCAGCTAACCggc                                                                                          >  1:300390/1‑148 (MQ=255)
                                                               ttCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAgg                                                                                       <  1:638500/150‑1 (MQ=255)
                                                                     ccGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATAcc                                                                                 <  1:1139838/150‑1 (MQ=255)
                                                                      cGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGTTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCa                                                                                >  2:1239501/1‑150 (MQ=255)
                                                                       ggTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAg                                                                               >  2:583307/1‑150 (MQ=255)
                                                                       ggTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAg                                                                               >  1:185157/1‑150 (MQ=255)
                                                                         tCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGc                                                                             <  1:494732/150‑1 (MQ=255)
                                                                            ttttCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTa                                                                          <  2:691757/150‑1 (MQ=255)
                                                                             tttCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAg                                                                         <  1:155518/150‑1 (MQ=255)
                                                                                 gTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTaca                                                                     >  1:1484229/1‑150 (MQ=255)
                                                                                 gTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTATAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTaca                                                                     >  1:833481/1‑150 (MQ=255)
                                                                                  ttGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTacac                                                                    <  2:417450/150‑1 (MQ=255)
                                                                                    gTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAg                                                                  <  2:1493685/150‑1 (MQ=255)
                                                                                     tAGTGGTACATCTCTTTTCCCTCTTTAAATAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCTCCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGc                                                                 <  1:1201688/150‑1 (MQ=255)
                                                                                     tAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGc                                                                 <  1:453734/150‑1 (MQ=255)
                                                                                       gTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCAc                                                               >  2:437784/1‑150 (MQ=255)
                                                                                         ggAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCTCTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTGCACAGCACTg                                                             >  2:210623/1‑150 (MQ=255)
                                                                                         ggAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTg                                                             >  1:509792/1‑150 (MQ=255)
                                                                                         ggAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGAACtt                                                             >  2:1556709/1‑149 (MQ=255)
                                                                                          taaCATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGa                                                            <  1:707022/149‑1 (MQ=255)
                                                                                          gAACATCTCTTTGCACGCTTTAAAAAAGCGTAAAGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTCCTAACCATTCACCACTCAACAAACAGCTAACCGTCAGGAATACCAGGCATAGTACACAGCACtta                                                            >  1:542492/1‑148 (MQ=255)
                                                                                           aaCATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGaa                                                           <  1:1667992/150‑1 (MQ=255)
                                                                                            aCATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCACCAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGaaa                                                          >  2:900727/1‑150 (MQ=255)
                                                                                             cATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAg                                                         >  2:1660753/1‑150 (MQ=255)
                                                                                             cATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAg                                                         <  2:909512/150‑1 (MQ=255)
                                                                                              aTCTCTTTGCACGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGt                                                        <  1:1387220/150‑1 (MQ=255)
                                                                                                ctctTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAg                                                      <  2:691127/150‑1 (MQ=255)
                                                                                                 tctTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGTGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACGACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGa                                                     <  2:397556/150‑1 (MQ=255)
                                                                                                    ttGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAAc                                                  <  1:249492/150‑1 (MQ=255)
                                                                                                     tGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAAcc                                                 <  2:1514815/150‑1 (MQ=255)
                                                                                                            tttAAATAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGGTCAGCCTTACTCTTTCCTGCTAACCATTCACGACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAg                                          <  1:645726/150‑1 (MQ=255)
                                                                                                             ttAAAAAATCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGc                                         <  2:1247661/150‑1 (MQ=255)
                                                                                                             ttAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCCGCTAACCGTCAGGAATACCAGGCTTAGTACCCAGCACTGAAAGTAGAAACCACTCCAGc                                         >  2:849351/1‑150 (MQ=255)
                                                                                                                 aaaaGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGCAACCACTCCAGCCAtt                                     >  2:883611/1‑150 (MQ=255)
                                                                                                                      cGTATCGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGtca                                <  1:692558/150‑1 (MQ=255)
                                                                                                                          aCGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCaa                            >  1:708331/1‑150 (MQ=255)
                                                                                                                            gtgtAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACACCCAGCTAACCGTCAGGAATACCAGGCTTAGTACAAAGCACTGAAAGTAGAAACAACTCCAGCCATTCGTCATCAAtt                          >  2:1337670/1‑150 (MQ=255)
                                                                                                                            gtgtAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAAtt                          >  2:590820/1‑150 (MQ=255)
                                                                                                                             tgtAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAAttt                         >  1:756736/1‑150 (MQ=255)
                                                                                                                              gtAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTg                        >  1:1469285/1‑150 (MQ=255)
                                                                                                                                aaaTGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTATCCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGaa                      >  2:400805/1‑150 (MQ=255)
                                                                                                                                aaaTGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGaa                      >  1:728368/1‑150 (MQ=255)
                                                                                                                                 ggaaagagTAAGGCTGAACCTTCATGTTCAACCTTACTCTTTCCTTCTAACCATTCACCACTCAACCACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCCCTCCAGCCATTCGTCATCAATTTGaac                     >  2:1602447/5‑150 (MQ=37)
                                                                                                                                 aaTGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCATGCTAACCAATCACCACTAAACAAACAGATTACCGTCAGTAATAACAGTCTTAGTACACATCACTGTAAGTATAAACCACTCGAGCCAGTCGTCATCACATTGaac                     >  2:291129/1‑150 (MQ=255)
                                                                                                                                  aTGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGaaca                    <  1:50463/150‑1 (MQ=255)
                                                                                                                                   tGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCCACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAACCCACTCCAGCCATTCGTCATCAATTTGaacac                   >  2:935977/1‑149 (MQ=255)
                                                                                                                                    gagaGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTTGTCATCAATTTGaacaac                  >  2:986866/1‑150 (MQ=255)
                                                                                                                                     agagTAAGGTTGAACTTGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGTATACCAGGCTTAGTACACAGCACTGAAAGTAGAAGCCACTACAGCCATTCGTCATCAATTTGaacaacc                 >  1:1076694/1‑149 (MQ=255)
                                                                                                                                     agagTAAGGCTGAACCTTAATCTGCAACCTTACTCTTTCCAGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGATTACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCCCTCAAGCCCTTCGTCATCAATTTGAaaaaca                 >  2:1175727/1‑150 (MQ=255)
                                                                                                                                      gagTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAacac                >  1:614701/1‑150 (MQ=255)
                                                                                                                                      gagTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAacac                <  1:437784/150‑1 (MQ=255)
                                                                                                                                       agTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGTAATACCAGGCTTAGTACACAGCACTGAAATTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAAcaca               <  1:1660753/150‑1 (MQ=255)
                                                                                                                                        gTAAGGCTGAACCTTCATGTTCAACCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACaa              >  1:426742/1‑150 (MQ=255)
                                                                                                                                          aaGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTACCCATTCACCACTCACCAACCAGCTACCCGTCAGGAATACCAGGCTTAGTACACAGCCCTGAAAGTAGACACCCCTCCAGCCATTCGTCATCAATTTGcacaacacaaga            >  1:1717682/1‑148 (MQ=255)
                                                                                                                                            ggTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTTGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACt          >  1:1510972/1‑150 (MQ=255)
                                                                                                                                             gTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCCGCTAACCGTCAGGAATACCAGGCTAAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACtc         >  2:1348457/1‑149 (MQ=255)
                                                                                                                                             gTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTg         >  2:1433159/1‑150 (MQ=255)
                                                                                                                                                gAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAAttt                         >  1:1699312/1‑131 (MQ=255)
                                                                                                                                                gAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAAttt                         <  2:1699312/131‑1 (MQ=255)
                                                                                                                                                gAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGAcc      >  2:543610/1‑150 (MQ=255)
                                                                                                                                                gAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGAcc      <  2:532774/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAAAACGCGGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGACCCACA  >  NC_000913/3663405‑3663702

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: