Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913941,7130TG51.9% 87.4 / 127.1 77V252G (GTG→GGG) dmsAdimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base T (6/31);  new base G (40/0);  total (46/31)
Fisher's exact test for biased strand distribution p-value = 7.25e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-12

ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTG  >  NC_000913/941569‑941858
                                                                                                                                                |                                                                                                                                                  
actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACt                                                                                                                                                <  2:778153/149‑1 (MQ=255)
  tcctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGAGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTTGGGTg                                                                                                                                                   >  2:1018336/1‑144 (MQ=255)
   cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTa                                                                                                                                           <  2:717818/151‑1 (MQ=255)
   cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACt                                                                                                                                            >  2:454844/1‑150 (MQ=255)
      ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCt                                                                                                                                       <  1:560959/151‑1 (MQ=255)
       cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCtt                                                                                                                                      <  1:100354/151‑1 (MQ=255)
       cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTA‑TCtt                                                                                                                                      >  1:184169/1‑151 (MQ=255)
        gcgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTA‑TCtt                                                                                                                                      >  2:515456/1‑150 (MQ=255)
           cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAAc                                                                                                                                  <  2:880818/151‑1 (MQ=255)
            aaaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTA‑TATTTaa                                                                                                                                   >  1:225733/1‑149 (MQ=255)
            aaaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGAGTTTCTA‑TATTGaa                                                                                                                                   >  1:879024/1‑149 (MQ=255)
             aaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTA‑TCTTTaa                                                                                                                                   >  2:232509/1‑148 (MQ=255)
              aTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAgg                                                                                                                               <  1:861048/151‑1 (MQ=255)
              aTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAgg                                                                                                                               <  1:384535/151‑1 (MQ=255)
              aTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAgg                                                                                                                               <  1:32831/151‑1 (MQ=255)
              aTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTA‑TCTTTAACAgg                                                                                                                               >  2:692182/1‑151 (MQ=255)
                  cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTATTT‑TCTTTAAAAAGCAc                                                                                                                            >  1:968778/1‑150 (MQ=255)
                    gAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTTGGGGGGCTTATTATTCTGGAAAAGGccgcc                                                                                                                          >  2:559821/1‑147 (MQ=255)
                    gAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTA‑TCTTGAAAAGGCACGcc                                                                                                                         >  2:34293/1‑151 (MQ=255)
                       ggTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGa                                                                                                                      <  1:389765/151‑1 (MQ=255)
                           tGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAAt                                                                                                                   <  1:684244/150‑1 (MQ=255)
                                tataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTA‑TCTTTAAAAAGCAAGCCAGAAAACTaa                                                                                                               >  1:688054/1‑149 (MQ=255)
                                tataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGACTTAACA‑ACTTGAAAAGGCACGCCAGaaaattaa                                                                                                               >  2:844318/1‑144 (MQ=255)
                                 ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGc                                                                                                            <  1:779478/151‑1 (MQ=255)
                                 ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCCTTCTA‑TCTTGAAAAGGCACGCCAGAAATCTaa                                                                                                               >  2:748437/1‑148 (MQ=255)
                                    ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTA‑TCTTTAAAAGGCACGCCAGAAAACTAATGccc                                                                                                          >  2:271640/1‑150 (MQ=255)
                                        cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATg                                                                                                     <  1:845262/151‑1 (MQ=255)
                                         ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGc                                                                                                        >  2:281195/1‑147 (MQ=255)
                                          gcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGat                                                                                                   <  2:474970/151‑1 (MQ=255)
                                           cggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGATTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTA‑TCTTTAACAGGCACGCCAGAAAACTAATGCCCGcct                                                                                                      >  1:199678/1‑145 (MQ=255)
                                            ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTTGGGGGAACTTCTA‑TCTTGAACAGGGACGCCAGAAATCTAATGCCCGCacg                                                                                                     >  2:680509/1‑145 (MQ=255)
                                            ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatca                                                                                                 <  1:629678/151‑1 (MQ=255)
                                            ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatca                                                                                                 <  2:960636/151‑1 (MQ=255)
                                               tGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTA‑TTTTTAAAAAGCAAGCCAGAAAATTAATGCCCGCCAGatc                                                                                                  >  1:567526/1‑147 (MQ=255)
                                                gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTA‑TTTTTAAAAGGCACGCCAGAAAACTAAAGCCCGCCTGat                                                                                                   >  1:30881/1‑145 (MQ=255)
                                                 ggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatcatc                                                                                            <  1:393481/151‑1 (MQ=255)
                                                    aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTA‑TATTGAACAAGGACGCCAGAAATCTAATGCCCGCCTGAACATCATCg                                                                                           >  2:647398/1‑149 (MQ=255)
                                                        ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcg                                                                                      <  1:851655/150‑1 (MQ=255)
                                                          cAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTA‑TCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcgc                                                                                     >  1:867589/1‑149 (MQ=255)
                                                                 ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGTGAACTACAA‑TCTTGTAAAGGAACGCCAGAAAATTTATGCCCGCATTATTCACATCCATCCgcgccttaacc                                                                            >  2:528860/1‑144 (MQ=255)
                                                                   gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGac                                                                          <  1:859355/151‑1 (MQ=255)
                                                                     cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGaca                                                                         >  2:870178/1‑150 (MQ=255)
                                                                         tatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTA‑TCTTTAACAGGCACGCCAGAAAACTAATGGCCGCATGATCATCATCGATACGCGCTCTACCGACACCgg                                                                     >  1:712819/1‑150 (MQ=255)
                                                                          atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGt                                                                    <  2:500460/150‑1 (MQ=255)
                                                                            cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTA‑TCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACCCCGGTGcc                                                                 >  1:615185/1‑151 (MQ=255)
                                                                            cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGcc                                                                 >  2:405457/1‑151 (MQ=255)
                                                                             gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCCACAACGGGGcc                                                                 >  2:608347/1‑150 (MQ=255)
                                                                             gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTA‑TCTTGAACAAGCACGCCAGAAAACTAATGCCCCCATGATCATCATCGATCCGCGCTATACCGACACCGGTGcc                                                                 >  2:424095/1‑150 (MQ=255)
                                                                              aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTA‑TTTTAAACAGGCACGACAGAAAACTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCGACCCCCGTg                                                                   >  1:808818/1‑147 (MQ=255)
                                                                               aaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTTCTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGGTCCGCGCTATACCGACACCGGTGCCggg                                                              >  1:217349/1‑151 (MQ=255)
                                                                                          tGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATg                                                    >  2:252322/1‑150 (MQ=255)
                                                                                          tGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTA‑TCTTTAAAAGGCACGCCAGAAAACTAATGGCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGgaga                                                      >  1:617877/1‑148 (MQ=255)
                                                                                               gTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTgg                                               >  2:569568/1‑150 (MQ=255)
                                                                                                tGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTA‑TCTTGAACAGGCACGCCCGAAATCTAATGCCCGCATGATCATCATCGGTCCGCGCTATACCGACACCCGTGCCGGGCGCGGAGATGAGTgg                                               >  1:471593/1‑149 (MQ=255)
                                                                                                tGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTA‑TCTTGAACAGGCACGCCCGAAATCTAATGCCCCCATGATCATCATCGATCCGCGCTATACCCACACCCGTGCCGGGCGCGAAGATGAGTgg                                               >  2:951612/1‑149 (MQ=255)
                                                                                                 gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGAt                                             <  1:271619/150‑1 (MQ=255)
                                                                                                    gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATccc                                          <  2:175509/150‑1 (MQ=255)
                                                                                                        ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAtt                                      <  2:650773/150‑1 (MQ=255)
                                                                                                        ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTc                                     <  1:515413/151‑1 (MQ=255)
                                                                                                         gTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTA‑TCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGGAGATGAGGGGATCCCTAtt                                      >  2:531435/1‑149 (MQ=255)
                                                                                                           aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAGTCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGt                                   <  1:68451/150‑1 (MQ=255)
                                                                                                           aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTc                                  <  1:570250/151‑1 (MQ=255)
                                                                                                            ataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTA‑TTTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCCATACCCACACCGGGGCCGGGCGCGAAGATGAGTGGGTCCCCATTTGTcc                                 >  1:49933/1‑151 (MQ=255)
                                                                                                            ataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATTATCATCGATCCGCGCTATACCGACACCGGTGGCGGGCGCGAAGATGAGTGGATCCCTATTCGTcc                                 >  2:861778/1‑151 (MQ=255)
                                                                                                              aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTAATC‑TTTTGAACAAGCACGCCCCAAATTTAATGCCCGCATGGTCATCATTGATCCCCGCTATACCGTCACCGGGGCCGGGCGCGCAGAAGAGTGGAACCCTAAtcctcc                                 >  1:615809/1‑149 (MQ=255)
                                                                                                               aCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTACTA‑TCTTGAACAAGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAAACCGACACCCGGGCCGGGCGGGAAGATGAGGGGATCCCTATTCGTCCggg                              >  1:279837/1‑151 (MQ=255)
                                                                                                               aCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCggg                              <  2:968877/151‑1 (MQ=255)
                                                                                                               aCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCCTCATCGATCCGCGCTATACCGGCACCGGGGCCGGGCGCGGAGATGAGTGGGTCCCCATTCGTCCggg                              >  1:119032/1‑151 (MQ=255)
                                                                                                                    ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAAc                           <  2:223200/149‑1 (MQ=255)
                                                                                                                    ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAg                         <  1:1004641/151‑1 (MQ=255)
                                                                                                                    ggCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCCCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCCATTCGTCCggg                              >  2:941231/1‑146 (MQ=255)
                                                                                                                    ggCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGaa                            >  2:1004010/1‑148 (MQ=255)
                                                                                                                      cGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGAt                       <  2:179018/151‑1 (MQ=255)
                                                                                                                          aCGCGAATGAGTGGCGGTGGGGGGGCTTACTA‑TCTTGAACACGCACGCCCGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCCGTGCCGGGGGCGAAGATGAGGGGATCCCTATTCGTCCGGGGACa                          >  1:212590/1‑144 (MQ=255)
                                                                                                                               aaTGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTg              <  1:307937/151‑1 (MQ=255)
                                                                                                                                          gTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCt   <  2:37788/151‑1 (MQ=255)
                                                                                                                                           tGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTg  >  1:322589/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTA‑TCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTG  >  NC_000913/941569‑941858

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