Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,663,547 T→C 10.6% intergenic (+20/+343) 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,5470TC10.6% 296.4 / 22.7 115intergenic (+20/+343)mdtF/gadWanaerobic multidrug efflux transporter, ArcA‑regulated/transcriptional activator of gadA and gadBC; repressor of gadX
Reads supporting (aligned to +/- strand):  ref base T (56/46);  new base C (3/9);  total (59/56)
Fisher's exact test for biased strand distribution p-value = 6.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-01

Only 100 of 116 total aligned reads displayed.

GGCGCGCAAAACGCGGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTG  >  NC_000913/3663399‑3663695
                                                                                                                                                    |                                                                                                                                                    
ggCGCGCAAAACGCGGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGtt                                                                                                                                                     <  1:1579525/150‑1 (MQ=255)
         aaCGCGGTAGGTACCGGCGTGATGGGCTGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGGGGAACAGCTCTTTGCCCGCTTTAAAAATGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGaa                                                                                                                                            <  1:1134179/150‑1 (MQ=255)
          aCGCGGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAg                                                                                                                                           <  2:1610822/150‑1 (MQ=255)
          aCGCGGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAg                                                                                                                                           >  1:909512/1‑150 (MQ=255)
           cgcgGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAgg                                                                                                                                          <  1:266254/150‑1 (MQ=255)
               gtaggtaCCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCa                                                                                                                                      >  1:775695/1‑150 (MQ=255)
                taggtaCCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAg                                                                                                                                     <  2:1316485/150‑1 (MQ=255)
                    taCCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGccct                                                                                                                                 <  2:829334/150‑3 (MQ=255)
                     aCCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTa                                                                                                                                >  2:118266/1‑150 (MQ=255)
                      ccGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTAc                                                                                                                               <  2:835611/150‑1 (MQ=255)
                      ccGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTAc                                                                                                                               <  2:674534/150‑1 (MQ=255)
                      ccGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTAc                                                                                                                               >  1:933786/1‑150 (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                                                                                                                            <  1:860814/150‑1 (MQ=39)
                         gCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGCTGAACCTTCATGTTCAACCTTActct                                                                                                                            <  1:1302946/150‑1 (MQ=39)
                         gCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGCTGAACCTTCATGTTCAACCTTActct                                                                                                                            <  1:764627/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)
                                 ggCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCt                                                                                                                    <  1:24949/150‑1 (MQ=255)
                                    gggATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAAc                                                                                                                 <  1:374636/150‑1 (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)
                                                          tGGCAATATACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGc                                                                                           >  2:509529/1‑150 (MQ=255)
                                                           ggCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCt                                                                                          >  2:1605138/1‑150 (MQ=255)
                                                                 ttACTTCGTTCCGTTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACGATTCACCACTCAACAACCAGCTAACCGt                                                                                    <  1:540215/150‑1 (MQ=255)
                                                                  tACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTc                                                                                   >  2:1373285/1‑150 (MQ=255)
                                                                  tACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACAATTCAACACTCAACAACCAGCTAACCggc                                                                                   >  1:300390/1‑148 (MQ=255)
                                                                     ttCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAgg                                                                                <  1:638500/150‑1 (MQ=255)
                                                                        gTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAAc                                                                                                                                                 >  2:213627/1‑82 (MQ=255)
                                                                        gTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAAc                                                                                                                                                 <  1:213627/82‑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)
                                                                                             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:909512/150‑1 (MQ=255)
                                                                                                   cATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAg                                                  >  2:1660753/1‑150 (MQ=255)
                                                                                                    aTCTCTTTGCACGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGt                                                 <  1:1387220/150‑1 (MQ=255)
                                                                                                      ctctTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAg                                               >  1:1499104/1‑150 (MQ=255)
                                                                                                      ctctTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAg                                               <  2:691127/150‑1 (MQ=255)
                                                                                                       tctTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGTGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACGACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGa                                              <  2:397556/150‑1 (MQ=255)
                                                                                                           tGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAAcc                                          <  2:1514815/150‑1 (MQ=255)
                                                                                                                  tttAAATAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGGTCAGCCTTACTCTTTCCTGCTAACCATTCACGACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAg                                   <  1:645726/150‑1 (MQ=255)
                                                                                                                     aaaaaaGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCa                                >  2:1676569/1‑150 (MQ=255)
                                                                                                                       aaaaGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGCAACCACTCCAGCCAtt                              >  2:883611/1‑150 (MQ=255)
                                                                                                                       aaaaGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAACCCACTCCAGCCAtt                              >  2:872103/1‑150 (MQ=255)
                                                                                                                            cGTATCGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGtca                         <  1:692558/150‑1 (MQ=255)
                                                                                                                             gTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGtcat                        >  2:338790/1‑150 (MQ=255)
                                                                                                                              tAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGtcatc                       >  2:1585500/1‑150 (MQ=255)
                                                                                                                               aaCGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGtcatca                      >  1:1615899/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)
                                                                                                                                      aaaTGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGaa               >  1:659661/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)
                                                                                                                                       aaTGAGAGTAAGGATGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGaac              <  1:1436778/150‑1 (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:437784/150‑1 (MQ=255)
                                                                                                                                              gTAAGGCTGAACCTTCATGTTCAACCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACaa       >  1:426742/1‑150 (MQ=255)
                                                                                                                                                aaGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTACCCATTCACCACTCACCAACCAGCTACCCGTCAGGAATACCAGGCTTAGTACACAGCCCTGAAAGTAGACACCCCTCCAGCCATTCGTCATCAATTTGcacaacacaaga     >  1:1717682/1‑148 (MQ=255)
                                                                                                                                                   gTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCCGCTAACCGTCAGGAATACCAGGCTAAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACtc  >  2:1348457/1‑149 (MQ=255)
                                                                                                                                                   gTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTg  >  2:1433159/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGCGCGCAAAACGCGGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTG  >  NC_000913/3663399‑3663695

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: