Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913941,7130TG45.9% 117.6 / 44.0 85V252G (GTG→GGG) dmsAdimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base T (6/40);  new base G (39/0);  total (45/40)
Fisher's exact test for biased strand distribution p-value = 3.26e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGA‑GTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941566‑941862
                                                                                                                                                   |                                                                                                                                                      
actactCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAACCGCGAATGAGTGGCGGTGGGGTGAc                                                                                                                                                     <  3:489397/151‑1 (MQ=255)
   actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACt                                                                                                                                                    <  3:1601333/149‑1 (MQ=255)
      cctccGCGCAAATTGCGTAAGGTTTGAACTATCCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTAc                                                                                                                                                 <  2:1008751/149‑1 (MQ=255)
      cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACt                                                                                                                                                <  3:1434811/150‑1 (MQ=255)
      cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTAc                                                                                                                                                 >  2:1427244/1‑149 (MQ=255)
          cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCtt                                                                                                                                           <  2:718322/151‑1 (MQ=255)
             gcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGaa                                                                                                                                        <  3:1060203/151‑1 (MQ=255)
             gcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATCtt                                                                                                                                           >  3:680530/1‑148 (MQ=255)
              cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAc                                                                                                                                       <  3:422147/151‑1 (MQ=255)
              cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCtttatatttt                                                                                                                                           >  3:195612/1‑139 (MQ=255)
               aaaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTaa                                                                                                                                        >  2:1377047/1‑149 (MQ=255)
                aaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAACa                                                                                                                                      >  2:813190/1‑150 (MQ=255)
                 aTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAgg                                                                                                                                    <  2:491906/151‑1 (MQ=255)
                   tGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGc                                                                                                                                   <  2:1361443/150‑1 (MQ=255)
                    gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCAc                                                                                                                                 <  3:23116/151‑1 (MQ=255)
                       gAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTAATATTTTTAAAAAGCACGc                                                                                                                               >  3:524682/1‑150 (MQ=255)
                           gTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGaa                                                                                                                          <  2:1590473/151‑1 (MQ=255)
                            tttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTGAAAAAGCACGGCCGaa                                                                                                                          >  3:30100/1‑150 (MQ=255)
                             ttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGaaa                                                                                                                         <  2:1476860/150‑1 (MQ=255)
                             ttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGaaa                                                                                                                         <  3:1059845/150‑1 (MQ=255)
                               gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATc                                                                                                                       <  3:1480179/150‑1 (MQ=255)
                                aaCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTa                                                                                                                     <  2:673781/151‑1 (MQ=255)
                                aaCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTa                                                                                                                     <  3:1107301/151‑1 (MQ=255)
                                 aCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAAAAAGCACGCCCGaaa                                                                                                                         >  2:924436/1‑146 (MQ=255)
                                    ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATg                                                                                                                  <  3:793784/150‑1 (MQ=255)
                                    ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGTGAACTACTATCTTTAACAAGCACGCCCGAAATCTAATGc                                                                                                                 >  2:218271/1‑151 (MQ=255)
                                           cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACCCCCGAAATCTAATGcccccc                                                                                                            >  2:863415/1‑146 (MQ=255)
                                           cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCCGAAATTTAAAGCCCGCacg                                                                                                          >  3:678379/1‑149 (MQ=255)
                                            ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGa                                                                                                         <  3:764872/151‑1 (MQ=255)
                                            ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATTTTGAACAAGCACGCCAGAAATATAATGCCCGc                                                                                                             >  2:2359/1‑147 (MQ=255)
                                               ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAATTTCTATTTTTGAAACGCACGGCCGAAAAATTATGccccgatg                                                                                                          >  3:1005043/1‑142 (MQ=255)
                                               ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAACTTCTATCTTTAACAGGCACGCCACAAAACTAATGCCCGCCTGatc                                                                                                       >  4:9364/1‑150 (MQ=255)
                                                gcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcat                                                                                                     <  3:1340845/151‑1 (MQ=255)
                                                  tGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTTGGGGGAGTTACGATCTTTAAAAAGCACGCCAGAAAACTAAAGCCCCCATGatcat                                                                                                     >  2:566974/1‑149 (MQ=255)
                                                    ggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatcatc                                                                                                 <  2:1031481/151‑1 (MQ=255)
                                                     gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTTAAAAGGCACGCCAGAAAACTAATTCCCGCATGAACATCATCg                                                                                                >  2:938165/1‑151 (MQ=255)
                                                         aTGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATcc                                                                                            <  2:711300/151‑1 (MQ=255)
                                                          tGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGGACTCCAGAAAACTAATGCCCGCATGGTCATCATCGATcc                                                                                            >  3:1235398/1‑150 (MQ=255)
                                                                  gCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCCTACTATATTGAACAGGCGCACCAGAAAACTAATGCCCCCACGAGGATCCTCCATCCCCGCtat                                                                                     >  2:25070/1‑149 (MQ=255)
                                                                     cGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACCTACTCTCTTTAAAAGGCACCCCCGGAAATTAATGCCCGCCAGATAATCCTTGATTCCCCCTAAAcc                                                                                  >  2:252730/1‑149 (MQ=255)
                                                                      gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTACTATTTTGAACAGGGACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCcata                                                                                    >  3:1502240/1‑146 (MQ=255)
                                                                       tCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCCTCGATCCGCGCTATACCGac                                                                               >  2:146473/1‑150 (MQ=255)
                                                                        cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATCTTGAAAAAGGAAGCCCGAAAACTAATGCCCGCCTGATCATCATCGAACCGCGCCATACCGa                                                                                >  3:296501/1‑148 (MQ=255)
                                                                         ggATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACAcc                                                                            <  2:795170/151‑1 (MQ=255)
                                                                          gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGAATTACTTTTTTGAAAAAGCACGCCACAAAACTAATACCACCATTATCATAATCGATCGGCGCTCTACCGACAACg                                                                           >  2:58057/1‑151 (MQ=255)
                                                                          gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAATTACTATCTTTAACAAGCACGCCAGAAATCTAATGCCCCCCTGATCATCATCGATCCCCGCTCTACCCACAcc                                                                            >  2:1130045/1‑150 (MQ=255)
                                                                               cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGTTTAATATATTTAAAAAGGAAGCCAGAAAATTAATGCCCGCCTGCTCATCATCTGGCCGCGCGATAACCACACCgggg                                                                        >  3:1210849/1‑147 (MQ=255)
                                                                                gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATTGATCCGCGCCATACCCACACCCGGGcc                                                                      >  3:119788/1‑150 (MQ=255)
                                                                                gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGCACGCCAGAAAACTAATGCCCGCCTGATCATCATCGATCCGCGCCATACCGACACCGGGGcc                                                                      >  2:1145415/1‑150 (MQ=255)
                                                                                 aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCgg                                                                    <  3:1203401/151‑1 (MQ=255)
                                                                                     cAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgcgc                                                                <  2:1248427/151‑1 (MQ=255)
                                                                                        tAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAAGCACGCCAGAAATCTAATGCCCCCATGATCATCATCGATCCGCGCTATACCGACACCCGTGCCGGGcgcg                                                               >  2:1629756/1‑149 (MQ=255)
                                                                                         aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGaa                                                             <  3:402115/150‑1 (MQ=255)
                                                                                         aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGaa                                                             <  3:463655/150‑1 (MQ=255)
                                                                                          aGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCAGAAAACTAAATCCCCCCCGATCCTCATCCATCCGCGCTATACCCACACCCGGGGCggg                                                                   >  2:1621476/1‑143 (MQ=255)
                                                                                            cTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGGAggt                                                          >  3:1011624/1‑148 (MQ=255)
                                                                                            cTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTATTATCTTTAACAGGCACGCCAGAAAAATAATGCCCGCCTGATAATCATCTAACCGCGCCATACCCACACCGGTGCCGGGCGCGGAGATg                                                         >  3:122173/1‑151 (MQ=255)
                                                                                             tGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTTATGCCCGCATGATCATCATCGATTCGCGCTATACCCACACCGGTGGCGGGggcg                                                               >  1:5434/1‑144 (MQ=255)
                                                                                             tGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAAACTAATGCCCGCCTGATCATCATCGATCCGCGCTATAACGACACCGGGGCCGGGGGCGAAGATg                                                         >  2:1174976/1‑150 (MQ=255)
                                                                                                     cTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCGACACCGGGGGCGGGGGCGAAGATGA‑GTGGGTcc                                               >  2:186193/1‑151 (MQ=255)
                                                                                                      tGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAAGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCCATAACGACACCCGTGCCGGGCGCGAAGAAGA‑GTGGATccc                                              >  2:1510785/1‑151 (MQ=255)
                                                                                                          tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGA‑GTGGATCCCTAt                                           <  3:1630150/150‑1 (MQ=255)
                                                                                                           ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGA‑GTGGATCCCTAtt                                          <  2:818822/150‑1 (MQ=255)
                                                                                                            gTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGA‑GTGGATCCCTATTCg                                        <  3:923179/151‑1 (MQ=255)
                                                                                                             taataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTAATATCTTTAACAAGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCCATCCGCGCTATACCCACACCCGTGGCGGGCGCGGAGATGA‑GTGGGTCCCTAtt                                          >  3:694800/1‑148 (MQ=255)
                                                                                                              aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGA‑GTGGATCCCTATTCGt                                       <  2:1145218/150‑1 (MQ=255)
                                                                                                              aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGAATTAATATCTTTTACAGGCACGCCACAAAACTAATGCCCCCATGATCATCATCTCTCCCCGCTATACCCCCCCCCGGGCCCGGGGCGAAGATGA‑GTGGGTCCCCAttcttc                                      >  3:26308/1‑151 (MQ=255)
                                                                                                               ataaCCCTGGCGAAACGCGAATGAGTGGCGGATGGGGGTTTTACTATCTTGAAAAGACACCCCAGCAAACTAATGTCCGCATGACCATCTTCGCAACGCGCCATACCGACACCCGGGCCGGGCGGGGAGATGA‑GTGGAGccccctttctcc                                     >  3:1390699/1‑142 (MQ=255)
                                                                                                                taaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAATTAATATCTTGAACAGGCCCGCCAGAAAACTAATGCCCGCATGATCATCATCGATACGCGCCATAACGACACCGGGGCCGGGGGCGGAGATGA‑GGGGATCCCCATTCGTcc                                     >  3:1532378/1‑150 (MQ=255)
                                                                                                                taaCCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACCATTTTTAACAGGCACGCCCGAAAACTAATGCCCCCCTGATCATCATCACTCCGCCCCATACCCCCACCCGCGCCGGGCGCGCAGGAGA‑GTGGGTGCCTATTCTTcc                                     >  3:1072565/1‑150 (MQ=255)
                                                                                                                 aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACCATCCTTAAAAAGCAACCCAAAAATCTAAAGCCCGCAAGATCAATAATGATCCGCCCCATACCCACACCGGTGCCCGGGGCGAAAATGAGGGGGATcccatttggtc                                      >  2:1514681/1‑141 (MQ=255)
                                                                                                                   cccTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATTTTGAACAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCCACACCCGTGCCGGGCGCGGAGATGA‑GTGGGTCCCTATTCGTCCggg                                  >  3:540443/1‑150 (MQ=255)
                                                                                                                       ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGA‑GTGGATCCCTATTCGTCCGGGAAc                               <  3:334648/149‑1 (MQ=255)
                                                                                                                       ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGA‑GTGGATCCCTATTCGTCCGGGAAc                               <  3:750139/149‑1 (MQ=255)
                                                                                                                       ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGA‑GTGGATCCCTATTCGTCCGGGAACAg                             <  3:84175/151‑1 (MQ=255)
                                                                                                                          gAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGCACGCCAGAAATCTAATGCCCCCCTGATCATCATCGATCCGCGCCATACCCACACCGGGGCCGGGCGCGAAAATGA‑GTGGATCCCTTTTCGTCCcgg                                  >  3:1425225/1‑143 (MQ=255)
                                                                                                                              cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGA‑GTGGATCCCTATTCGTCCGGGAACAGATGCCGc                      <  2:497705/151‑1 (MQ=255)
                                                                                                                                     gAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGA‑GTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGtt               <  3:814307/151‑1 (MQ=255)
                                                                                                                                     gAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGA‑GTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGtt               <  2:645220/151‑1 (MQ=255)
                                                                                                                                      aGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGA‑GTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTa              <  2:620613/151‑1 (MQ=255)
                                                                                                                                             gTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGA‑GTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCt       <  2:266862/151‑1 (MQ=255)
                                                                                                                                              tGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGA‑GTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTg      >  2:187830/1‑151 (MQ=255)
                                                                                                                                                gggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGA‑GTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGc    <  3:1353717/151‑1 (MQ=255)
                                                                                                                                                 ggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGA‑GTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCg   <  2:213730/151‑1 (MQ=255)
                                                                                                                                                 ggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGA‑GTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCg   >  2:200370/1‑151 (MQ=255)
                                                                                                                                                   tGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGA‑GTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGt  >  3:438049/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
ACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGA‑GTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941566‑941862

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: ‑