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

CGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GG‑AACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGGTCTG  >  NC_000913/801973‑802267
                                                                                                                                                   |                                                                                                                                                      
cGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GG‑AACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAa‑ccacc                                                                                                                                                  <  1:1067586/151‑1 (MQ=255)
 gAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GG‑AACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCCATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCAccg                                                                                                                                                 >  1:376976/1‑151 (MQ=255)
       gCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GG‑AACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTAt                                                                                                                                           >  2:906916/1‑151 (MQ=255)
       gCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GG‑AACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTAt                                                                                                                                           <  2:1181855/151‑1 (MQ=255)
             gTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GG‑AACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCg                                                                                                                                     >  2:1130418/1‑151 (MQ=255)
                cTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GG‑AACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGtt                                                                                                                                  >  1:632674/1‑151 (MQ=255)
                    ttGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GG‑AACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTActc                                                                                                                              >  1:1189830/1‑151 (MQ=255)
                    ttGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GG‑AACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTActc                                                                                                                              >  1:129586/1‑151 (MQ=255)
                     tGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GG‑AACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTActct                                                                                                                             >  1:888270/1‑151 (MQ=255)
                        tttttCCCTGGGGGCCTTTGTGGGGCTGGTTATCAAGCTTTTCCG‑GA‑ACCCTGTCGTCCTGTTAATTGCCGTTGCTCCCTAAATGGTGGTGGTCGGTAATTTACCCGAGGGTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCtt                                                                                                                            <  2:560023/147‑1 (MQ=255)
                         aTTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GG‑AACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTc                                                                                                                           >  1:629010/1‑149 (MQ=255)
                           ttACCTGGCGGACATTGTGGGGATGGTTATCAAGCATTTCCC‑GG‑AACATGTCGTACTGTTAATTGACGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGAGCGTTTAAAA‑CCACCGTCGTATTAAGCGGTAActct                                                                                                                             >  1:847182/1‑145 (MQ=255)
                           ttACCGGGCGGCCTTTGTGGGGCTGGTTATCAACCCTTCCCC‑GG‑AACCTGTCGTTTTGTTAATTCCGGTTTCTGCCTCATTGGGGGTGTTGGGTAACTTATCCGCCGGTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAgg                                                                                                                        <  1:1276115/150‑1 (MQ=255)
                              ccgggcggcCTTTGGGGGGCTGGTTATCAAGCCTTTCCC‑GG‑ACCCTGTCGTACTTTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGt                                                                                                                       <  2:308013/145‑1 (MQ=255)
                              ccgggcggcCTTTGGGGGGCTGGTTACCAAGCCTTCCCC‑GG‑ACCCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTACCCGACGGTGCGTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTa                                                                                                                      <  2:1268356/146‑1 (MQ=255)
                              ccTGGCGGCGATTGGGGGGCTGTTTATCAACCCTTCCCG‑GG‑AACCTGTCGTACTTTTAATTGCTGTTGCCCCCTCAAGGGTGGTGGTGGGTAATTTATCCGGGGGTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTa                                                                                                                      <  1:895614/149‑1 (MQ=255)
                                 ggcgccCTTTGGGGGGCTGGTTATCAACCTTTTCCC‑GG‑ACCCTGTCGTACTGTTATTTGCCGTTGCTCCCTCATTGGGGTTGGTGGTTAACTTATCCGACGGTGCTTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTacc                                                                                                                    <  2:1012624/148‑1 (MQ=255)
                                    ggcCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GG‑AACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctgg                                                                                                              >  2:516960/1‑151 (MQ=255)
                                        aTTGTGGGGCTGGTTATCAAGCCTTTCCC‑GG‑AACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctggctgg                                                                                                          >  2:400741/1‑151 (MQ=255)
                                           gtgGGGCGGTTTAGCACTCTTTCCCC‑GA‑ACCCTGTCGTACTTTTATTTGCCGTTGCTGCCTCAAGGGTGGTGGTCGGTATCTTCCCCGACGTTGCGTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGgtg                                                                                                        <  1:1254773/150‑1 (MQ=255)
                                                gCTGGTTATCAAGCTTTCCCC‑GG‑ACCCTGTCGTATTGTTAATTGCCGTTGCTGCCTAAATGGTGGTGGTCGTTAACTTACCCGACGGTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTtctc                                                                                                   <  2:1012670/150‑1 (MQ=255)
                                                     ttATCAAGCCTTCCCG‑GG‑AACTTTTCGTATTGTTAATGGCCGTTGCTGCCTAAGTGTGGGTGGTCGGTAATTTATCCGAGGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCg                                                                                                  <  1:1132434/146‑1 (MQ=255)
                                                       aTCAAGCCTTTCCC‑GG‑ACCCTGTCTTACTTTTATTGCCCGTTGCTGCCTCAATGGTGGTGGTCGGTAATTTATCCGGCGGTGCGTTAAAAC‑CCGCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGtt                                                                                             <  1:715582/149‑1 (MQ=255)
                                                       aTCAAGCCTTTCCC‑GG‑ACCCTGTCGTCTTTTTAATTGCCTTTGCTCCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTa                                                                                           <  2:554859/151‑1 (MQ=255)
                                                       aTCAACCCTTCCCC‑GG‑ACCCTGTCGTACTTTTAATTCCCGTTGCTGCCTCAATGGTGGGGGTCGGTACCTTACCCGCGGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTa                                                                                           <  2:475188/151‑1 (MQ=255)
                                                        tCAAGCCTTTCCC‑GG‑AACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTAc                                                                                          >  2:203618/1‑151 (MQ=255)
                                                          acgCCTTCCCC‑GG‑ACCCTGTCGTACTTTTAATTCCCGTTGCGGCCTCAGGGGTGGTGGTCGGTACCTTACCCGCCGGTGCTTTAAAAC‑CCGCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTAc                                                                                          <  2:1237605/147‑1 (MQ=255)
                                                              cTTTCCC‑GGAAACCTGTCGTTTTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGTTAACTTACCCGGGGGTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCtt                                                                                       <  1:891586/149‑1 (MQ=255)
                                                              cTTTCCC‑GG‑AACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTa                                                                                      >  2:1153668/1‑149 (MQ=255)
                                                                ttCCCGGG‑AACCTGTCGTACTTTTATTTGCCGTTGCTGCCTCAAGGGGGGTGGTGGGTATCTTCTCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTaa                                                                                     <  1:846484/149‑1 (MQ=255)
                                                                  ccc‑GG‑CACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGCGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGcgcc                                                                                >  1:397325/1‑151 (MQ=255)
                                                                   cc‑GG‑ACCCGGTCGTCCTGTTAAGTGCCGTTGCTGCCTCAAGGGTGGGGGTGGGTACTTTGCCCGGGGGTGCTTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGcgccg                                                                               <  2:888929/151‑1 (MQ=255)
                                                                   cc‑GG‑AACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGcgccg                                                                               >  2:700943/1‑151 (MQ=255)
                                                                         accctgccttccttttAATTGCCGTTGCTCCCTCAAGGGGGGTGGTCGGTAACTTACCCGACGGTGCGTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAt                                                                            <  2:436761/136‑1 (MQ=255)
                                                                             tggcgtggttttAATAGCTGTTCCTCCCTCAATGGTGGTGGTCGGTAACTTCGCCGGCGGTGTGTTAAAAC‑CCTCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAttt                                                                          <  2:218318/138‑1 (MQ=255)
                                                                              gtcgtttttttATTTGCCGTTGCGCCCTAAGTGGTGGTGGTCTTTATCTTAGCCGTGGGTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtgt                                                                        <  2:131365/140‑1 (MQ=255)
                                                                                cGTACTTTTATTGGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGac                                                                     <  2:400844/150‑1 (MQ=255)
                                                                                 gTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTACCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGacc                                                                    >  2:18163/1‑150 (MQ=255)
                                                                                  tACTGTTAATTGCCGTTGCTGCCTCAAGGGTGGTGGTCGGTAATTTATCCGCCGGTGCTTTAAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtgtg                                                                       <  2:164863/146‑1 (MQ=255)
                                                                                       ttAATTGCCGTTGCTGCCTCATTGGTGGTGGTGGTTAATTTACCCGACGTTGGTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCgg                                                               <  1:1213508/149‑1 (MQ=255)
                                                                                       ttAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTACCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCgg                                                               <  2:946782/149‑1 (MQ=255)
                                                                                        tAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGttt                                                            >  1:764595/1‑151 (MQ=255)
                                                                                        tAAGTGCCGTTGCTCCCTAAGTGGGGGTGGTCGGTATCTTTCCCAACGGTGGTTTTAAAC‑CCCCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCTCCGGt                                                              <  1:327693/149‑1 (MQ=255)
                                                                                         aaTTGCCTTTGCTGCCTCATGGGGGGGGGTCGTTACCTTATCCGCCGGTGCTTTAAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGtt                                                             <  2:65681/149‑1 (MQ=255)
                                                                                         aaTTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTACCTTATCCGACGGTGCGTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTa                                                           <  2:131376/151‑1 (MQ=255)
                                                                                           ttGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAACGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAgg                                                         >  1:343761/1‑151 (MQ=255)
                                                                                              ccGTTGCTGCCTCATTGGTGGGGGTTGGTAACTTACCCGACGGTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTa                                                       <  2:1158976/150‑1 (MQ=255)
                                                                                              ccGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTACCCGACGGTGCGTTAAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGt                                                        <  2:1224441/149‑1 (MQ=255)
                                                                                               cGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTaa                                                      >  1:624126/1‑150 (MQ=255)
                                                                                                   gctgcCTCAATGGTGGTGGGCGGTAACTTATCCGACGGTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGt                                                  <  2:970111/150‑1 (MQ=255)
                                                                                                     tgcCCCAGTGGGGGGGGTCGGTAATTTCTCCTACGGT‑CGTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCTCCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTa                                                       <  2:1146562/143‑1 (MQ=255)
                                                                                                        cTCAAGGGTGGTGGTCGTTACTTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAt                                                <  2:955538/147‑1 (MQ=255)
                                                                                                         tCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCt                                           >  2:282756/1‑151 (MQ=255)
                                                                                                           atgggtggtggtCGTAAGTTTATCCGGCGGTGCGTTTAAAG‑CCCCCGCCGTATTAAGCGGTTACTCTTCAGGTTCCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTa                                                 <  1:1258756/140‑1 (MQ=255)
                                                                                                           aagggtggtggtCGGTAACTTATCCGAGGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGcc                                            <  2:717838/145‑1 (MQ=255)
                                                                                                           aaTGGGGGTGGTCGGTAACTTACCCGGCGGTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTa                                          <  2:806333/150‑1 (MQ=255)
                                                                                                              gggggtggtCGGTATCTTACCCGACGGTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATct                                       <  2:392916/147‑1 (MQ=255)
                                                                                                              gggggtggtCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATct                                       <  2:402726/147‑1 (MQ=255)
                                                                                                                tggtggtCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgct                                    >  1:949201/1‑151 (MQ=255)
                                                                                                                 ggtggtCGGTAACTTATCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgctg                                   <  2:757777/151‑1 (MQ=255)
                                                                                                                   tggtCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGAt                                 >  1:326157/1‑151 (MQ=255)
                                                                                                                       cGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGt                             >  1:1202101/1‑151 (MQ=255)
                                                                                                                        ggtaatttccccgggggTGCTTTTAAAA‑CCCCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTa                            <  2:503682/136‑1 (MQ=255)
                                                                                                                        ggTATCTTACCCGAGGGTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGt                             <  2:1225463/150‑1 (MQ=255)
                                                                                                                           aaCTTATCCGCCGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTa                            <  2:997034/148‑1 (MQ=255)
                                                                                                                           aaCTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTaaa                          >  1:117583/1‑150 (MQ=255)
                                                                                                                              ttATCCGCGGGTGTTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATc                       <  1:601774/150‑1 (MQ=255)
                                                                                                                                aTCCGACGGTGCGTTCAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGt                    >  2:843218/1‑151 (MQ=255)
                                                                                                                                  ccggggttGCGTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAAt                        <  1:217022/138‑1 (MQ=255)
                                                                                                                                        gTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAAcaccac            >  2:170573/1‑151 (MQ=255)
                                                                                                                                         tGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACg           >  2:547426/1‑151 (MQ=255)
                                                                                                                                          gCGTTAAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAAcaccac            <  2:761562/149‑1 (MQ=255)
                                                                                                                                              ttAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctggg      >  1:1190563/1‑151 (MQ=255)
                                                                                                                                               tAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgg       >  1:586220/1‑149 (MQ=255)
                                                                                                                                                aaac‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtc    <  2:415752/151‑1 (MQ=255)
                                                                                                                                                 aaa‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtct   >  1:31972/1‑151 (MQ=255)
                                                                                                                                                  aa‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtctg  >  2:607183/1‑151 (MQ=255)
                                                                                                                                                  aa‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtctg  >  1:423589/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
CGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCC‑GG‑AACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGGTCTG  >  NC_000913/801973‑802267

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