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

TAAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGGTCTGG  >  NC_000913/801969‑802268
                                                                                                                                                      |                                                                                                                                                       
tAAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTaaaa                                                                                                                                                         <  2:1002558/150‑1 (MQ=255)
  aGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAa‑cca                                                                                                                                                     <  1:835534/151‑1 (MQ=255)
   gCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAa‑ccac                                                                                                                                                    >  4:441547/1‑151 (MQ=255)
    cGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAa‑ccacc                                                                                                                                                   <  4:95179/151‑1 (MQ=255)
     gAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCAccg                                                                                                                                                  <  2:221576/151‑1 (MQ=255)
       aCTGGCGTGGGTGCTTTTTGGGATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑ACAccgcc                                                                                                                                                >  2:302133/1‑151 (MQ=255)
            cGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTAtt                                                                                                                                           >  2:578174/1‑151 (MQ=255)
             gTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTAtt                                                                                                                                           >  3:363602/1‑150 (MQ=255)
               gggTGTTTTTTGGTATTTCCCTGGCGGCCATTGTGGGGCTGGTTATCAACCTTTTCCC‑GGAACCTGTCTTCTTTTTATTTCCCGTTGCTGCCTCAATGGTGGTGGTGGTTAACTTACCCGACGGTGCTTTAAAAC‑CCACCGCCGTATTa                                                                                                                                          <  2:961771/149‑1 (MQ=255)
                   gCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCgg                                                                                                                                     >  2:583249/1‑150 (MQ=255)
                     tttttGGTATTTCCCTGGCGGCCATTGTGGGGCTGTTTATCAACCTTTCCCC‑GGACCCTGTCGTATTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTACCCGACGGTGCTTTAAAAC‑CCACCGCCGTATTAAGCgg                                                                                                                                     <  2:984636/148‑1 (MQ=255)
                     ttttgGGTATTTCCCTGGCGGCCTTTGTGGGGCTGTTTATAAAGCCTTCCCC‑GGACCCTGTCGTACTGTAAATTGCCGTTGCTGCCTCAAGGGTGGTGGTCGTTAATTTACCCGACGGTGCGTTAAAAC‑CCACCGCCGTATTAAGCGGtt                                                                                                                                   <  2:585294/150‑1 (MQ=255)
                      ttttGGTTTTTACCGGGGGGCTTTTGTGGGGCTGTTTATCAAGCCTTTCCC‑GGACCCTGTCGTACTGTTATTTGCCTTTGCTGCCTCAATGGTGGTGGTGGGTAACTTATCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTa                                                                                                                                  <  3:260416/150‑1 (MQ=255)
                         tGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTActct                                                                                                                              >  1:891570/1‑151 (MQ=255)
                              tttgcgggggggcCGTTGAGGGGCTGTTTATAAGGCATTCCCG‑GGATCCTGTCTTACTTTTTATTCCCGTTGCTGCCCCAGGGGTGGTGGTCGTTAACTTACCCGACGGTGCGTTTACAA‑CCGCCGCCGTATTAAGCGGTTACTCTTc                                                                                                                            <  1:407472/138‑1 (MQ=255)
                              tttccgtGGCGGCTATTGTGGGGCTGGTTATCAAGCCTTCCCG‑GGACCCTTTCGTACTTTTAATTGCCTTTGCTGCCTCAATGGTGGTGTTCGTTACCTTATCCGCCGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCa                                                                                                                           <  3:63661/143‑1 (MQ=255)
                              tttCCCTGGCGGCCATTGGGGGGCTGGTTATCAAGCCTTTCCC‑GGACCCTGTCTTACTTTTATTTGCCGTTGCTGCCTCAATGGTGGTGGTCGTTAATTTACCCGACGGTGCGTTAAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCa                                                                                                                           <  2:523536/149‑1 (MQ=255)
                               ttcccgggcggcCTTGGTGGGGTTGGTTATCACGCCTTTCCC‑GGACCCTGTCGTCCTGTTAATTGCCTTTGCTGCCTCAAGGGGGGTGGTCGTTACTTTACCCGAGGGTG‑GTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAgg                                                                                                                         <  3:21338/144‑1 (MQ=255)
                               ttACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGt                                                                                                                        >  3:282318/1‑151 (MQ=255)
                                    tGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTaccac                                                                                                                   >  2:1199762/1‑151 (MQ=255)
                                    tGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTaccac                                                                                                                   >  1:845485/1‑151 (MQ=255)
                                    tGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTaccac                                                                                                                   >  2:64407/1‑151 (MQ=255)
                                    tGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTaccac                                                                                                                   >  1:591571/1‑151 (MQ=255)
                                        ggcCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctgg                                                                                                               >  2:1047564/1‑151 (MQ=255)
                                          ccATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctggct                                                                                                             >  2:417142/1‑151 (MQ=255)
                                           cATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctggctg                                                                                                            >  3:71872/1‑151 (MQ=255)
                                             ttGTGGGGCTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGgt                                                                                                          >  3:161553/1‑151 (MQ=255)
                                              tttgGGGCTGTTTATCATCCCTTTCCC‑GGACCCTTTCGTCTTGTTATTTGCGGTGCCTCCCTCATTGGTGTTGGTCGGTATTTTCTCCGCCGTTCCGTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGgt                                                                                                          <  1:119800/148‑1 (MQ=255)
                                                tgGGGCTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGtt                                                                                                       >  1:880358/1‑151 (MQ=255)
                                                 ggggCTGTTTATAAAGCCTTTCCC‑GGAACCTGTCGTACTTTTATTTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTATTTTACCCGAGGGTGCGTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGtt                                                                                                       <  2:722381/150‑1 (MQ=255)
                                                  gggTTGTTTATCAAGCCTTTCCC‑GGACCCTGTCGTACTGTTAATTGCGGTTGCTGCCTAAATGGTGGTGGTCGTTACTTTATCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGtttc                                                                                                     <  1:877937/151‑3 (MQ=255)
                                                  gggTTGTTTATCAA‑CCTTTCCCGGGACCCTGTCTTACTTTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTGGTTAATTTACCCGGCGGTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTtct                                                                                                     <  2:175803/151‑1 (MQ=255)
                                                  gggCTGTTTATCAACCTTTTCCC‑GGACCCTGTCTTACTGTTAATTGCCGTTGCGGCCTAAATGGTGGTGGGCGGTACTTTGCCCGAGGGGGCTTTAAAAA‑CCTCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTtct                                                                                                     <  1:906694/151‑1 (MQ=255)
                                                   ggCTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAATTTATCCGACGGTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTtctc                                                                                                    <  4:260516/151‑1 (MQ=255)
                                                   ggCTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTtctc                                                                                                    >  2:1165139/1‑151 (MQ=255)
                                                     cTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCgg                                                                                                  >  1:1251014/1‑151 (MQ=255)
                                                     cTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCgg                                                                                                  >  2:944793/1‑151 (MQ=255)
                                                      tGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGc                                                                                                 >  2:1217702/1‑151 (MQ=255)
                                                         ttATCAAGCCTTCCCC‑GGACCCTGTCGTTCTGTTAATTGCGGTTGCTGCCTCAATGGTGGTGGTCGGTATCTTATCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCg                                                                                                   <  2:503883/146‑1 (MQ=255)
                                                           aTCAAGCTTTTCCC‑GGACCCTGTCTTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGGGGTCGGTATTTTACCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGttt                                                                                             <  2:245118/150‑1 (MQ=255)
                                                           aTCAAGCCTTCCCC‑GGACCCTGTCGTAGTGTTAATGGCCTTTGCTGCCTAAAGGTGGGGGGTGGTTACCTTACCCGACGGTGCGTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTa                                                                                            <  2:1269357/151‑1 (MQ=255)
                                                                 ccTTTCCC‑GGACCCTGTCGTACTTTTAATTCCCTTTGCTGCCTCAATGGTGGTGGTCGGTAATTTACCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTa                                                                                       <  1:167348/150‑1 (MQ=255)
                                                                 ccTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTa                                                                                       >  4:157532/1‑150 (MQ=255)
                                                                   ttcccc‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTAAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTa                                                                                       <  2:391147/145‑1 (MQ=255)
                                                                    ttCCC‑GGAACCTGTCGTACTGTTAATTCCCGTTTCTCCCAAATTGGGGGTGTTCGATAACTTATCCGGCGGTGCTTTAAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTa                                                                                       <  2:977429/147‑1 (MQ=255)
                                                                     tCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgcgc                                                                                  >  2:990042/1‑151 (MQ=255)
                                                                      ccc‑GGAACCTGTCGTACTTTTAATTGCCGTTGCTGCCTAAATGGTGGTGGTCGGTAATTTATCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgcg                                                                                   <  4:158060/149‑1 (MQ=255)
                                                                        ccgaacccTGTCTTACTGTTATTTGCCGTTGCTCCCTAAAGGGTGGTGGTCGTAATTTTACCCGCCGGTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGcgccg                                                                                <  1:1237292/145‑1 (MQ=255)
                                                                        c‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGcgccgc                                                                               >  4:267980/1‑151 (MQ=255)
                                                                            aaCCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAtt                                                                            >  2:1107035/1‑151 (MQ=255)
                                                                            aaCCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGGGTTCTCGGCGTTTACCTTAAGCGCCGCAt                                                                             >  2:828961/1‑150 (MQ=255)
                                                                              ctcgtggTATTGTTAATGGCCGTTGTTCCCTATGTGGTGGTGGTGGATATCTTTCCCGAGGGTGAGTTAAAAC‑CCCCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAttt                                                                           <  2:1058843/144‑1 (MQ=255)
                                                                              ccTGTCGTACTTTTAATGGCCGTTGCTGCCTCAATGGTGGTGGTGGGTAACTTATCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAttt                                                                           <  4:133383/150‑1 (MQ=255)
                                                                               cTGTCGTTCTGTTAATGCCCTTTGCTGCCTCAATGGTGGTGTTGGGTAACTTACCCGGGGGTGCTTTTAACA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAttt                                                                           <  1:405807/149‑1 (MQ=255)
                                                                                tGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtgtg                                                                        >  2:115748/1‑151 (MQ=255)
                                                                                 gTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtgtg                                                                        >  1:1071355/1‑150 (MQ=255)
                                                                                  tCGTTCTTTTAATTGCCGTTGCTGCCTAAATGGTGGTGGTCGGTAACTTACCCGACGGTGCGTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtgtg                                                                        <  1:1202073/149‑1 (MQ=255)
                                                                                  tCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGac                                                                      >  1:1117150/1‑151 (MQ=255)
                                                                                    gTACTGTTAGTTGCCTTTGCTCCCTCAAGGGTGGTGGTCGTTAATTTATCCGACGGTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGacca                                                                    <  2:1270625/151‑1 (MQ=255)
                                                                                     tACTGTTAATTGCCTTTCCTGCCTCATTGGTGGTGGTCGTTACCTTATCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGa                                                                       <  3:27213/147‑1 (MQ=255)
                                                                                        tGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCgg                                                                >  3:239695/1‑151 (MQ=255)
                                                                                          ttAATTGCCGTTGCTCCCTCAAGGGTGGTGGTCTGTAACTTACCCGACGGGGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCg                                                                 <  2:245672/148‑1 (MQ=255)
                                                                                            aaTTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGttt                                                             >  4:156110/1‑150 (MQ=255)
                                                                                            aaTGGCCGTTGCTGCCACAATGGTGGTGGTCGGTAACTTAGCCGTGGGTGCTTTTAAAC‑CCCCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGtt                                                              <  1:583363/149‑1 (MQ=255)
                                                                                              ttGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAgg                                                          >  2:355482/1‑151 (MQ=255)
                                                                                               tGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGt                                                         >  2:786769/1‑151 (MQ=255)
                                                                                                gCCGTTGCTGCCTCAGTGGTGGTGGTGGGTAACTTATCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGt                                                         <  3:442367/150‑1 (MQ=255)
                                                                                                 cctttGCTGCCTCAATGGTGGTGTTGGGTAATTTACCCGACGGTGCTTTAAAAC‑CCGCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAgg                                                          <  1:885712/145‑1 (MQ=255)
                                                                                                    ttgcgccctgatgggtggtggtCGTTACCTTATCCGACGGTGTTTTGAAAA‑CCTCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAAc                                                     <  1:302180/137‑1 (MQ=255)
                                                                                                         gcCTCAATGGGGGTGGTGGTTAACTTATCCGACGGTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAtt                                                <  4:161555/150‑1 (MQ=255)
                                                                                                          ccTCAATGGTGGTGGTCGGTAATTTATCCGACGGTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAtt                                                <  1:381117/149‑1 (MQ=255)
                                                                                                               aTGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATc                                         >  1:18888/1‑151 (MQ=255)
                                                                                                                 ggtggtggtCGTTAATTTACCCGACGGGGCTTTAAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCCCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctg                                       <  3:259382/151‑1 (MQ=255)
                                                                                                                 ggtggtggtCGGTAATTTACCCGAGGGTGCTTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGct                                             <  3:307636/145‑2 (MQ=255)
                                                                                                                 ggtggtggtCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATct                                        <  3:139506/150‑1 (MQ=255)
                                                                                                                    ggtggtGGGTACCTTATCCGACGGTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgct                                     <  1:528370/150‑1 (MQ=255)
                                                                                                                    ggtggtCGGTAACTTATCCGACGGTGTGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgct                                     <  4:426171/150‑1 (MQ=255)
                                                                                                                       ggtCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGAtt                                 >  1:583946/1‑151 (MQ=255)
                                                                                                                        gtCGTTAACTTCCCCGCGGTTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCGTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGAtt                                 <  2:554482/150‑1 (MQ=255)
                                                                                                                           ggTAATTTACCCGACGGTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGt                              <  2:1256829/150‑1 (MQ=255)
                                                                                                                              aaCTTACCGGGCGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTaaa                           <  1:930094/150‑1 (MQ=255)
                                                                                                                              aaCTTACCCGCCGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTgg                               <  1:448015/146‑1 (MQ=255)
                                                                                                                               aCTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAAt                         >  1:468913/1‑151 (MQ=255)
                                                                                                                                 ttctCCGACGGTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATc                        <  2:1226282/147‑1 (MQ=255)
                                                                                                                                 ttACCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTaaaa                          <  1:308863/148‑1 (MQ=255)
                                                                                                                                 ttACCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTaaaa                          <  1:614808/148‑1 (MQ=255)
                                                                                                                                    tCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGt                     >  3:281782/1‑150 (MQ=255)
                                                                                                                                        aCGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAacac                >  1:588686/1‑151 (MQ=255)
                                                                                                                                                ttaaaac‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctg         <  2:470263/147‑1 (MQ=255)
                                                                                                                                                 ttAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctggg       >  3:331003/1‑151 (MQ=255)
                                                                                                                                                  tAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggt      >  4:195/1‑151 (MQ=255)
                                                                                                                                                   aaac‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtc     <  2:933529/151‑1 (MQ=255)
                                                                                                                                                      a‑ccaccGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtctgg  >  2:799547/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
TAAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGGTCTGG  >  NC_000913/801969‑802268

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