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

Only 100 of 119 total aligned reads displayed.

CGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGC  >  NC_000913/941564‑941860
                                                                                                                                                     |                                                                                                                                                   
cGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTg                                                                                                                                                    <  3:165784/151‑1 (MQ=255)
 gACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGa                                                                                                                                                   <  2:106645/151‑1 (MQ=255)
     actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACtt                                                                                                                                                <  4:320962/150‑1 (MQ=255)
       tcctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCtt                                                                                                                                                >  4:1895155/1‑148 (MQ=255)
           ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAt                                                                                                                                           <  2:277061/149‑1 (MQ=255)
           ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCttt                                                                                                                                           >  4:527657/1‑147 (MQ=255)
              cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGa                                                                                                                                      <  4:740317/151‑1 (MQ=255)
               gcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCtt                                                                                                                                        >  3:337408/1‑148 (MQ=255)
                cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTaa                                                                                                                                     >  2:169871/1‑150 (MQ=255)
                cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTATATCTTTaa                                                                                                                                     >  4:624056/1‑150 (MQ=255)
                cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATTTTTaa                                                                                                                                     >  4:886489/1‑150 (MQ=255)
                 aaaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACa                                                                                                                                   <  1:1668585/151‑1 (MQ=255)
                  aaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTCTCTTGAACa                                                                                                                                   >  4:566595/1‑150 (MQ=255)
                  aaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTATTATTTTGaaaa                                                                                                                                   >  1:1508837/1‑148 (MQ=255)
                    ttGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTTCTATCTTTAACaag                                                                                                                                 >  4:1318652/1‑148 (MQ=255)
                     tGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGc                                                                                                                                <  4:850066/150‑1 (MQ=255)
                       cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACg                                                                                                                             <  1:1189833/151‑1 (MQ=255)
                          aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCa                                                                                                                          <  1:85976/151‑1 (MQ=255)
                             gTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGaa                                                                                                                       <  3:241916/151‑1 (MQ=255)
                                tGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAAt                                                                                                                     <  4:1586194/150‑1 (MQ=255)
                                   aCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTaa                                                                                                                 <  4:232922/151‑1 (MQ=255)
                                     tataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTATTTTTAACAAGCACGCCCGAAAACTaa                                                                                                                 >  4:360470/1‑149 (MQ=255)
                                      ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATg                                                                                                               <  1:729982/150‑1 (MQ=255)
                                        aCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATCTTGAAAAAGGCCGCCAAAAATCTaaagtcc                                                                                                            >  1:911050/1‑146 (MQ=255)
                                        aCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATTTTTAACAAGGACCCCCGAAAACTAAAGccc                                                                                                            >  4:581268/1‑151 (MQ=255)
                                        aCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGATTAATAACTTGAAAAGGCACGCCAGAAAaataattgcc                                                                                                            >  4:883832/1‑147 (MQ=255)
                                         ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTACTATCTTTAAAAAGCAACCCCGAAAACTaacacccc                                                                                                           >  1:698350/1‑145 (MQ=255)
                                           tACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAAACCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTTATATCTTGAAAAAGCACCCCAGaaaaacaaatgccgc                                                                                                          >  2:11030/1‑138 (MQ=255)
                                            aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGc                                                                                                          <  1:174747/149‑1 (MQ=255)
                                            aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATCTTGAAAAAGGACGCCAGAAAAATAATACCCg                                                                                                           >  3:148487/1‑148 (MQ=255)
                                             cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGc                                                                                                          <  4:1803713/148‑1 (MQ=255)
                                             cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCAt                                                                                                        <  1:92491/150‑1 (MQ=255)
                                             cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATg                                                                                                       <  2:13803/151‑1 (MQ=255)
                                                     gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatca                                                                                                <  1:1302561/150‑1 (MQ=255)
                                                     gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatc                                                                                                 <  4:750195/149‑1 (MQ=255)
                                                      ggCAGATGGCAACAGCCCGCCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatcatc                                                                                              <  4:1096953/151‑1 (MQ=255)
                                                        cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGa                                                                                            <  1:1468983/151‑1 (MQ=255)
                                                        cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatcatc                                                                                              >  4:369832/1‑149 (MQ=255)
                                                         aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGAt                                                                                           <  1:246773/151‑1 (MQ=255)
                                                            tGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAAGGCCGCCCGAAATCTAATGCCCGCATGATCATCATCGATcc                                                                                         >  2:119463/1‑150 (MQ=255)
                                                             ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcg                                                                                        <  4:1291849/150‑1 (MQ=255)
                                                             ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcg                                                                                        <  4:394702/150‑1 (MQ=255)
                                                             ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcg                                                                                        <  4:833214/150‑1 (MQ=255)
                                                                 aCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCta                                                                                   <  2:410655/151‑1 (MQ=255)
                                                                   aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCtata                                                                                 <  4:393921/151‑1 (MQ=255)
                                                                   aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCtata                                                                                 <  1:1216953/151‑1 (MQ=255)
                                                                   aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAAGCACGCCCGAAATCTTATGGCCGCCTGATCATCATCGATCCGCGCta                                                                                   >  4:951117/1‑149 (MQ=255)
                                                                      ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAAGCACGCCCGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAcc                                                                               >  1:10908/1‑150 (MQ=255)
                                                                      ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAATTACTATCTTTAAAAGGCCCCCCAGAAAACTAATGCCCGCATGAACATCATCGATCCGCGCTATAcc                                                                               >  4:1398092/1‑150 (MQ=255)
                                                                       cGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCCGAAAACTAATGCCCGCATGATCATCATCGATCCCCGCTATAcc                                                                               >  4:143972/1‑149 (MQ=255)
                                                                         tCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTTTCTTGAACAGGCACGCCAGAAAAATAAAGCCCGCATGATCATCATCGATCCGCGCTATACCGac                                                                            >  3:476669/1‑150 (MQ=255)
                                                                          cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCCTTCTTTTTTGAACAAGCACGCCAGAAAACTAAAGACCGCCTGAACAGCATCGAGACGGGCCATAACGaca                                                                           >  4:1123323/1‑150 (MQ=255)
                                                                           ggATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTAATATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCCACAcc                                                                         >  4:1216097/1‑151 (MQ=255)
                                                                            gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTAATATCTTGAACAGGCGCGCCCGAAATCTAATTCCCCCATGATCATCATCGATCCGCGCTATACCCACAcc                                                                         >  4:1716560/1‑150 (MQ=255)
                                                                                   aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCgg                                                                 <  1:829043/151‑1 (MQ=255)
                                                                                    aaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGGTTACTATCTTGAAAAAGCACGCCAGAAAACTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCGACACCGGGGGCggg                                                                >  1:1129052/1‑151 (MQ=255)
                                                                                     aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgc                                                               <  4:1562345/151‑1 (MQ=255)
                                                                                     aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTTTTTTGAAAAGGCACGCCAGAAAACTAATGCCCCCCTGATAATCATCTAACCGCGCTTTACCGACACccgggccggg                                                                >  2:331789/1‑150 (MQ=255)
                                                                                         gTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGa                                                           <  3:12194/151‑1 (MQ=255)
                                                                                           aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGCACGCCCGAAAACTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCCACACCGGGGCCGGGggcg                                                            >  4:858815/1‑148 (MQ=255)
                                                                                            aGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTATTATCTTTAAAAGGCACGCCCGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCCACACCCGGGCCGGgcg                                                              >  1:165153/1‑145 (MQ=255)
                                                                                            aGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCGTTCTATCTTTAAAAAACACGCCAGAAATATAAAGGCCGCCTGTTAATAATTGATCCGCGCTTTACCGACAAccgggccggggggg                                                            >  3:184872/1‑143 (MQ=255)
                                                                                             gCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGAt                                                       <  1:1890616/151‑1 (MQ=255)
                                                                                                ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGa                                                     <  4:1583567/150‑1 (MQ=255)
                                                                                                ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGa                                                     <  1:1918312/150‑1 (MQ=255)
                                                                                                 gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGTCTCGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGt                                                   <  3:498195/151‑1 (MQ=255)
                                                                                                   aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGCACGCCCGAAATCTAATACCCGCATGATCATCATCTATCCCCGCTATACCGCCACCGGGGCCGGGCGCGGAGATGAGTgg                                                 >  1:1895742/1‑151 (MQ=255)
                                                                                                    gTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGa                                                <  3:87013/151‑1 (MQ=255)
                                                                                                     tGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGT‑GCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATc                                              >  1:1723329/1‑151 (MQ=255)
                                                                                                      gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATc                                              <  3:234008/151‑1 (MQ=255)
                                                                                                       cTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATc                                              <  1:1542753/150‑1 (MQ=255)
                                                                                                          tttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTACTATCTTGAACAGGCACGCCAGAAAACTAATGCCCGCATGATCATCAACGATCCCCGCTATACCCCCACCCGTGCCGGGGGCGAAGATGAGTGGGTccc                                            >  2:547166/1‑149 (MQ=255)
                                                                                                          tttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCt                                           <  1:527446/150‑1 (MQ=255)
                                                                                                           ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGAGTTACTATCTTTAAAAAGCACGCCCGAAAACTAAAGCCCGCACGATCATCCTCGATCCGCGCTATACCGCCACCCGGGCCGGGCGCGCAGATGAGGGGAGcccca                                          >  1:602405/1‑148 (MQ=255)
                                                                                                           ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACCATCTTTAAAAGGCACGACAGAAAACTAATGCCCGCATGATCATCATCGGTCCCCGCTAAACCGAACACGGGGCCGGGCGCGGAGATGAGTGGATCCCTAt                                         >  4:1792324/1‑151 (MQ=255)
                                                                                                            tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCCGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATAACGACACCGGGGCCGGGGGCGAAGATGAGTGGGTCCCTAtt                                        >  2:496705/1‑151 (MQ=255)
                                                                                                                aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGt                                     <  1:566362/150‑1 (MQ=255)
                                                                                                                aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCAGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTc                                    <  1:972629/151‑1 (MQ=255)
                                                                                                                   aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCg                                  <  1:1642765/150‑1 (MQ=255)
                                                                                                                   aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCg                                  <  4:1313098/150‑1 (MQ=255)
                                                                                                                    aCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCggg                                <  4:943152/151‑1 (MQ=255)
                                                                                                                       cTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTAATATCTTTAAAAAGCCCGCCCGAAAACTAAAGCCCCCATGAACATCAACGATCCGCGCTATACCGACACCCGGGCCGGGCGCGGAGAAGAGTGGATCCCCATTTGTCCgggga                              >  4:1733622/1‑148 (MQ=255)
                                                                                                                           cGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGAt                         <  4:1265055/151‑1 (MQ=255)
                                                                                                                             aaaCGCGAATGAGTGGCGGTGGGGGGGCTTAATATCTTTAACAGGCACGCCAGAAATCTAATGCCCCCATGATCATCATCGATCCCCGCCATAACGACACCGGGGCCGGGCGCGAAGATGAGGGGATCCCTATTCGTCCGGGGAAAGa                          >  1:1177910/1‑148 (MQ=255)
                                                                                                                             aaaCGCGAATGAGTGGCGGTGGGGGGACTTACTTTTTTTAACAAGCACGCCAGAAAATCAATGCCCGCCTGATTATCATCGATTCGCGCTATAACCACACCCGGGCCGGGCGCGGAGAAGAGCGGATCCCTATTCGTCCggg                                >  1:449357/1‑142 (MQ=255)
                                                                                                                               aCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCg                     <  3:401990/151‑1 (MQ=255)
                                                                                                                                cgcgAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGc                    >  2:121471/1‑151 (MQ=255)
                                                                                                                                 gcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCa                   <  4:941083/151‑1 (MQ=255)
                                                                                                                                  cgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCAc                  <  1:886106/151‑1 (MQ=255)
                                                                                                                                     aTGAGTGGCGGTGGGGGGACTTTATATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCCACACCGGTGCCGGGCGCGGAGATGAGTGGGTCCCCATTCGTCCGGGAACAGATGCCGCCCTgg               >  3:450999/1‑151 (MQ=255)
                                                                                                                                       gAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGtt             <  1:1575591/151‑1 (MQ=255)
                                                                                                                                       gAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTgg               >  2:372557/1‑149 (MQ=255)
                                                                                                                                         gTGGCGGTGGGGTGACTTACTATCTTGAGCAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTaa           <  3:49004/151‑1 (MQ=255)
                                                                                                                                         gTGGCGGTGGGGGGACTTACTATCTTTAAAAGGCACGCCCGAAATCTAAAGCCCGCATGATCATTATTGATCCGCGCTATACCGACACCCGGGCCGGGGGCGAAGATGAGTGGATCCCCCTTTGTCCGGGGACAGATGCCGCACTgggt             >  3:288704/1‑147 (MQ=255)
                                                                                                                                           ggCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACg         <  4:615122/151‑1 (MQ=255)
                                                                                                                                           ggCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACg         <  1:581025/151‑1 (MQ=255)
                                                                                                                                           ggCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACg         <  1:143926/151‑1 (MQ=255)
                                                                                                                                             cGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGt       <  1:38901/151‑1 (MQ=255)
                                                                                                                                             cGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGt       >  1:31851/1‑151 (MQ=255)
                                                                                                                                                  gggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGc  <  2:298001/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGC  >  NC_000913/941564‑941860

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