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

AAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACC‑TGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGGTCTGG  >  NC_000913/801970‑802268
                                                                                                                                                      |                                                                                                                                                       
aaGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACC‑TGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAa‑cc                                                                                                                                                      <  1:1240964/151‑1 (MQ=255)
   cGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAAAC‑TGTCGTACTGTTAATTGCCGTTGCTGCCTCAAGGGTGGTGGTCGGT‑AACTTATCCGCCGGTGCGTTTAAAt‑ccacc                                                                                                                                                   <  2:1265980/151‑1 (MQ=255)
     aaCTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACC‑TGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCAccgc                                                                                                                                                 <  2:985393/151‑1 (MQ=255)
          gCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACC‑TGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTAt                                                                                                                                            <  2:984004/151‑1 (MQ=255)
           cGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACC‑TGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTAtt                                                                                                                                           >  1:1241601/1‑151 (MQ=255)
           cGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACC‑TGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTAtt                                                                                                                                           >  1:1429644/1‑151 (MQ=255)
            gTGGGTGTTTTTTGGTATTTCCCTGGCGGCCATGGTGGGGCTGGTTATCAAGCCTTTCCCGGAACC‑TGTCGTACTGTTAATTGCCGTTGCTCCCTCAAGGGTGGTGGTCGGT‑AACTTACCCGACGGTGCGTTAAAAA‑CCACCGCCGTAtt                                                                                                                                           <  2:1293515/150‑1 (MQ=255)
              gggTGTTTTTTGGTATTTACCTGGCGGCCTTTGTGGGGCTGGTTATCAAGCCTTCCCCGGAACC‑TGTCTTACTGTTAATTGCCGTTGCTCCCTCAATGGTGGTGGTCGGT‑AACTTACCCGCCGGTGCGTTAAAAC‑CCACCGCCGTATTa                                                                                                                                          <  2:832444/149‑1 (MQ=255)
              gggTGTTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTCCCCGGACCC‑TGTCGTACTGTTAATTGCCGTGGCTGCCTCAATGGTGGTGGTCGGT‑AACTTACCCGACGGTGCGTTTAAAC‑CCACCGCCGTAtt                                                                                                                                           <  1:1320135/148‑1 (MQ=255)
                  gCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACC‑TGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGt                                                                                                                                    >  1:563129/1‑151 (MQ=255)
                      tttGTTTTTTCCCTGGCGGCCTTTGTGGGGCTGGTTATCAAGCTTTCCCGGGACCC‑TGTCGTACTGTTAATTCCCGTTGCTCCCACAATGGTGGTGTTGGGT‑AATTTATCCGAGGGTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTAc                                                                                                                                 <  2:129249/150‑1 (MQ=255)
                        tgttttttCCCTGGGGGCTTTTGTGGGGCTGGTTATCAAG‑CTTTCCCCGAACCTTTTCGTACTTTTATTTGCTTTTGCTCCCTCAATGGTGGTGGTCGGT‑AATTTATCCGACGGTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTActc                                                                                                                               <  2:1317118/145‑1 (MQ=255)
                        tGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACC‑TGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTActct                                                                                                                              >  1:1416271/1‑151 (MQ=255)
                        tGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACC‑TGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTActct                                                                                                                              >  1:327782/1‑151 (MQ=255)
                        tGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACC‑TGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTActct                                                                                                                              >  1:228012/1‑151 (MQ=255)
                          gTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACC‑TGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTc                                                                                                                            >  1:774379/1‑151 (MQ=255)
                             tttCCCGGGCGCCCTTTGTGGGGTTGTTTACCAAGCTTTCCCCGGACCC‑GGTCTTACTGTAATTTGCCGTTGCTGCCCCAATGGGGGTGGTCGGT‑AATTTACCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCa                                                                                                                           <  2:813497/149‑1 (MQ=255)
                              ttACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACC‑TGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGt                                                                                                                        >  2:778106/1‑151 (MQ=255)
                                  cTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACC‑TGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTacc                                                                                                                     >  2:291440/1‑150 (MQ=255)
                                           aTTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACC‑TGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctggctgg                                                                                                           >  2:6908/1‑151 (MQ=255)
                                             tgtgGGGCTGGTTATCAAGCCTTTCCCGGAACC‑TGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGgtg                                                                                                         >  2:510141/1‑151 (MQ=255)
                                                ggggTTGGTTATCAAGCCTTTCCCGGACCC‑TGTCGTATTGTTAATTGCCTTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTACCCGACGGTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGtt                                                                                                       <  1:97203/150‑1 (MQ=255)
                                                  ggCTGGTTATAAAGCCTTCCCCGGACCC‑TGTCGTATTGTTAATTGCGGTTGCTGCCTCAATGGGGGTGGTGGGT‑AACTTACCCGGGGGGGTTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTtctc                                                                                                    <  2:526971/151‑1 (MQ=255)
                                                    cTGGTTATCAGGCCTTTCCCGGACCC‑TGTCGTACTGTTAATTGCCGTTGCTCCCTAAATGGGGGTGGTCGGT‑ACCTTACCCGACGGTGCTTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCg                                                                                                   <  1:342589/150‑1 (MQ=255)
                                                      ggTTATCAAGGCTTTCCCGTAACC‑TGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGCTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCATCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCg                                                                                                >  2:607565/1‑151 (MQ=255)
                                                        ttATCAAGCCTTTCCCGGAACC‑TGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGtt                                                                                              >  2:628947/1‑151 (MQ=255)
                                                             aaGCCTTCCCCGGACCC‑TGTCGTACTGTTATTGGCCGTTGCTCCCTCAATGGTGGTGGTCGGT‑AACTTACCCGACGGTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCt                                                                                         <  1:496187/151‑1 (MQ=255)
                                                                 cTTTCCCGGAACC‑TGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTa                                                                                       >  1:332493/1‑149 (MQ=255)
                                                                  tttCCCGGAACC‑TGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAg                                                                                     >  2:913607/1‑150 (MQ=255)
                                                                   ttCCCGGACCC‑TGTCTTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgc                                                                                    <  1:145492/150‑1 (MQ=255)
                                                                   ttCCCGGAACC‑TGTCGTATTGTTAATGGCCGTTGCTGCCTAAAGGGGGGTGGTCGTT‑AATTTATCCGACGGTGCGTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgc                                                                                    <  2:1376120/150‑1 (MQ=255)
                                                                   tcccgggACCC‑TGTCTTACTGTTAATTGCCGTTGGTGCCTCAAGGGTGGGGGTGGGT‑AACTTACCCGAGGGTGCGTTAAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgc                                                                                    <  1:936046/145‑1 (MQ=255)
                                                                         gaccc‑TGTCGTACTGTTAATGGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAC‑CCGCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAt                                                                             <  1:6908/148‑1 (MQ=255)
                                                                          accc‑TTTCGTCCTTTAAATGGCCGTTGTTGCCTCAATGGTGGTGGTCGTA‑AATTTACCCGCCGGTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAtt                                                                            <  1:420387/149‑1 (MQ=255)
                                                                          aaCC‑TGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAtt                                                                            >  1:814956/1‑151 (MQ=255)
                                                                                gTCGTATTGTTATTTGCCGTTGCTGCCTAAAGGGGGGTGGTCGTT‑AACTTATCCGACGGTGCGTTAAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtg                                                                          <  2:1110359/148‑1 (MQ=255)
                                                                                    tACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGaccac                                                                   >  1:1277914/1‑151 (MQ=255)
                                                                                     aCTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGaccacc                                                                  >  2:1146517/1‑151 (MQ=255)
                                                                                      cTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCg                                                                 >  1:1073185/1‑151 (MQ=255)
                                                                                         ttAATTGCCGTTGCTCCCTAAAGGGTGGTGGTCGGT‑AACTTACCCGACGGGGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGacca                                                                    <  2:670471/145‑1 (MQ=255)
                                                                                          tAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AATTTATCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGtt                                                              <  1:1079961/150‑1 (MQ=255)
                                                                                           atttGCCGTGGCTGCCTCAATGGTGGTGGTGGGT‑AACTTATCCGACGGTGCGTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTa                                                            <  1:731068/149‑1 (MQ=255)
                                                                                              ggcggtggctgcCCCTAGGGTGGTGG‑CGGTAAACTGCCCCGACGGTGCTTTAAACC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAggc                                                         <  2:608940/144‑2 (MQ=255)
                                                                                                ccGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGCCGGTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGt                                                         <  1:1315686/149‑1 (MQ=255)
                                                                                                ccGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTa                                                        >  1:269046/1‑150 (MQ=255)
                                                                                                    tgctgcCTCAGGGGGGGTGTTCGGT‑AATTTACCCGCCGGTGCTTTAAAAC‑CCGCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTaaa                                                      <  2:904262/148‑1 (MQ=255)
                                                                                                    tgctgcCTCAATGGTGGTGGTCGGT‑AATTTACCCGGCGGTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAgg                                                          <  2:332488/144‑1 (MQ=255)
                                                                                                     gctgcCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGt                                                   >  2:1165987/1‑150 (MQ=255)
                                                                                                      ctgcCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACg                                                    <  1:778121/148‑1 (MQ=255)
                                                                                                      ctgcCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAt                                                 >  1:514492/1‑151 (MQ=255)
                                                                                                        gcCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTg                                               >  1:266436/1‑151 (MQ=255)
                                                                                                         ccTAAATGGGGGTGGTGGGT‑AACTTACCCGGCGGTGCTTTTAACC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAtt                                                <  2:893895/149‑1 (MQ=255)
                                                                                                           tCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCt                                            >  1:534888/1‑151 (MQ=255)
                                                                                                             aaggggggtggtCGGT‑AACTTACCCGGCGGTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAt                                                 <  1:1384230/138‑1 (MQ=255)
                                                                                                             aaTGGTGGTGGTGGTT‑AATTTATCCGCCGGTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAtt                                                <  2:904962/145‑1 (MQ=255)
                                                                                                             aaTGGTGGTGGTCGTT‑ATCTTACCCGACGGTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTg                                               <  2:1133747/146‑1 (MQ=255)
                                                                                                             aaTGGTGGTGGTCGGT‑AACTTATCCGAGGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAt                                                 <  2:1077662/144‑1 (MQ=255)
                                                                                                                ggtggtggtCGGT‑ATTTTACCCGACGGTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATct                                        <  2:781380/150‑1 (MQ=255)
                                                                                                                ggtggtggtCGGT‑AACTTATCCGACGGTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATct                                        <  2:1242754/150‑1 (MQ=255)
                                                                                                                ggtgggggtCGGT‑AACTTACCCGAGGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCt                                            <  2:1371710/146‑1 (MQ=255)
                                                                                                                      gtccGGT‑ACTTTATCCGCCGGT‑CTTTTAAAACCCTCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATct                                        <  1:628955/141‑1 (MQ=255)
                                                                                                                         cGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGt                              >  1:1277862/1‑151 (MQ=255)
                                                                                                                         cGGT‑AACTTACCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTgg                               <  2:123517/150‑1 (MQ=255)
                                                                                                                                 ttaccggACGGTGCGTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATc                        <  2:1085651/144‑1 (MQ=255)
                                                                                                                                 ttATCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTaaaa                          <  2:1165082/148‑1 (MQ=255)
                                                                                                                                  tATCCGACGGGCCGTTAAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATc                        <  1:417950/149‑1 (MQ=255)
                                                                                                                                   aTCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGt                     >  1:1403167/1‑151 (MQ=255)
                                                                                                                                    tCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGt                     >  2:1045031/1‑150 (MQ=255)
                                                                                                                                    cccgagggcggttttAACC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTa                    <  2:1145847/139‑1 (MQ=255)
                                                                                                                                     ccaagggTGCGTTAAAAC‑CCTCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTa                    <  2:213558/145‑1 (MQ=255)
                                                                                                                                          ggTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAAcacc               <  1:204467/150‑1 (MQ=255)
                                                                                                                                          ggTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAAcacc               >  1:442154/1‑150 (MQ=255)
                                                                                                                                            tGCGTTTAAAA‑CCACCGCTGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACg            >  1:1017530/1‑151 (MQ=255)
                                                                                                                                              cGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGct          >  1:377175/1‑151 (MQ=255)
                                                                                                                                              cGTTTAAAA‑CCACCGCCATATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGct          >  1:84148/1‑151 (MQ=255)
                                                                                                                                               ttttAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACg            <  1:946516/147‑1 (MQ=255)
                                                                                                                                                   aaac‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtc     <  2:220093/151‑1 (MQ=255)
                                                                                                                                                   aaaa‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtc     >  2:1377332/1‑151 (MQ=255)
                                                                                                                                                      a‑ccaccGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtctgg  >  2:1315503/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
AAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACC‑TGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGT‑AACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGGTCTGG  >  NC_000913/801970‑802268

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