Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,663,564 G→A 6.3% intergenic (+37/+326) 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,5640GA 6.3% 506.1 / 10.3 157intergenic (+37/+326)mdtF/gadWanaerobic multidrug efflux transporter, ArcA‑regulated/transcriptional activator of gadA and gadBC; repressor of gadX
Reads supporting (aligned to +/- strand):  ref base G (76/71);  new base A (7/3);  total (83/74)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

AGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGACCCACATTCCCGTAATC  >  NC_000913/3663416‑3663713
                                                                                                                                                    |                                                                                                                                                     
 ggtaCCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGcc                                                                                                                                                     >  2:1436008/1‑150 (MQ=255)
   taCCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCtt                                                                                                                                                   >  1:2099396/1‑150 (MQ=255)
   taCCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCtt                                                                                                                                                   >  2:183707/1‑150 (MQ=255)
    aCCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATAAGAGTAAGGCTGAACCTTCATGTTCAACCTTa                                                                                                                                                  <  1:658821/150‑1 (MQ=255)
     ccGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTAc                                                                                                                                                 >  2:391610/1‑150 (MQ=255)
     ccGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTAc                                                                                                                                                 <  2:2027833/150‑1 (MQ=255)
      cGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCTGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGCTGAACCTTCATGTTCAACCTTAct                                                                                                                                                <  2:2057470/150‑1 (MQ=255)
         cGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCtt                                                                                                                                             >  2:1766768/1‑150 (MQ=255)
          gTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCttt                                                                                                                                            >  1:2177904/1‑150 (MQ=255)
            gATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTcc                                                                                                                                          <  1:867269/150‑1 (MQ=255)
             cggggcgggcGGCTTGCCGCAACAGTGCTGGCAATTTACTTCTGTCCGGTCATTTTCGTTGTAGTGGACCATCTCTTTGCCCGCTTTAAAAAACCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCt                                                                                                                                         <  2:2081132/148‑1 (MQ=255)
              tGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACTTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTg                                                                                                                                        >  2:332974/1‑150 (MQ=255)
                ggCGGGATGTTTGCCGCAACAGTGCTGGCTATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGACCATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCTGCCTTACTCTTTCCTGCt                                                                                                                                      >  2:140494/1‑150 (MQ=255)
                   gggATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAAc                                                                                                                                   <  1:498871/150‑1 (MQ=255)
                     gATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTAGTAGTGGAACATCTCTTTGCCCGCATTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCa                                                                                                                                 >  1:373272/1‑150 (MQ=255)
                       tGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCAtt                                                                                                                               <  1:1983405/150‑1 (MQ=255)
                        gTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTc                                                                                                                              <  2:1975750/150‑1 (MQ=255)
                           tGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTcacc                                                                                                                           <  2:1163979/150‑1 (MQ=255)
                            gCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTcacca                                                                                                                          >  1:1028084/1‑150 (MQ=255)
                                 aaCAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTcaa                                                                                                                     >  1:970670/1‑150 (MQ=255)
                                 aaCAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTcaa                                                                                                                     <  2:167949/150‑1 (MQ=255)
                                  aCAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTcaac                                                                                                                    <  2:1902931/150‑1 (MQ=255)
                                     gTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCaacaac                                                                                                                 >  2:2032932/1‑150 (MQ=255)
                                     gTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCaacaac                                                                                                                 <  2:886660/150‑1 (MQ=255)
                                      tGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAAcc                                                                                                                <  2:442234/150‑1 (MQ=255)
                                        cTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAg                                                                                                              <  1:1002793/150‑1 (MQ=255)
                                         tGGCAATTTAATTCGTTCAGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCATGCTAACCATTCACCACTCAACAACCAGc                                                                                                             >  2:1008444/1‑150 (MQ=255)
                                           gCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAGATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTa                                                                                                           <  2:1078724/150‑1 (MQ=255)
                                           gCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTa                                                                                                           >  2:1904320/1‑150 (MQ=255)
                                            cAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCCCCACTCAACAACCAGCTaa                                                                                                          >  1:1861340/1‑150 (MQ=255)
                                               tttACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCg                                                                                                       >  2:123352/1‑150 (MQ=255)
                                                   cTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAg                                                                                                   >  1:765612/1‑150 (MQ=255)
                                                    ttCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCTTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACCACCAGCTCACCGTCAgg                                                                                                  >  1:364400/1‑150 (MQ=255)
                                                    ttCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAgg                                                                                                  <  1:410350/150‑1 (MQ=255)
                                                     tCGTTCCTGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGATGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGa                                                                                                 <  1:2123520/150‑1 (MQ=255)
                                                     tCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTATGCCCGCTTTAAAAAAGCGTAACCTGTAAATGAGAGTAAGCTTGAACATGAAGCTTCAGCCTTACTCTTTCCTGCTAACCATTAACCACTCAACAAACAGCTAACCTTCAgca                                                                                                 >  1:1974680/1‑148 (MQ=255)
                                                        ttCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATa                                                                                              >  2:904902/1‑150 (MQ=255)
                                                           cGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACAATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCa                                                                                           >  1:308158/1‑150 (MQ=255)
                                                               cTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCATCTAACCGTCAGGAATACCAGGCt                                                                                       >  2:595887/1‑150 (MQ=255)
                                                               cTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCt                                                                                       >  1:1942354/1‑150 (MQ=255)
                                                                    tCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTa                                                                                  <  2:96210/150‑1 (MQ=255)
                                                                    tCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTa                                                                                  <  1:1360614/150‑1 (MQ=255)
                                                                    tCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTa                                                                                  <  1:646899/150‑1 (MQ=255)
                                                                     cGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTac                                                                                 >  2:2033897/1‑150 (MQ=255)
                                                                     cGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTac                                                                                 >  1:2193982/1‑150 (MQ=255)
                                                                     cGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTac                                                                                 <  1:353919/150‑1 (MQ=255)
                                                                              ggAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTg                                                                        <  2:2199419/150‑1 (MQ=255)
                                                                              ggAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTg                                                                        >  1:2189063/1‑150 (MQ=255)
                                                                              ggAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTg                                                                        >  2:448579/1‑150 (MQ=255)
                                                                                aaCATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGGTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGaa                                                                      <  1:1779934/150‑1 (MQ=255)
                                                                                  cATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACCCAGCACTGAAAg                                                                    >  1:1127413/1‑150 (MQ=255)
                                                                                  cATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAg                                                                    <  2:1449980/150‑1 (MQ=255)
                                                                                   aTCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGt                                                                   <  2:767902/150‑1 (MQ=255)
                                                                                   aTCTCTTTGCCCGCTTTAAAAAAGCGTAACGCGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGACTTAGTACACAGCACTGAAAGt                                                                   <  2:1095039/150‑1 (MQ=255)
                                                                                    tctcTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTa                                                                  <  2:26862/150‑1 (MQ=255)
                                                                                    tctcTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTa                                                                  <  2:1044154/150‑1 (MQ=255)
                                                                                      tctTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGa                                                                >  1:1098322/1‑150 (MQ=255)
                                                                                          tGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAAcc                                                            >  1:2085721/1‑150 (MQ=255)
                                                                                          tGCACGCTTTAAAAAATCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGGTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCATGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAAcc                                                            <  1:1567749/150‑1 (MQ=255)
                                                                                           gCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCa                                                           >  2:1593502/1‑150 (MQ=255)
                                                                                                  ttAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCAATCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGc                                                    <  1:1164697/150‑1 (MQ=255)
                                                                                                    aaaaaaGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCa                                                  >  2:679457/1‑150 (MQ=255)
                                                                                                    aaaaaaGCGTAACGTGTAAATGAGAGTAAGGTTGAACAAGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCa                                                  >  1:32841/1‑150 (MQ=255)
                                                                                                     aaaaaGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCAt                                                 <  2:238321/150‑1 (MQ=255)
                                                                                                      aaaaGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACACCACTGAAAGTAGAAACCACTCCAGCCAtt                                                >  1:701682/1‑150 (MQ=255)
                                                                                                           cGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGtca                                           <  1:1658411/150‑1 (MQ=255)
                                                                                                             tAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGtcatc                                         >  1:689035/1‑150 (MQ=255)
                                                                                                               aCGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCaa                                       >  2:1628039/1‑150 (MQ=255)
                                                                                                               aCGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCaa                                       >  2:1486197/1‑150 (MQ=255)
                                                                                                                 gtgtAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAAtt                                     <  1:391610/150‑1 (MQ=255)
                                                                                                                 gtgtAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAAtt                                     <  2:2087513/150‑1 (MQ=255)
                                                                                                                  tgtAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAAttt                                    >  1:1502720/1‑150 (MQ=255)
                                                                                                                   gtAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTg                                   <  2:427305/150‑1 (MQ=255)
                                                                                                                    tAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGa                                  <  2:1705294/150‑1 (MQ=255)
                                                                                                                    tAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGa                                  <  2:920844/150‑1 (MQ=255)
                                                                                                                    caggaaagagTAAGGCTGAACCTTCATGTTCAACCTTACTCTTTCCTGCTAACCATTCACCCCGCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACCCCGCACTGAAGGTAGCAACCACTCCAGCCATACGTCATCAATTTGa                                  >  2:1583494/7‑150 (MQ=25)
                                                                                                                     aaaTGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACCGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGaa                                 >  2:368923/1‑150 (MQ=255)
                                                                                                                      ggaaagagTAAGGCTGATCCTTCATGTTCAACCTTACTCTTTCCTGCTAAACATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGcac                                >  1:1112102/5‑150 (MQ=38)
                                                                                                                      ggaaagagTAAGGCTGAACCTTCATGTTCAACCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATACGTCATCAATTTGaac                                >  1:1844862/5‑150 (MQ=38)
                                                                                                                      aaTGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGaac                                <  2:1776077/150‑1 (MQ=255)
                                                                                                                      aaTGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGaac                                <  2:1171658/150‑1 (MQ=255)
                                                                                                                      aaTGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGaac                                <  1:247802/150‑1 (MQ=255)
                                                                                                                      aaTGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGGCATCAATTTGaac                                >  1:1373710/1‑150 (MQ=255)
                                                                                                                        tGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCCTTCACCACTCAACACCCAGCTAACCGTCAGGAATACCAGGCTTAGTAAACAGCACGGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGaacaa                              >  2:1087724/1‑150 (MQ=255)
                                                                                                                         aagagTAAGGCTGATCCTTCATGTTCAACCTTACTCTTTCCTGCTAACAATACACCACTCAACACCCAGCTAACAGTCAGGAATACCAGGCTTAGTTCACAGCACTGAAAGAAGAAGCCACTCCCGCCATTCGTCATAAATTTGaacaac                             >  1:1826577/2‑150 (MQ=255)
                                                                                                                         aagagTAAGGCTGAACCTTCATGTTCAACCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGaacaac                             >  1:843121/2‑150 (MQ=39)
                                                                                                                          agagTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCATACCATGCACCTCTCACCAACCAGCTAACCGTCAGGAATACCAGCCTTAGTACACAGCACAGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAacaaca                            <  2:355235/150‑1 (MQ=255)
                                                                                                                           gagTAAGGCTGAACCTTCATGTTCAACCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAacac                           >  1:1270360/1‑150 (MQ=255)
                                                                                                                             ttAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTAAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACAACTCCAGCCATTCGTCATCAATTTGAACAACACaa                         >  1:520806/2‑150 (MQ=255)
                                                                                                                             gTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACaa                         >  1:2064638/1‑150 (MQ=255)
                                                                                                                             gTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACAGAAAGTAGAAACCACTCCAGCCATTCGGCATCAATTTTAACAACACaa                         >  2:384713/1‑150 (MQ=255)
                                                                                                                             gTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCAGCTAACAATTCAACACTCAACAACCAGATAACAGTCAGGAATACCAGGCTTAGTAAACAGAAATGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAGCACaa                         >  2:522093/1‑150 (MQ=255)
                                                                                                                                 ggTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACt                     >  2:583467/1‑150 (MQ=255)
                                                                                                                                    tGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGAc                  <  2:221354/150‑1 (MQ=255)
                                                                                                                                     gAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGAcc                 <  1:23831/150‑1 (MQ=255)
                                                                                                                                           gAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTACCCGTCAGGAATACCAGGCTTAGTACACAGCACTCAAAGTAGAAACCACTCCAGACAATCGTCATCAATTAGAACAACACAATACTGACCCACAtt           >  2:919134/1‑150 (MQ=255)
                                                                                                                                               gTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGACCCACATTCCCg       <  1:1436008/150‑1 (MQ=255)
                                                                                                                                                ttCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGACCCACATTCCCGt      <  1:448579/150‑1 (MQ=255)
                                                                                                                                                ttCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGACCCACATTCCCGt      <  2:1098322/150‑1 (MQ=255)
                                                                                                                                                    gCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACACCACTGAAAGTACAAACCACTCCAGCCATTCTTCATCAATTTGAACAACACAATACTGACACACATTCCCGTAATc  >  1:1349977/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGACCCACATTCCCGTAATC  >  NC_000913/3663416‑3663713

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: