Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913802,1180AC47.6% 61.1 / 17.9 63T78P (ACC→CCC) ybhIputative DASS family tricarboxylate or dicarboxylate transporter
Reads supporting (aligned to +/- strand):  ref base A (22/11);  new base C (0/30);  total (22/41)
Fisher's exact test for biased strand distribution p-value = 3.78e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGGTCTGG  >  NC_000913/801977‑802268
                                                                                                                                              |                                                                                                                                                       
cTGGCGTGGGGGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCACGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCAccgcc                                                                                                                                                >  2:758266/1‑150 (MQ=255)
   gCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTAt                                                                                                                                            >  1:135790/1‑151 (MQ=255)
      tGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTa                                                                                                                                          >  1:881928/1‑150 (MQ=255)
          tGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCgg                                                                                                                                     >  2:307805/1‑151 (MQ=255)
             tttttGTTATTTCCCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTCCCCGGAC‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAAGGGGGGTGGTCGGTAATTTACCCGACGGTGCGTTAAAAC‑CCACCGCCGTATTAAGCGGtt                                                                                                                                   <  2:54790/150‑1 (MQ=255)
              tttgGGTATTTCCCTGGCGGCCTTTGTGGGGCGGGTTATCAA‑CCTTTCCCCGAACCCGGTCTTATTGTTAATGGCCTTTGCTGCCTCAATGGTGGGGGTCGGTAATTTACCCGCCGGTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGtt                                                                                                                                   <  1:234232/149‑1 (MQ=255)
                 tggtttttccgggggggcAATTGTGGGCGTGTTAACCAAGCTTTTCCCGGAA‑CTTTTCGTCTTTTTAATTCCCGTTGCTCCCTCAATGGTGTTGGTGGTTAACTTAGCCGACGGGGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTActc                                                                                                                               <  2:917679/135‑1 (MQ=255)
                 ggttttttCCCGGGGGGCCTTTGGGGGGCTTGTTATAAGGCCTTCCCCGGAC‑CTTTTCGTACTTTTAATTGCCGTTCCTCCCTAATGGGTGGTGGTGGGTACTTTTCCCGACGGTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTActc                                                                                                                               <  2:1171811/145‑1 (MQ=255)
                    tATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTc                                                                                                                            >  1:115864/1‑150 (MQ=255)
                       ttcccgggcggcCATTGTGGGGCTGTTTAACAAGCTTTCCCCGGAA‑CCTGTCGTACTGTTAATGGCCGTTGCTGCCTAAATGGGGGGGGTGGGTAACTTACCCGACGGTGCGTTAAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAgg                                                                                                                         <  1:155165/144‑1 (MQ=255)
                       ttACCTGGCGGCGTTTGGGGGGTTGTTTATCAAGCTTTTCCCGGAC‑CCTTTGGTATTGTAAATTGCTGTTGCTGCCTAAAGGGTGGTGGTCGTTACCTTACCGGACGTTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAgg                                                                                                                         <  1:582837/150‑1 (MQ=255)
                          ccgggcggcCTTTGTGGGGCTGGTTATCAAGCTTTCCCCGGAC‑CCTGTCGTACTGTTAATTGCCGTTGCTCCCTCAATGGTGGTGGTCGGTAATTTACCCGACGGGGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTa                                                                                                                       <  1:1116233/146‑1 (MQ=255)
                          ccTGGCGCCCTTTTGGGGGCTGTTTAACAAGCTTTTCCGGAAC‑CCTGTCGTCTTGTTAATTGCGTTTGCTGCCCCAATGGTGGTGTTCGGTAATTTACCCGCCGGTGTTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGt                                                                                                                        <  1:574692/148‑1 (MQ=255)
                            tGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTaccac                                                                                                                   >  2:882617/1‑151 (MQ=255)
                             ggcggcCTTTGTGGGGCTGGTTATCAAGCCTTTCCCGGAA‑CCTTTCGTACTTTTAATTGCCGTTGCTGCCTAAAGGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTaccac                                                                                                                   <  1:1048852/150‑1 (MQ=255)
                                       gtgGGGCTGGTTATCAAGCTTTCCCCGGAC‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAATTTACCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGgt                                                                                                          <  1:256230/149‑1 (MQ=255)
                                       gtgGGGCTGGTTATAAACCTTTTCCCGGAC‑CCTGTCTTATTTTTATTTGCCTTTGCTGCCTCAATGGTGGTGGTCGGTACCTTACCCGACGGTGCTTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGgtgt                                                                                                        <  1:996349/151‑1 (MQ=255)
                                        tgGGGCTGGTTATCAAGCCTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTCCCACCTGGCTGGTGtt                                                                                                       >  2:242876/1‑151 (MQ=255)
                                         ggggTTGGTTATCAAGCTTTCCCCGGAA‑CTTGTCTTATTGTTAATTGCCGTTGCTGCCTAAATGGTGGTGGTCGGTAATTTATCCGACGGTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGtt                                                                                                       <  1:135459/150‑1 (MQ=255)
                                              tGGTTATCAAGCCTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGACGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTtctc                                                                                                    >  2:795645/1‑148 (MQ=255)
                                                   aTCAAGCCTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGttt                                                                                             >  2:399578/1‑150 (MQ=255)
                                                      atgcttttCCGAGAT‑CCTGTCGTATTGTTAATTACCTTCGCTGCCACAATGGTGGTGGTGGAAATCTTAGCCGTGGGTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTAcc                                                                                          <  2:268800/145‑1 (MQ=255)
                                                      aatccTTTCCCGAAA‑CCGGCCGTTTTTTAAAGCGCCGTTGCGGCCTCATTGGGGTGGGTCGGTAACTTATCCGACGGTGCTTTAAAAC‑CCTCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGttt                                                                                             <  1:117902/144‑1 (MQ=255)
                                                         ccTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTa                                                                                       >  2:133098/1‑150 (MQ=255)
                                                           ttccccGGAA‑CCTTTCGTTCTTTAAATTGCCTTTGCTCCCTCAATGGTGGTGGTGGGAAACTTACCCGCCGGTGCTTTAAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCtt                                                                                        <  1:576432/144‑1 (MQ=255)
                                                             tCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgcgc                                                                                  >  1:385241/1‑151 (MQ=255)
                                                              cccGGAC‑CTTGTCGTACTTTTAATTCCCTTTGCTCCCTAAAGGGTGGTGGTCGGTAACTTATCCGACGGGGTTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCtt                                                                                        <  2:219045/144‑1 (MQ=255)
                                                              cccGGAC‑CCTGTGGTACTTTTATTTGCCTTTGCCGCCTCAAGGGTGGTGGTGGGTATTTTATCGGCCGGTGTTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgc                                                                                    <  2:1058433/148‑1 (MQ=255)
                                                               ccGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGcgccg                                                                                >  2:584112/1‑151 (MQ=255)
                                                                    acccTTTCGTACTGTTATTGGCCTTTGCTGCCACAATGGGGGTGGTCGGTATCTTATCCGAGGGTCTTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCa                                                                              <  1:539995/147‑1 (MQ=255)
                                                                       cTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtgt                                                                         >  1:849749/1‑151 (MQ=255)
                                                                            gTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGacc                                                                     >  1:501278/1‑150 (MQ=255)
                                                                             tACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGaccac                                                                   >  1:726645/1‑151 (MQ=255)
                                                                              actttaatttGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAATTAACCCGCCGGTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGaccac                                                                   <  1:509861/142‑1 (MQ=255)
                                                                                tGTAAATGGCCGTTGCGCCCTAAAGGGGGGGGTGCGGTAACTTTCCCCACGTTGCTTTAAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGacca                                                                    <  2:236949/147‑1 (MQ=255)
                                                                                         ccGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTa                                                        >  2:968003/1‑150 (MQ=255)
                                                                                         ccGTTGCTGCCTCAAGGGTGGGGGTCGTTAACTTACCCGACGGTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGt                                                         <  1:290163/149‑1 (MQ=255)
                                                                                            ttGCTCCCTCAAGGGGGTTGGTCGGTAATTTATCCGGGGGTGCTTTAAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAAc                                                     <  2:583930/150‑1 (MQ=255)
                                                                                             tgctgcCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGt                                                   >  1:241880/1‑151 (MQ=255)
                                                                                                tgcCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAtt                                                >  1:474396/1‑151 (MQ=255)
                                                                                                 gcCTAAAGGGGGGTGTTGGTTAACTTATCCGACGGTGCGTTAAAAA‑CCTCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGTCCACCGGTTTAGGTAAAc                                                     <  1:627812/145‑1 (MQ=255)
                                                                                                  ccTAAATGGTGGTGGTGGGTACTTTCTCCGAGGGGGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAtt                                                <  2:387393/149‑1 (MQ=255)
                                                                                                   ctaattggtgggggTGGTTAACTTACCCGAGGGTGCGTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGt                                                   <  1:897521/140‑1 (MQ=255)
                                                                                                      aaTGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTAt                                          >  2:627420/1‑151 (MQ=255)
                                                                                                      aaTGGTGGTGGTCGGTAACTTACCCGACGGTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTa                                           <  1:201588/150‑1 (MQ=255)
                                                                                                         ggtggtggtCGGTATCTTACCCGACGGTGCGTTAAAAC‑CCGCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATct                                        <  1:882677/150‑1 (MQ=255)
                                                                                                         gggggtggtCGTTATCTTACCCGACGGTGGTTTAAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctg                                       <  2:1154830/148‑1 (MQ=255)
                                                                                                               ggttggtttcTTACCCGAGGGTGGTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgc                                      <  2:1105178/137‑1 (MQ=255)
                                                                                                                gtgggtattttACCCGGCGGTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTg                                <  1:584140/142‑1 (MQ=255)
                                                                                                                  cGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGt                              >  2:862028/1‑151 (MQ=255)
                                                                                                                      aaCTTTCCCGTCGGTGCTTAAAAAC‑CCCCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGAt                                  <  2:1177818/143‑1 (MQ=255)
                                                                                                                         ttACCCGACGGTGTTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGt                              <  1:339756/144‑1 (MQ=255)
                                                                                                                          ttacccgcgggtgttttAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATc                        <  2:209215/136‑1 (MQ=255)
                                                                                                                           acccGACGGTGCGTTAAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCgg                      <  1:806405/148‑1 (MQ=255)
                                                                                                                           aTCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCgg                      >  1:114505/1‑150 (MQ=255)
                                                                                                                             ccgacgtgtcGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTa                    <  2:137733/141‑1 (MQ=255)
                                                                                                                             ccGAGGGTGTTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTaa                   <  2:1169137/151‑1 (MQ=255)
                                                                                                                             ccGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGt                     <  2:455227/149‑1 (MQ=255)
                                                                                                                             ccGACGGTGCGTTTAAAA‑CCGCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTa                    <  2:778510/150‑1 (MQ=255)
                                                                                                                                aCGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAacac                >  2:1139694/1‑151 (MQ=255)
                                                                                                                                  ggTGCGTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAAcacca              <  2:447237/151‑1 (MQ=255)
                                                                                                                                     gctttgaaac‑CCCCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTaa                   <  1:793139/137‑1 (MQ=255)
                                                                                                                                      gctttAAAACCCTCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACg            <  2:1168071/148‑1 (MQ=255)
                                                                                                                                        tttAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACg            <  2:201313/147‑1 (MQ=255)
                                                                                                                                          aaaac‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggt      <  2:683213/150‑1 (MQ=255)
                                                                                                                                           aaaa‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtc     >  1:394129/1‑151 (MQ=255)
                                                                                                                                              a‑ccaccGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtctgg  >  1:1055368/1‑151 (MQ=255)
                                                                                                                                              |                                                                                                                                                       
CTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGGTCTGG  >  NC_000913/801977‑802268

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