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

TAAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCTGTCGTACTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGGTCTGG  >  NC_000913/801969‑802268
                                                                                                                                                        |                                                                                                                                                       
tAAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCTGTCGTACTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTaaaa                                                                                                                                                         >  1:444803/1‑150 (MQ=255)
   gCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCTGTCGTACTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAa‑ccac                                                                                                                                                    <  1:393847/151‑1 (MQ=255)
    cGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCTGTCGTACTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAa‑ccacc                                                                                                                                                   >  1:480206/1‑151 (MQ=255)
           gCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCTGTCGTACTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTAt                                                                                                                                            <  1:2357/151‑1 (MQ=255)
               gggTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCTGTCGTACTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTa                                                                                                                                          >  2:411032/1‑149 (MQ=255)
                    cTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCTGTCGTACTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGtt                                                                                                                                   >  2:723499/1‑151 (MQ=255)
                       tttGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCTGTCGTACTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTAct                                                                                                                                >  2:528318/1‑151 (MQ=255)
                            tATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCTGTCGTACTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTc                                                                                                                            >  1:5058/1‑150 (MQ=255)
                             aTTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCTGTCGTACTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTc                                                                                                                            >  2:118035/1‑149 (MQ=255)
                              tttCCCTGGCGGCCTTGGTGGGGTTGTTTATAAAGCCTTT‑CCCGGACCTTGTCGTCCTGTTAATTGCCGTT‑GCTCCCTCAATGG‑TGGTGGTCGGTAACTTATCCGCCGGTGCGTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAg                                                                                                                          <  1:194503/150‑1 (MQ=255)
                              tttCCCTGGCGCCCTTTGGGGGGCTGTTAATC‑AGCCTTTCCCCGGACCCTGTCGTCTTTTTAATTTCTGTT‑TCTGCCCCATTGG‑TGGTGGTCGGTAACTTATCCGAGGGTG‑GTTTAAAACCCCCCGCCGTATTAAGCGGTTACTCTTCa                                                                                                                           <  2:504189/149‑1 (MQ=255)
                                tACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCTGTCGTACTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGt                                                                                                                        >  1:215943/1‑150 (MQ=255)
                                           cATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCTGTCGTACTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctggctg                                                                                                            >  1:827427/1‑151 (MQ=255)
                                             ttGGGGGGTTGGTTATCAACCTTTT‑CCCGAACCCTGTCTTACTTTTAATTGCCGTT‑GCTGCCTCAATGG‑GGGTGGTCGGTACCTTACCCGACGGGGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctggctgg                                                                                                           <  1:910946/150‑1 (MQ=255)
                                              tgtgGGGCTGGTTATCAAGCCTTT‑CCCGGAACCTGTCGTACTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGTGGTTACTCTTCAGGTACCACCTGGCTAgtg                                                                                                         >  2:21630/1‑151 (MQ=255)
                                               gtgGGGCTGGTTATCAAGCTTTT‑CCCGGAACCTGTCGTACTGTTAATTGCCGTT‑GCTGCCTCAAGGG‑TGGTGGTCGGTAATTTATCCGCCGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGgtg                                                                                                         <  1:319811/150‑1 (MQ=255)
                                               gtgGGGCTGGTTATAAACCTTTT‑CCCGGATCCTGTCGTATTTTTAATTCCCTTT‑GCTCCCTCATTGG‑TGGTGGTGGTTACTTTACCCGGCGTTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGgtg                                                                                                         <  1:281967/150‑1 (MQ=255)
                                                     cTGGTTATCAAGCCTTT‑CCCGGAACCTGTCGTACTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCgg                                                                                                  >  2:143910/1‑151 (MQ=255)
                                                     cTGGTTATCAAGCCTTT‑CCCGGAACCTGTCGTACTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCgg                                                                                                  >  1:1041576/1‑151 (MQ=255)
                                                      tGGTTATCAAGCCTTT‑CCCGGACCCTGTCGTTCTTTTAATTGCCGTT‑GCGGCCT‑AATGGTTGGTGGTGGGTAATTTTTCCGCGGTTGCGTTTAAAG‑CTGCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCg                                                                                                   <  2:697163/149‑1 (MQ=255)
                                                       ggTTATCAAGCCTTT‑CCCGGAACCTGTCGTACTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCg                                                                                                >  2:14827/1‑151 (MQ=255)
                                                       ggTTATCAAGCCTTT‑CCCGGAACCTGTCGTACTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCg                                                                                                >  1:684112/1‑151 (MQ=255)
                                                             caatcattt‑CCCGAACCCTTTCATACTGTTAATTGCCGTT‑GCTGCCACAATGG‑TGGTGGTCGTTAATTTATCCGAGGGTGCGTTTAAAA‑CCCCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTAc                                                                                           <  2:1072038/144‑1 (MQ=255)
                                                              acgctttt‑CCCGGAACCGGTCGTACTGTTAATTGCCGTT‑GCTCCCTCAATGG‑TGGGGGTCGGTAATTTATCCGAGGGTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCt                                                                                         <  2:871697/146‑1 (MQ=255)
                                                              accccTTT‑CCCGGACCCTGTCTTCTTGTTATTTGCCGTT‑GCTGCCTCAAGGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGtt                                                                                              <  2:469898/143‑1 (MQ=255)
                                                              aaGCTTTT‑CCCGGAACTTGTCGTACTGTTAATTGCCTTT‑TCTGCCTCAATGG‑TGGTGGTCGGTAATTTATCCGAGGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCggggaa                                                                                           <  1:542045/149‑6 (MQ=255)
                                                              aaGCCTTT‑CCCGGAACCTGTCGTACTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTAcc                                                                                          <  1:991888/150‑1 (MQ=255)
                                                                 ccTTT‑CCCGGAACCTGTCGTACTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTa                                                                                       >  1:396830/1‑150 (MQ=255)
                                                                   ttt‑CCCGGACCCTTTCTTACTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTaa                                                                                      <  1:527309/149‑1 (MQ=255)
                                                                    ttccccGGACCCTGTCGTACTTTTAATTGCCGTT‑GCTGCCTAAATGG‑TGGTGGTCGGAATCTTACCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTaa                                                                                      <  1:991156/146‑1 (MQ=255)
                                                                          ggAACCTGTCGTACTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGcgccgc                                                                               >  1:226068/1‑150 (MQ=255)
                                                                            acccTGTCGTACTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGGGGTCGGTAATTTACCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAtt                                                                            <  1:100711/149‑1 (MQ=255)
                                                                              ccTGTCGTACTTTTAATTGCCTTT‑GCTGCCTCAATGG‑TGGTGGTCGTTAATTTATCCGACGGTGCGTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAttt                                                                           <  2:640456/150‑1 (MQ=255)
                                                                              ccTGTCGTACTGTTAATTGCCTTT‑GCTCCCTCAAGGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAt                                                                             <  1:91050/148‑1 (MQ=255)
                                                                                  tCGTACTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGac                                                                      >  2:293576/1‑151 (MQ=255)
                                                                                   cgttttGTTAATTGCCGTTCTCTCCCT‑AGTGG‑TGGTGATCGATATTTTACCCGATGTTTATTTAAAAA‑CCCCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGac                                                                      <  2:811491/145‑1 (MQ=255)
                                                                                    gtattttaaaTTGCCGTT‑GCTGCCGAATTGG‑TGGTGGTCGATAATTTACCCGAGGGTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGacc                                                                     <  1:231212/142‑1 (MQ=255)
                                                                                     tACTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGTTCGGTAACTTACCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGacc                                                                     <  2:15349/149‑1 (MQ=255)
                                                                                       cTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCg                                                                 >  2:745134/1‑151 (MQ=255)
                                                                                           tAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGttt                                                             >  1:255414/1‑151 (MQ=255)
                                                                                           tAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGttt                                                             >  1:406704/1‑151 (MQ=255)
                                                                                                gCCGTT‑GCTGCCTCAATGG‑TGGGGGTCGGTAATTTATCCGAGGGTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGt                                                         <  1:28786/150‑1 (MQ=255)
                                                                                                    tt‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACg                                                    >  1:185022/1‑151 (MQ=255)
                                                                                                        ctgcCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAt                                                 >  1:461375/1‑151 (MQ=255)
                                                                                                           ccTCAATGG‑TGGTGGTCGGTAACTTATCCGCCGGTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTa                                                  <  2:342127/147‑1 (MQ=255)
                                                                                                           ccTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGc                                              >  1:433831/1‑151 (MQ=255)
                                                                                                               aaTGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTAt                                          >  2:416219/1‑151 (MQ=255)
                                                                                                               aaTGG‑TGGTGGTCGGTAACTTATCCGACGGGGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCt                                            <  1:1057372/149‑1 (MQ=255)
                                                                                                                         ggtCGGTATTTTATCCGACGGTGCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGAtt                                 <  2:684109/151‑1 (MQ=255)
                                                                                                                         ggtCGGTACCTTACCCGACGGTGCGTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGa                                   <  1:227995/149‑1 (MQ=255)
                                                                                                                            cGGTAACTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGt                              >  2:671600/1‑151 (MQ=255)
                                                                                                                              gTAACTTACCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTg                                <  1:675419/147‑1 (MQ=255)
                                                                                                                                aaCTTATCCGACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTaa                            <  2:491768/149‑1 (MQ=255)
                                                                                                                                 aCTTATCCGACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAAt                         >  2:141892/1‑151 (MQ=255)
                                                                                                                                 aCTTATCCGACGGAACGTTTAAAA‑CCACCGCCGTATTAAGCGGCTACTCTTCAGCTACCACCTGGCTGCCGCTCTCGGCGTTTACCTTACGTGCAGCCTTTGTGCCCACCCGTCTAGGTAAACGTATAACCTATCTGCTGATTGGt                              >  2:641734/1‑146 (MQ=255)
                                                                                                                                   ttttCCGTCGGTGCGTTTAAAA‑CCTCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTaaaa                          <  1:916768/145‑1 (MQ=255)
                                                                                                                                   ttagcggaaggTTCGTTAAAAC‑CCCCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATc                        <  2:849951/141‑1 (MQ=255)
                                                                                                                                       ccgccgGTGCGTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGt                     <  1:951511/149‑1 (MQ=255)
                                                                                                                                       ccGAGGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTa                    <  2:1015324/150‑1 (MQ=255)
                                                                                                                                         gACGGTGCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGCACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAac                  <  1:206473/150‑1 (MQ=255)
                                                                                                                                         gACGGTGCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAac                  >  1:145675/1‑150 (MQ=255)
                                                                                                                                            ggTGCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAAcacca              <  1:411034/151‑1 (MQ=255)
                                                                                                                                                cGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGct          >  2:843026/1‑151 (MQ=255)
                                                                                                                                                    tAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggt      >  1:721535/1‑151 (MQ=255)
                                                                                                                                                       aa‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtctg   >  2:489042/1‑151 (MQ=255)
                                                                                                                                                        a‑ccaccGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtctgg  >  1:842943/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                       
TAAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCTGTCGTACTGTTAATTGCCGTT‑GCTGCCTCAATGG‑TGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAA‑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 < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑