Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913941,7130TG40.5% 131.3 / 26.3 74V252G (GTG→GGG) dmsAdimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base T (4/40);  new base G (30/0);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 9.89e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941563‑941862
                                                                                                                                                      |                                                                                                                                                      
gCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGt                                                                                                                                                        <  1:428922/151‑1 (MQ=255)
  gACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGa                                                                                                                                                      <  2:799891/151‑1 (MQ=255)
        tcctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCtt                                                                                                                                                   >  1:747585/1‑148 (MQ=255)
        tcctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCtt                                                                                                                                                   >  2:541519/1‑148 (MQ=255)
         cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTAc                                                                                                                                                 <  2:336712/149‑1 (MQ=255)
          ctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTAt                                                                                                                                              >  1:395232/1‑151 (MQ=255)
                 cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTTGGGGGGGTTTCTTTCTTgga                                                                                                                                        >  2:1091288/1‑148 (MQ=255)
                 cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAc                                                                                                                                       <  1:217066/151‑1 (MQ=255)
                 cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAc                                                                                                                                       <  1:1253017/151‑1 (MQ=255)
                  aaaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACa                                                                                                                                      <  2:249475/151‑1 (MQ=255)
                  aaaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACa                                                                                                                                      <  1:1206752/151‑1 (MQ=255)
                  aaaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTaa                                                                                                                                        >  1:409236/1‑149 (MQ=255)
                    aTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAgg                                                                                                                                    <  2:1057587/151‑1 (MQ=255)
                       gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCa                                                                                                                                  <  1:311135/150‑1 (MQ=255)
                        cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAAGCAc                                                                                                                                 >  1:1191938/1‑150 (MQ=255)
                           aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGcc                                                                                                                              <  1:595326/150‑1 (MQ=255)
                            aGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAg                                                                                                                            <  2:935960/151‑1 (MQ=255)
                                 tGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAAt                                                                                                                       <  1:1234557/150‑1 (MQ=255)
                                      tataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGGACGCCCG‑AAAACTaa                                                                                                                   >  2:1068417/1‑149 (MQ=255)
                                       ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGGCGTATTATATTTGACAAGCACGccccaaaacttaag                                                                                                                  >  1:222835/1‑138 (MQ=255)
                                        taCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATTTTTAAAAAGCACGCCCG‑aaaaatta                                                                                                                   >  2:1000249/1‑142 (MQ=255)
                                        taCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATCTTTAAAAAGCACGCCCA‑AAAATTAATGcc                                                                                                               >  2:1005816/1‑151 (MQ=255)
                                         aCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGccc                                                                                                              <  1:541518/151‑1 (MQ=255)
                                         aCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTCCTATTTTTGAAAGGCACGcccccaaactgaaagc                                                                                                                >  1:241148/1‑136 (MQ=255)
                                            tACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTTTCTTTTAAAAGGACGCCCG‑AAAAATAAAGccccc                                                                                                            >  2:726941/1‑148 (MQ=255)
                                             aCGGCGGCTGGGCAGATGGCAACAGCCCGTCCGATATCGAAAACAGTAAGCTGGGAGTGCTGGTTGGAAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGATTATTATCTTTAACAAGAACCCCAC‑aaaaccaaataacgc                                                                                                            >  1:180561/1‑137 (MQ=255)
                                                 cggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTACTATCTTTAAAAAGCACGCCCGAAAATCTaaaccccc                                                                                                             >  1:902077/1‑139 (MQ=255)
                                                     tGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTGAACAAGCACGCCCG‑AAATCTAATGCCCGCATGa                                                                                                        >  1:1030294/1‑145 (MQ=255)
                                                      gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGatcatc                                                                                                   <  2:143016/149‑1 (MQ=255)
                                                       ggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGatcatcatc                                                                                                <  1:1237542/151‑1 (MQ=255)
                                                         cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGa                                                                                              <  2:40094/151‑1 (MQ=255)
                                                                      cccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAcc                                                                                 <  1:835573/151‑1 (MQ=255)
                                                                      cccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAcc                                                                                 <  1:1069247/151‑1 (MQ=255)
                                                                      cccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAcc                                                                                 <  2:534826/151‑1 (MQ=255)
                                                                       ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCAG‑AAAACTAATGCCCCCATGATCATCATCGATCCCCCCTATAcc                                                                                 >  2:1162005/1‑150 (MQ=255)
                                                                         gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGac                                                                              <  1:381573/151‑1 (MQ=255)
                                                                         gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGac                                                                              <  2:141494/151‑1 (MQ=255)
                                                                         gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGac                                                                              <  2:409237/151‑1 (MQ=255)
                                                                           cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATCTTTAACAGGCACGCCCG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGacac                                                                            >  1:808548/1‑151 (MQ=255)
                                                                             gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCg                                                                          <  2:348720/151‑1 (MQ=255)
                                                                              atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAAGCACGCCCG‑AAAACTAATGCCCGCATGATCATCATCGATCCCCGCTATACCGACAcccg                                                                         >  1:927395/1‑149 (MQ=255)
                                                                              atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTAATATTTTGAAAAAGCACGCCAG‑AAAACTAATGCCCGCATGATCAACATCGCACCGCGCCATACCGAAAcccg                                                                         >  1:1128155/1‑149 (MQ=255)
                                                                                atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGTGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTg                                                                       <  1:10131/151‑1 (MQ=255)
                                                                                atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGt                                                                        <  2:1192630/150‑1 (MQ=255)
                                                                                     aaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGTGGCTTACTATCTTTAACAGGCACGCCCG‑AAATATAATGCCCCCATGCTCATCATCGTTCCGCGCTATACCGACACCCGTGCCggg                                                                  >  1:739731/1‑151 (MQ=255)
                                                                                            aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAg                                                           <  1:1165990/151‑1 (MQ=255)
                                                                                             aGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCCCGCCCG‑AAATCTAATGCCCGCATGATTATTATCGGTCCGCGCCATACCCACACCGGTGGCggg                                                                  >  1:558560/1‑143 (MQ=255)
                                                                                             aGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGGCCG‑AAATATAATTCCCCCGTGATCATCATCGGTCCCCGCTATACCCACACccgggccggggg                                                                >  2:808640/1‑143 (MQ=255)
                                                                                               cTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTTCTATCTTGAAAACGCAAGCCCG‑AAAACTAATGCCCGCATGAACCTTATCGATCCGCGCTATACCGACACCGGGGCCCGgcg                                                                >  2:994711/1‑143 (MQ=255)
                                                                                                 ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATg                                                        <  1:726939/149‑1 (MQ=255)
                                                                                                    aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCCGGGCCGGGCGCGAAGATGAGTgg                                                   >  1:1229675/1‑151 (MQ=255)
                                                                                                     gTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGa                                                  <  1:355627/151‑1 (MQ=255)
                                                                                                       gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATc                                                <  2:596684/151‑1 (MQ=255)
                                                                                                            ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAAAAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGGAGATGAGGGGATccc                                              >  1:932193/1‑148 (MQ=255)
                                                                                                             tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAt                                           <  2:31128/150‑1 (MQ=255)
                                                                                                              ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAAGCACGCCCG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAtt                                          >  2:1210512/1‑150 (MQ=255)
                                                                                                                taataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGGCCGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAACGACACCGGGGCCGGGGGCGAAGATGAGGGGGTCCCTATTCGt                                       >  2:939554/1‑151 (MQ=255)
                                                                                                                 aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCTCGAAGATGAGTGGATCCCTATTCGTc                                      <  1:372152/151‑1 (MQ=255)
                                                                                                                 aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAAGCACGCCCG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCCACACCCGGGGCGGGGGCGAAAATGGGGGGGACCCCATTCGt                                       >  1:962421/1‑150 (MQ=255)
                                                                                                                   taaCCCTGGCGAAACGCGAATGAGTGGCGGTTGGGGGGCCTTCTATCTTTTAAAGGCCCGCCGG‑AAATCTAATGCCCGCGTGATTAACAACGAACCCCCCCAAAACCACACCGGTGCCGGGCGCGAAGGTGAGGGGGTTCCCAttt                                         >  2:274368/1‑145 (MQ=255)
                                                                                                                    aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCgg                                   <  1:728084/151‑1 (MQ=255)
                                                                                                                     aCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCggg                                  <  2:923465/151‑1 (MQ=255)
                                                                                                                      cccTGGCGAAACGCGAATGAGTGGCGGGGGGGTGGCCTTATTTTTTGTAAAGGGACGCCAG‑GAAACTAAAGGGCCCATGATTTTCATCGGTCCGCGCTATAACGGCACCGGGGGCGGGGGCGGAGGTGGGTGGGACCCTCTTCGTCCggg                                  >  1:102748/1‑150 (MQ=255)
                                                                                                                        cTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACCTACTTTCTTTAACAGGGACGCCAG‑AAATCTAATGCCCGCCTGATCATTATCGATCCGCGCTATACCCACACCCGTGCCGGGCGCGGAGATGAGTGGATCCCTATTCCTCCgggga                                >  2:75968/1‑148 (MQ=255)
                                                                                                                             gAAACGCGAATGAGTGGCGGTTGGGGGGCTTACCTTATTTAACACGGAAGCCCC‑AAAACCAATTCCCGGATTATCATCAACGAGTCTCCCGATACCCACACCGGGGCCGGGGGCGGAGATGCGTGGATCCCCATTCCTCCgggca                                >  2:173758/1‑143 (MQ=255)
                                                                                                                             gAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGAACCCTATTCGTCCGGGGACAGAt                           >  2:477583/1‑150 (MQ=255)
                                                                                                                              aaaCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATg                          <  1:818066/150‑1 (MQ=255)
                                                                                                                               aaCGCGAATGAGGGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGTGAACAGATGc                         <  1:765507/150‑1 (MQ=255)
                                                                                                                                   cgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCAc                    <  1:409748/151‑1 (MQ=255)
                                                                                                                                     aaTGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTg                  <  1:866160/151‑1 (MQ=255)
                                                                                                                                         aGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTa              <  1:558072/151‑1 (MQ=255)
                                                                                                                                          gTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTaa             <  2:1072037/151‑1 (MQ=255)
                                                                                                                                            ggCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAAc            <  2:941187/150‑1 (MQ=255)
                                                                                                                                            ggCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACg           <  1:1257580/151‑1 (MQ=255)
                                                                                                                                               ggTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTc        >  1:1035522/1‑151 (MQ=255)
                                                                                                                                                      tGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGt  >  1:1239916/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
GCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941563‑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: ‑