Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913802,1180AC46.8% 105.5 / 15.0 77T78P (ACC→CCC) ybhIputative DASS family tricarboxylate or dicarboxylate transporter
Reads supporting (aligned to +/- strand):  ref base A (29/12);  new base C (0/36);  total (29/48)
Fisher's exact test for biased strand distribution p-value = 6.92e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGGTC  >  NC_000913/801976‑802265
                                                                                                                                                 |                                                                                                                                                   
aCTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGTGCGTTTAAAACCAccgcc                                                                                                                                             <  2:1023157/151‑1 (MQ=255)
  tGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGt                                                                                                                                           >  2:348783/1‑151 (MQ=255)
     cGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTAtt                                                                                                                                        >  1:498658/1‑151 (MQ=255)
      gTGGGTGTTTTTTGTTTTTTCCCTGGGGGCCATTGTGGGGCTGGTTATCAAGCTTTT‑CCC‑GGACCCTGTCTTACTTTTAATTGCCGTTGCTGCCTAAATGGTGGTGGTCGGTAAC‑TTACCCGACGGTGCTTTTAAACCCACCGCCGTAtt                                                                                                                                        <  1:1005678/150‑1 (MQ=255)
       tGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTATTa                                                                                                                                       >  1:111629/1‑150 (MQ=255)
        gggTGTTTTTTGGTATTTCCCTGGCGGCCTTTGTGGGGCTGTTTATCAACCTTTT‑CCC‑GGACCCTTTCGTACTGTTAATTGCCGTTGCTGCCTCAAGGGTGGTGGTCGGTAAT‑TTACCCGACGGTGCTTTTAAACCCACCGCCGTATTaa                                                                                                                                      <  1:634898/150‑1 (MQ=255)
              tttttGGTATTTCCCTGGGGGCCTTTGGGGGGCTGGTTACCAAGCCTTT‑CCC‑GGACCCTGTCGTATTGTTAATTGCCGTTGCTCCCTCAAGGGTGGTGGTCGGTACT‑TTACCCGACGGTGCTTTAAAAACCACCGCCGTATTAAGCgg                                                                                                                                  <  1:712758/148‑1 (MQ=255)
                tttGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTAc                                                                                                                              >  1:437734/1‑150 (MQ=255)
                   ggTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCtt                                                                                                                          >  1:700543/1‑151 (MQ=255)
                     tATTTCCCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTC‑CCC‑GGACCCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTACCCGAGGGTGCTTTTAAAACCACCGCCGTATTAAGCGGTTActct                                                                                                                           <  2:347718/148‑1 (MQ=255)
                       tttCCCGGGCGCCCATTGTGGGGCTGGTTATCAA‑CCTTT‑CCCGGGACCCTGTCGTCCTGTTAATTGCCGTTGCTGCCTCAAGGGTGGTGGTCGGTAAC‑TTACCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCtt                                                                                                                          <  2:135290/147‑1 (MQ=255)
                       tttACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAgg                                                                                                                      >  2:20026/1‑151 (MQ=255)
                       tttACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAgg                                                                                                                      >  2:550163/1‑151 (MQ=255)
                              ggcggcCATTGGGGGGCTGGTTATCAAGCCTTC‑CCC‑GGACCCTGTCGTTCTTTTAATTCCCGTTGCTGCCTCAAGGGTGGTGGTCGGTATC‑TTACCCGACGGTGCGTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTaccac                                                                                                                <  2:237463/150‑1 (MQ=255)
                              ggcgccCATTGTGGGGCTGTTTACAAAGCCTTC‑CCC‑GGACCCTGTCTTACTTTTAATTGCCGTTGCTCCCTCAATGGGGGTGGCGGGTAAT‑TTACCCGGCGGTGCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTacca                                                                                                                 <  1:535165/149‑1 (MQ=255)
                               gcggcattgggggggCGGTTTATCAAGCTTTC‑CCC‑GGACCCTGTGGTCCTTTTAATTGCCTTTGCTGCCTAAAGGGGGGTGGTGGGTACC‑TTACCCGACGTTGCGTTTAAACCCCCCGCCGTATTAAGCGGTTACTCTTCAGGTacca                                                                                                                 <  2:1011723/137‑1 (MQ=255)
                                  gcCATTGGGGGGCTGGTTATCAAGCCTTT‑CCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTACCCGACGGTGCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctgg                                                                                                            <  1:542777/150‑1 (MQ=255)
                                                ggTTATCAAGCCTTC‑CCC‑GGAACCTGTCGTACTTTTAATTGCGGTTGCTGCCTCAATGGGGGTGGTCGGTATT‑TTATCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGc                                                                                              <  1:867202/150‑1 (MQ=255)
                                                 gTTATCAAGCCTTT‑CCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGt                                                                                            >  2:430889/1‑151 (MQ=255)
                                                  ttATCAACCTTTT‑CCC‑GGACCCTGTCTTACTTTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCg                                                                                             <  1:940158/149‑1 (MQ=255)
                                                     tCAAGCCTTT‑CCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTAc                                                                                        >  1:759536/1‑151 (MQ=255)
                                                       aacccTTT‑CCG‑GGACCCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAAGGGTGGTGGTCGTTAAC‑TTATCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTAc                                                                                        <  2:606060/146‑1 (MQ=255)
                                                        aaccTTTCCCC‑GGAACCTGTTTTACTTTTATTTGCCTTTGCTGCCTCAATGGTGGTGGTCGGTAAT‑TTATCCGACGGTGCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCt                                                                                      <  2:440385/149‑1 (MQ=255)
                                                         gCCTTT‑CCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCtt                                                                                     >  1:512539/1‑150 (MQ=255)
                                                           cTTT‑CCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTa                                                                                    >  2:708732/1‑149 (MQ=255)
                                                            ttt‑CCG‑GGACCCTGTCTTCTTGTTATTTGCCGTTGCTGCCTCAAGGGTGGTGGTCGTTAAC‑TTATCCGACGTTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTaa                                                                                   <  2:1161791/149‑1 (MQ=255)
                                                             tt‑CCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgcg                                                                                >  2:185845/1‑151 (MQ=255)
                                                             tt‑CCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgcg                                                                                >  1:759937/1‑151 (MQ=255)
                                                             tt‑CCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgcg                                                                                >  2:943575/1‑151 (MQ=255)
                                                                ccc‑GGCTCCTGTCGTTTTGTTAATGCCGGTTGCTTCCTCAATGG‑GGGGGATAGTAACTTTACCCGGTTGGGCGTTTAAACCCGCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgct                                                                                <  1:1056129/149‑2 (MQ=255)
                                                                ccc‑GGACCCTGTCGTCCTTTTAATTGCCGTTGCTCCCTCAATGGTGGTGGTCGGTAAT‑TTACCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgcgc                                                                               <  1:558967/150‑1 (MQ=255)
                                                                ccc‑GGACCCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTATT‑TTACCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgc                                                                                 <  1:1082680/148‑1 (MQ=255)
                                                                ccc‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAAGGGGGGTGGTCGGTAAT‑TTACCCGAGGTTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgcgc                                                                               <  2:579547/150‑1 (MQ=255)
                                                                      acccTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAT‑TTATCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAt                                                                          <  1:335666/148‑1 (MQ=255)
                                                                      aaCCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAtt                                                                         >  1:246054/1‑151 (MQ=255)
                                                                      aaCCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAtt                                                                         >  2:738897/1‑151 (MQ=255)
                                                                           gTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTACCCGACGGGGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtgt                                                                      <  2:1199890/149‑1 (MQ=255)
                                                                             cGTTCTTTTAATGGCCGTTGCTGCCTCACTGGGGGTGTTCGTTATT‑TTATCGGCCGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtgt                                                                      <  1:703172/147‑1 (MQ=255)
                                                                                    ttATTTCCCTTTGCTCCCCAAATGGTGTTGGTCGTTATC‑TTACCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGaccac                                                                <  1:358781/146‑1 (MQ=255)
                                                                                    ttAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGacca                                                                 <  1:965518/145‑1 (MQ=255)
                                                                                    ttAATTCCCGTGGCTGCCTCAAGGGTGGTGGTCGGTAAC‑TTTTCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCgg                                                             <  2:543468/149‑1 (MQ=255)
                                                                                    ttAATGGCCGTTGCTGCCTAAAGGGTGGTGGTCGGTAAC‑TTACCCGACGGTGCTTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCgg                                                             <  2:104842/149‑1 (MQ=255)
                                                                                     tAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGttt                                                          >  2:36085/1‑151 (MQ=255)
                                                                                              tttccgcCTAAAGGGTGGTGGTCGTAAAC‑TTATCCGACGTTGTTTTAAACACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTaaa                                                   <  1:995596/144‑1 (MQ=255)
                                                                                              ttgcgcccTAAAGGGTGGTGGTCGGTAAC‑TTACCCGAGGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAg                                                        <  1:1182832/138‑1 (MQ=255)
                                                                                                gctgcCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGt                                                <  2:1164327/150‑1 (MQ=255)
                                                                                                   gcCTCAATGGTGGGGGTGGGTAAC‑TTATCCGGCGGTGCTTTTAAACCCCCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAt                                              <  1:916011/149‑1 (MQ=255)
                                                                                                     cTCATTGGGGGTGGTCGTTACT‑TTACCCGAGGGTGCTTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAtt                                             <  1:260905/148‑1 (MQ=255)
                                                                                                        agtggtggtggTTGGAAAT‑TAATCCGTCGGGGCTTTAAAAACCTCCGCCGTATTAAGCGGTTACGCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTAt                                       <  1:730225/149‑1 (MQ=255)
                                                                                                        aaTGGTGGTGGTCGGTAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTAt                                       >  2:785828/1‑151 (MQ=255)
                                                                                                        aaTGGGGGTGGTGGGTCAT‑TTACCCGGCGGTGTTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTg                                            <  1:833494/146‑1 (MQ=255)
                                                                                                        aaTGGGGGTGGTCGGTAAC‑TTACCCGACGGTGCGTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTAt                                       <  1:550137/151‑1 (MQ=255)
                                                                                                           ggtggtggtTGGTACC‑TTATCCGGCGGTGCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTa                                        <  1:711773/147‑1 (MQ=255)
                                                                                                           ggtggtggtGGGTAAC‑TTACCCGACGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCt                                         <  2:613376/146‑1 (MQ=255)
                                                                                                           ggggggggtCGTTACC‑TTATCCGACGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATct                                     <  2:1150321/144‑1 (MQ=255)
                                                                                                            gtggtggtCGGTAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgc                                   >  1:1947/1‑151 (MQ=255)
                                                                                                              ggtggtGGGTACC‑TTATCCGACGGTCCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgc                                   <  2:842658/149‑1 (MQ=255)
                                                                                                              ggtggtCGGTAAT‑TTACCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgc                                   <  1:904237/149‑1 (MQ=255)
                                                                                                              ggtggtCGGTAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgctg                                 >  1:748229/1‑151 (MQ=255)
                                                                                                                  ggcGGTAAC‑TTATCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGAt                               <  1:148952/147‑1 (MQ=255)
                                                                                                                   tCGGTAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTgg                            >  2:194728/1‑151 (MQ=255)
                                                                                                                    cGGTAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGt                           >  1:387733/1‑151 (MQ=255)
                                                                                                                     ggTAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGt                           >  2:127116/1‑150 (MQ=255)
                                                                                                                       tAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTaa                         >  2:163227/1‑150 (MQ=255)
                                                                                                                         acctttcccGACGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTaaa                        <  1:147218/143‑1 (MQ=255)
                                                                                                                         aC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAAt                      >  1:24915/1‑151 (MQ=255)
                                                                                                                              acccGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGt                  <  1:908525/149‑1 (MQ=255)
                                                                                                                               tCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGt                  >  2:546166/1‑150 (MQ=255)
                                                                                                                                ccgccgGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGt                  <  2:608744/149‑1 (MQ=255)
                                                                                                                                ccGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTa                 <  1:1060150/150‑1 (MQ=255)
                                                                                                                                ccGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTACGCGCCGCATTTGTGACCACCGGTTTAGGTAAAGGTATTGCCTATCTGCTGATTGGGAAAATCGGTa                 >  2:629790/1‑150 (MQ=255)
                                                                                                                                     ggTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAacac             <  2:111846/149‑1 (MQ=255)
                                                                                                                                     ggTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAAcacc            <  2:191937/150‑1 (MQ=255)
                                                                                                                                       tGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACg         >  1:981098/1‑151 (MQ=255)
                                                                                                                                        gCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAAcaccac          <  1:696834/149‑1 (MQ=255)
                                                                                                                                           tttAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACg         <  2:512563/147‑1 (MQ=255)
                                                                                                                                           tttAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgg     >  2:549470/1‑151 (MQ=255)
                                                                                                                                           ttaaaaaCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGct       <  2:673995/146‑1 (MQ=255)
                                                                                                                                              aaaaCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtc  >  1:703821/1‑151 (MQ=255)
                                                                                                                                              aaaaCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtc  >  2:243826/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGGTC  >  NC_000913/801976‑802265

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