Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,663,558 G→T 6.3% intergenic (+31/+332) 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,5580GT 6.3% 461.4 / 16.1 153intergenic (+31/+332)mdtF/gadWanaerobic multidrug efflux transporter, ArcA‑regulated/transcriptional activator of gadA and gadBC; repressor of gadX
Reads supporting (aligned to +/- strand):  ref base G (75/68);  new base T (7/3);  total (82/71)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.99e-01

Only 100 of 153 total aligned reads displayed.

GTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGACCCACATTC  >  NC_000913/3663414‑3663705
                                                                                                                                                |                                                                                                                                                   
ttaggtaCCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTTGCAATTTAATTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCATTGCCCGCTATAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCa                                                                                                                                                <  1:174954/149‑1 (MQ=255)
  aggtaCCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGGCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGc                                                                                                                                              <  2:224716/150‑1 (MQ=255)
     taCCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCtt                                                                                                                                           >  2:183707/1‑150 (MQ=255)
       ccGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTAc                                                                                                                                         <  2:2027833/150‑1 (MQ=255)
       ccGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTAc                                                                                                                                         >  2:391610/1‑150 (MQ=255)
        cGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGTCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTAct                                                                                                                                        <  1:565781/150‑1 (MQ=255)
        cGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCTGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGCTGAACCTTCATGTTCAACCTTAct                                                                                                                                        <  2:2057470/150‑1 (MQ=255)
          gCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGGGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGCTGAACCTTCATGTTCAACCTTActct                                                                                                                                      <  1:1951150/150‑1 (MQ=39)
           cGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCtt                                                                                                                                     >  2:1766768/1‑150 (MQ=255)
              gATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTcc                                                                                                                                  <  1:867269/150‑1 (MQ=255)
                 gggCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAACAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGc                                                                                                                               <  2:11907/150‑1 (MQ=255)
                 gggCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGc                                                                                                                               <  1:2140853/150‑1 (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:483163/150‑1 (MQ=255)
                              gCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTcacca                                                                                                                  >  1:1028084/1‑150 (MQ=255)
                                cGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACt                                                                                                                >  2:588995/1‑150 (MQ=255)
                                   aaCAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTcaa                                                                                                             >  1:970670/1‑150 (MQ=255)
                                   aaCAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTcaa                                                                                                             <  2:411879/150‑1 (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)
                                         gCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGACCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCa                                                                                                       <  1:671320/150‑1 (MQ=255)
                                          cTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAg                                                                                                      <  1:1002793/150‑1 (MQ=255)
                                           tGGCAATTTAATTCGTTCAGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCATGCTAACCATTCACCACTCAACAACCAGc                                                                                                     >  2:1008444/1‑150 (MQ=255)
                                             gCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTa                                                                                                   >  2:1904320/1‑150 (MQ=255)
                                              cAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCCCCACTCAACAACCAGCTaa                                                                                                  >  1:1861340/1‑150 (MQ=255)
                                                   tACTTCGTTCCGGTCTTTTTCGGTGTAGTGGAACATCTCTTTGCCCGCGTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCGTTACTCTGTCCGGCTAACCATTCCCCACTCCACAGCCAGCTAACCGTc                                                                                             >  2:2072937/1‑150 (MQ=255)
                                                     cTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAg                                                                                           >  1:765612/1‑150 (MQ=255)
                                                      ttCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAgg                                                                                          <  1:410350/150‑1 (MQ=255)
                                                       tCGTTCCTGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGATGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGa                                                                                         <  1:2123520/150‑1 (MQ=255)
                                                          ttCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATa                                                                                      >  2:904902/1‑150 (MQ=255)
                                                            ccGGTATTTTTCGTTGTAGTGGAACATCTATTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATAcc                                                                                    <  2:373272/150‑1 (MQ=255)
                                                             cGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACAATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCa                                                                                   >  1:308158/1‑150 (MQ=255)
                                                              ggTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAg                                                                                  <  1:275603/150‑1 (MQ=255)
                                                               gTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAgg                                                                                 <  1:1442366/150‑1 (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                                                                         <  1:353919/150‑1 (MQ=255)
                                                                       cGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTac                                                                         >  2:2033897/1‑150 (MQ=255)
                                                                            tAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGc                                                                    >  1:2092359/1‑150 (MQ=255)
                                                                             aGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCa                                                                   >  1:493028/1‑150 (MQ=255)
                                                                                ggAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTg                                                                >  1:2189063/1‑150 (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:1044154/150‑1 (MQ=255)
                                                                                      tctcTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTa                                                          <  2:26862/150‑1 (MQ=255)
                                                                                       ctctTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAg                                                         <  2:765612/150‑1 (MQ=255)
                                                                                       ctctTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAg                                                         >  1:221354/1‑150 (MQ=255)
                                                                                       ctatttGCCCGCTTTAAAAAAGCGTTACGTGTAAATGAGATTAATGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAAGCATTCACCACTCAAAAACCAGCTAACCTTCAGGAATACCAGGCTTAGTACACAGCAGTGAAAGTAg                                                         <  1:341337/147‑1 (MQ=255)
                                                                                        tctTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGa                                                        >  1:1098322/1‑150 (MQ=255)
                                                                                         ctTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGaa                                                       <  2:272708/150‑1 (MQ=255)
                                                                                         ctTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGaa                                                       <  1:723498/150‑1 (MQ=255)
                                                                                         ctTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGaa                                                       <  2:1561412/150‑1 (MQ=255)
                                                                                          tttGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGaaa                                                      <  1:1637903/150‑1 (MQ=255)
                                                                                            tGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAAcc                                                    >  1:2085721/1‑150 (MQ=255)
                                                                                            tGCACGCTTTAAAAAATCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGGTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCATGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAAcc                                                    <  1:1567749/150‑1 (MQ=255)
                                                                                               ccGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACt                                                 <  1:183707/150‑1 (MQ=255)
                                                                                                cGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTc                                                >  2:1156138/1‑150 (MQ=255)
                                                                                                  cTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTACTACACAGCACTGAAAGTAGAAACCACTCCa                                              >  1:1696381/1‑150 (MQ=255)
                                                                                                    ttAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCAATCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGc                                            <  1:1164697/150‑1 (MQ=255)
                                                                                                       aaaaaGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCAt                                         >  2:1046929/1‑150 (MQ=255)
                                                                                                       aaaaaGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCAt                                         >  2:435801/1‑150 (MQ=255)
                                                                                                          aaGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCg                                      <  2:2099396/150‑1 (MQ=255)
                                                                                                              gTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACCGCACTGAAAGAAGAAACCACTCCAGCCATTCGtcat                                  >  2:1437737/1‑150 (MQ=255)
                                                                                                               tAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGtcatc                                 >  1:689035/1‑150 (MQ=255)
                                                                                                                 aCGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCaa                               >  2:1628039/1‑150 (MQ=255)
                                                                                                                   gtgtAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAAtt                             <  1:391610/150‑1 (MQ=255)
                                                                                                                   gtgtAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAAtt                             <  2:2087513/150‑1 (MQ=255)
                                                                                                                     gtAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTGGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTg                           <  2:133747/150‑1 (MQ=255)
                                                                                                                      tAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGa                          <  2:920844/150‑1 (MQ=255)
                                                                                                                      tAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGa                          <  2:1705294/150‑1 (MQ=255)
                                                                                                                      caggaaagagTAAGGCTGAACCTTCATGTTCAACCTTACTCTTTCCTGCTAACCATTCACCCCGCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACCCCGCACTGAAGGTAGCAACCACTCCAGCCATACGTCATCAATTTGa                          >  2:1583494/7‑150 (MQ=25)
                                                                                                                        ggaaagagTAAGGCTGATCCTTCATGTTCAACCTTACTCTTTCCTGCTAAACATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGcac                        >  1:1112102/5‑150 (MQ=38)
                                                                                                                        aaTGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGaac                        <  2:1171658/150‑1 (MQ=255)
                                                                                                                        aaTGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGaac                        <  2:1776077/150‑1 (MQ=255)
                                                                                                                        aaTGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGGCATCAATTTGaac                        >  1:1373710/1‑150 (MQ=255)
                                                                                                                          tGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCCTTCACCACTCAACACCCAGCTAACCGTCAGGAATACCAGGCTTAGTAAACAGCACGGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGaacaa                      >  2:1087724/1‑150 (MQ=255)
                                                                                                                          tGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGaacaa                      >  2:1623857/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)
                                                                                                                               gTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCAGCTAACAATTCAACACTCAACAACCAGATAACAGTCAGGAATACCAGGCTTAGTAAACAGAAATGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAGCACaa                 >  2:522093/1‑150 (MQ=255)
                                                                                                                                      tGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGAc          <  2:221354/150‑1 (MQ=255)
                                                                                                                                       gAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGAcc         <  1:23831/150‑1 (MQ=255)
                                                                                                                                         aCATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGACCca       <  1:2186172/150‑1 (MQ=255)
                                                                                                                                            tGAAGGTTAAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACAACTCCAGCCATTCGTCATCAATTTGAACAACACTATACTGTCCCACAt    >  2:345066/1‑150 (MQ=255)
                                                                                                                                             gAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTACCCGTCAGGAATACCAGGCTTAGTACACAGCACTCAAAGTAGAAACCACTCCAGACAATCGTCATCAATTAGAACAACACAATACTGACCCACAtt   >  2:919134/1‑150 (MQ=255)
                                                                                                                                              aaGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGACCCACATTc  >  1:2016282/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGACCCACATTC  >  NC_000913/3663414‑3663705

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: