Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913941,7150AG43.5% 130.0 / 11.6 86T253A (ACT→GCT) dmsAdimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base A (12/36);  new base G (37/0);  total (50/36)
Fisher's exact test for biased strand distribution p-value = 5.69e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTAC  >  NC_000913/941568‑941864
                                                                                                                                                   |                                                                                                                                                      
tactCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCtt                                                                                                                                                     >  2:1042525/1‑151 (MQ=255)
 actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACt                                                                                                                                                      <  1:422524/149‑1 (MQ=255)
 actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTa                                                                                                                                                    <  2:397398/151‑1 (MQ=255)
   tcctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTa                                                                                                                                                    >  1:233787/1‑149 (MQ=255)
   tcctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCt                                                                                                                                                  >  2:664063/1‑151 (MQ=255)
   tcctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCATACt                                                                                                                                                  >  1:446864/1‑151 (MQ=255)
    cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTAc                                                                                                                                                   <  2:14501/149‑1 (MQ=255)
    cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACt                                                                                                                                                  >  2:757403/1‑150 (MQ=255)
      tccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTa                                                                                                                                                 >  1:243693/1‑149 (MQ=255)
       ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAt                                                                                                                                                <  1:747371/149‑1 (MQ=255)
        cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTAATTTCtt                                                                                                                                             >  1:697893/1‑151 (MQ=255)
         gcgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTg                                                                                                                                            <  1:1357381/151‑1 (MQ=255)
           gcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGaa                                                                                                                                          <  1:1276126/151‑1 (MQ=255)
           gcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTTTCtt                                                                                                                                             >  2:416526/1‑148 (MQ=255)
             aaaTTGCGGAAGGTTTGAACTATCCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACa                                                                                                                                        <  1:1313000/151‑1 (MQ=255)
              aaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACa                                                                                                                                        <  2:21820/150‑1 (MQ=255)
               aTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAgg                                                                                                                                      <  1:885169/151‑1 (MQ=255)
                 tGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGGCTTACTTTCTTTAACAgg                                                                                                                                      >  1:135485/1‑149 (MQ=255)
                   cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACg                                                                                                                                  <  2:601672/151‑1 (MQ=255)
                   cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGATTTTCTTACTTGAACAggta                                                                                                                                    >  1:230793/1‑147 (MQ=255)
                     gAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGcc                                                                                                                                <  1:254194/151‑1 (MQ=255)
                      aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTACTCTTTTTTACAAGCACGcc                                                                                                                                >  2:984794/1‑150 (MQ=255)
                         gTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCCGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTAATATCTTTCAAAAGCACGCCAg                                                                                                                              >  1:498031/1‑149 (MQ=255)
                          tttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAAAAGGGACGcc                                                                                                                                >  1:152849/1‑146 (MQ=255)
                                cTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTTCTATTTTTTAAAAGCACGCCCGaaaaataa                                                                                                                      >  2:355937/1‑145 (MQ=255)
                                        aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTTCTATATTTAAAAAGCAAGCCCGAAATCTAATAccc                                                                                                                >  2:476600/1‑147 (MQ=255)
                                          ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTACTATCTTTAAAAGGCACGCCCGAAAACTTAAGCCCGCacg                                                                                                           >  2:1314549/1‑148 (MQ=255)
                                           gcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTTTCTTTGACAGGAACTGCCGAAAACTTATGccc                                                                                                                >  1:712384/1‑144 (MQ=255)
                                             ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGGACGCCAGAAATCTAATGGCCGCATGat                                                                                                         >  1:1243986/1‑149 (MQ=255)
                                             ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATTTTGAACAGGCACGCCAGAAATCTAATGccccc                                                                                                              >  1:440208/1‑142 (MQ=255)
                                                tGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACTTTCTATTTTGGACAGGCACGCCAGAAAAATAAAGCCCGCAAGatcat                                                                                                      >  2:573033/1‑149 (MQ=255)
                                                 gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatca                                                                                                    <  2:256429/150‑1 (MQ=255)
                                                  ggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatcatc                                                                                                  <  1:1320403/151‑1 (MQ=255)
                                                     aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGAt                                                                                               <  2:422760/151‑1 (MQ=255)
                                                     aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATCTTTAAAAGGCACGCCAGAAAACTAATGCCCGCATGatca                                                                                                       >  2:529353/1‑143 (MQ=255)
                                                      gATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATc                                                                                              <  1:644648/151‑1 (MQ=255)
                                                      gATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCATGCCAGAAATCTAATGCCCGCATGATCATCATCTAt                                                                                               >  2:1006469/1‑150 (MQ=255)
                                                        tGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTACTATTTTTAACAGGCACCCCCGAAAACTAATGCCCGCATGATCATCATCGATcc                                                                                             >  2:1182462/1‑150 (MQ=255)
                                                         ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcgc                                                                                           <  2:520698/151‑1 (MQ=255)
                                                           cgaCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcgcg                                                                                          <  2:231327/148‑1 (MQ=255)
                                                            aaCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTCTTAACTTTGAAACGAACCCCAGAAAACTAAAGCCCCCAATatcatcaaccttccctccc                                                                                        >  1:1288835/1‑139 (MQ=255)
                                                             aCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAAAAGGGACGCCCGAAAACTAATGCCCCCATGATCATCATCGATCCccgc                                                                                         >  1:1447582/1‑149 (MQ=255)
                                                               aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATTTTTAACAGGCACGCCAGAAAACTAATGCCCGCCCGATCATCatccatcc                                                                                             >  2:825197/1‑143 (MQ=255)
                                                                 cccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATTTTGAAAAAGCACGCCCGAAATTTAATTGCCCCCAGATAATCATCGAACCCCGCCATAcc                                                                                   >  2:51218/1‑151 (MQ=255)
                                                                 cccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTTATATCTTTAACAAGCACGCCCGAAATATTATAACCCCAAGgtcatcatc                                                                                                  >  2:1201301/1‑136 (MQ=255)
                                                                       ggATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATCTTGAAAAAGCACGCCAGAAATATAATGGCCGCCAGATCATCATTGATCCGCGCTATACCGacac                                                                              >  2:840992/1‑150 (MQ=255)
                                                                         atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCgg                                                                           <  1:757418/151‑1 (MQ=255)
                                                                         atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAGGCACGCCCGAAAACTAATTCCCCCATGATCATCATTGATACGCGCTATAACcacaac                                                                             >  2:4904/1‑147 (MQ=255)
                                                                           atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGt                                                                          <  1:1116865/150‑1 (MQ=255)
                                                                                 aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAAACTAATGCCCCCATGATCATCATCGATCCGCGCTATACCCACACCCGTGCCggg                                                                    >  1:95633/1‑150 (MQ=255)
                                                                                  aCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgcg                                                                  <  1:526574/151‑1 (MQ=255)
                                                                                  aCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTATTATCTTTAAAAAGCACGCCAGAAAACTAAATCCCCCCTGATCATCATCGATCCGCGCTATAACGCCACCCGTGCCggggg                                                                  >  2:900645/1‑149 (MQ=255)
                                                                                     gTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGTGAGTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCACGATCATCATCGATACGCGCTATACCCGCAAccgggccggg                                                                    >  2:887488/1‑146 (MQ=255)
                                                                                       aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAGTTACTATCTTGAAAAAGGACGCCAGAAAACTAATGCCCGCATGGTCATCAACGGTTCGCGCTATACCGGCACCCGGGCCGGGcgcg                                                                >  2:839648/1‑148 (MQ=255)
                                                                                         gCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGAt                                                           <  1:969856/151‑1 (MQ=255)
                                                                                          cTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTACTATCTTGAACAGGCACCCCAGAAAACTAATGCCCGCATGATCATCATCGATTCGCGCTATACCGACACCGGTGCCGGGcgcg                                                                >  2:829083/1‑145 (MQ=255)
                                                                                                gTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGa                                                    <  1:416533/151‑1 (MQ=255)
                                                                                                gTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATCTTTAACAAGCACGCCAGAAAACTAATGCCCGCATGATTATCCTCGAGCCCCGGTCTACCGCCACCCGGGGCGGGCGCGGAGATGAgggg                                                     >  2:751162/1‑150 (MQ=255)
                                                                                                    tGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTCGGGGGAGTTACTATCTTTCACAGGCACCCCTGAAATCCGAATCCCGCCTGCACCTCCTCGGAACCCGCCATACCGGCCCCGGGGCCCGGCGCGAAGAAGAGTggggccc                                                >  2:6435/1‑146 (MQ=255)
                                                                                                     gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCt                                               <  1:647969/151‑1 (MQ=255)
                                                                                                      tttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGTCTTATTATCTTGAACAGGCACGCCAGAAAACTAATGCCCCCATGAACATAATCGATCCGCGCTATACCGACACCGGGGGCGGGCGCGGAGATGAGGGGATccc                                                >  2:1356958/1‑149 (MQ=255)
                                                                                                         ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTc                                           <  2:1046464/151‑1 (MQ=255)
                                                                                                            aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGt                                         <  1:111716/150‑1 (MQ=255)
                                                                                                            aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTACTTTCTTTAACAGGCACGCCAGAAAAATAATGCCCGCAAGAACATCCTCAATCCGCGCCATACCCACACCCGGGCCGGGGGCGAAGATGAGGGGATccccctt                                            >  1:785475/1‑143 (MQ=255)
                                                                                                               aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTTAACAGGGCCGCCCGAAAACTAATGGCCGCCTGATCATCATCGATCCGCGCCATACCGACACCCGTGCCGGGGGCGAAGATGAGGGGGTCCCCATTCTTcccg                                     >  2:736719/1‑149 (MQ=255)
                                                                                                                  ccTGGCGAAACGCGAATGAGTGGCGGCGGGGGGGCTTATTTTCCTTAACAGGGACGCCAGAAAAATAATGCCCGCCCGATTATCATCGATTCGCGCTATAACCCCACCCGGGCCGGGGGCGAAGATGAGGGGATCCCTATTTGTCCgg                                     >  1:410651/1‑148 (MQ=255)
                                                                                                                     ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAg                               <  1:1159764/151‑1 (MQ=255)
                                                                                                                     ggCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAAAAGGCACGCCAGAAAACTAAAGCCCGCCTGATCATCATCGATCCGCGCTATACCGACACCCGGGCCGGGCGCGGAGATGAGGGGATCCCTATTCGTCCgggga                                  >  2:1368398/1‑146 (MQ=255)
                                                                                                                        gAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAAAAAGCACGCCAGAAATCTAATGCCCGCATGATCATTATCGATCCGCGCCATAACCCCACCCGGGGCGGGCGCGGAGATGAGTGGGTCCCTATTCGTCCCGGAACa                                >  2:895347/1‑147 (MQ=255)
                                                                                                                        gAAACGCGAATGAGTGGCGGTGGGGGGGATTTCTATCTTTAACAAGCACCCCAGAAATCTAATGCCCGCATGATCATCCTCGTTCCCCCCTATACCCACCCCCGTGCCGGGCGCGGAGATGAGGGGATCCCCCTTTGTCCCGGGACAGa                              >  2:968916/1‑149 (MQ=255)
                                                                                                                                 aTGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCCCGCCAGAAATCTTATGCCCGCATGATCATCATCGATCCGCGCTATACCCACACCGGGGCCGGGGGCGGAGATGAGTGGGTCCCTATTCGTCCGGGGACAGATGCCGCCCTg                    >  1:1416243/1‑150 (MQ=255)
                                                                                                                                  tGAGTGGCGGTGGGGGGGCTTACTATCTTGAAAAAGCACGCCAGAAATATAATGCCCCCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGGTCCCTATTCGTCCGGGGACAGATGCCGc                        >  1:808198/1‑145 (MQ=255)
                                                                                                                                   gAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGt                  <  2:1258027/150‑1 (MQ=255)
                                                                                                                                    aGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTa                <  2:243690/151‑1 (MQ=255)
                                                                                                                                         cGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGt           <  1:900105/151‑1 (MQ=255)
                                                                                                                                         cGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGt           <  2:250544/151‑1 (MQ=255)
                                                                                                                                         cGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGt           >  2:1460225/1‑151 (MQ=255)
                                                                                                                                            tGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTg        >  2:816721/1‑151 (MQ=255)
                                                                                                                                             ggggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTgg       <  1:1042550/151‑1 (MQ=255)
                                                                                                                                             ggggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTgg       <  2:233784/151‑1 (MQ=255)
                                                                                                                                             ggggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTg        <  2:464738/150‑1 (MQ=255)
                                                                                                                                              gggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGc      <  2:1443180/151‑1 (MQ=255)
                                                                                                                                               ggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCg     <  1:749630/151‑1 (MQ=255)
                                                                                                                                                 tGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGt    >  2:49414/1‑150 (MQ=255)
                                                                                                                                                 tGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGt    >  1:353209/1‑150 (MQ=255)
                                                                                                                                                  gACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTAc  <  1:984819/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
TACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTAC  >  NC_000913/941568‑941864

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