Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,707,4410CA26.3% 110.2 / 10.4 58A121S (GCG→TCG) rseCprotein RseC
Reads supporting (aligned to +/- strand):  ref base C (34/8);  new base A (0/15);  total (35/23)
Fisher's exact test for biased strand distribution p-value = 2.22e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACG‑T‑CCTCGCATTTGTTATGCAAAATGCAACAAAGCCA‑GTGAAAT‑CACTGGCTCGCGTCTTCCGAAG‑ATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAACACGCGG  >  NC_000913/2707193‑2707691
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                          
gTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACG‑T‑CCTCGCATTTGTTATGCAAAATGCAACAAAGCCA‑GTGAAAT‑CACTGGCTCGCGTCTTCCGAAG‑ATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCa                                                                                                                                                                                                                                                           >  1:300142/1‑250 (MQ=255)
         aaTAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACG‑T‑CCTCGCATTTGTTATGCAAAATGCAACAAAGCCA‑GTGAAAT‑CACTGGCTCGCGTCTTCCGAAG‑ATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAAc                                                                                                                                                                                                                                                  <  1:79705/250‑1 (MQ=255)
            aGAAACTCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCCGGAAACG‑T‑CCTCGCATTTGTTATGCAAAATGCAACAAAGCCA‑GTGAAAT‑CACTGGCTCGCGTCTTCCGAAG‑ATGCTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGACAGGCTGCCATTCTGCCCGGGCAGCAAACTTACG‑CCAGAAAGCC‑CGCGGAATCAGGAACCCGc                                                                                                                                                                                                                                              >  2:372804/1‑251 (MQ=255)
               aaCGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACG‑T‑CCTCGCATTTGTTATGCAAAATGCAACAAAGCCA‑GTGAAAT‑CACTGGCTCGCGTCTTCCGAAG‑ATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCa                                                                                                                                                                                                                                            >  1:441916/1‑250 (MQ=255)
                             cactacaACATGAAGGGGAGAAAATGTTCATACCGTATGGATGTTGGTATC‑GGAAACCCT‑CCTCGCTTTTTTAAGGCAAAATGAAAAAAAGCCA‑GAGAAAT‑CACTGGCCCGCGTCTTCAGAAG‑ATGTTTCAAACCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGAGGATCGGTTGCCATTCGCCCCGGGAAGCAACCTTACG‑GGAGTA‑GCCGCGCGAAATCAGGACCCCGCCAATACCACCGAGGATc                                                                                                                                                                                                                             <  2:441916/251‑1 (MQ=255)
                                     caTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACG‑T‑CCTCGCATTTGTTATGCAAAATGCAACAAAGCCA‑GTGAAAT‑CACTGGCTCGCGTCTTCCGAAG‑ATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGcaccac                                                                                                                                                                                                                      >  1:483259/1‑250 (MQ=255)
                                      aTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACG‑T‑CCTCGCATTTGTTATGCAAAATGCAACAAAGCCA‑GTGAAAT‑CACTGGCTCGCGTCTTCCGAAG‑ATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACAt                                                                                                                                                                                                                    >  2:174494/1‑251 (MQ=255)
                                         gAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACG‑T‑CCTCGCATTTGTTATGCAAAATGCAACAAAGCCA‑GTGAAAT‑CACTGGCTCGCGTCTTCCGAAG‑ATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAAt                                                                                                                                                                                                                 >  1:19249/1‑251 (MQ=255)
                                                ggAAATGTTCATGCCGTATGGATTGTGGTATCTGGAAACG‑T‑CCTCGCATTTGTTATGCAAAATGCAACAAAGCCA‑GTGAAAT‑CACTGGCTCGCGTCTTCCGAAG‑ATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCg                                                                                                                                                                                                           >  1:391918/1‑250 (MQ=255)
                                                  aaaTGTCCATACCGTATGGATTTTGGTATCTGGAAACG‑T‑CCTCGCTTTTGTTATGAAAAATGCAACAAAGCAA‑GTGAAAT‑CACTGGCTCGCTTCTTCCGAAT‑ATTTTACAAATCGCACAGGGCCAGGCGGCAAGGCCACGCTTAAGATGTTCGGTTGCCATTCTCCCCGGGCAGCAAACTGCCG‑GGAGGA‑CCCGCGCGAAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACg                                                                                                                                                                                                        <  1:175093/251‑1 (MQ=255)
                                                          aTACCGTATGGATTGTGGTATCTGAAACCC‑T‑CCTCGCTTTTTTTAGGAAAAAAGCAACAAGGCCA‑GTAAAAC‑CCCTGTCCCGCTTCCTCCGAAG‑TTTTTTCAAAGCGCACAGGCCCCGGGGGCAGTGCCACGCTTAGGAGGTTGGGTCGCCATCCGGCCCGGGCAGCAACCTTCCG‑CGAGTA‑CCCGCGCCCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGc                                                                                                                                                                                                <  1:213498/251‑1 (MQ=255)
                                                             ccGTATGGATTGTGTTATCTG‑AAACC‑TGCCTCGATTTTTTTAGGCAAAATCAAACAAGGCGG‑GGAAAAT‑CACGGGCTCGGGTCT‑CCGAAGGTTGTTCAAAATCGCAAAAGGCCGGGCGGCAATGCCCCGCTAAGGATGACCGGTGTCCTTTCTGCCCGGGCAGCAACTTTGCG‑CGTGGG‑GCCGCGCGAAACCAGAACCCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCaaa                                                                                                                                                                                             <  1:248163/251‑1 (MQ=255)
                                                                tATGGGTTGTGTTTTCTGAAAACG‑T‑CTTCGCATTTTTTATGAAAAAGGCAACAAAGC‑A‑GGGGAATACCCGGGCTCGCGTCTCCAAAAG‑ATGTTTAAAACCGCAAAAGGCCGGGCGGAAATCCCACGTTTAAGAGGACCGGTGGCAATTCTGCCGGGACAAAAAACTTACG‑CGGGGA‑GCCGCGCGAAAGCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAaga                                                                                                                                                                                           <  2:421418/250‑1 (MQ=255)
                                                                tATGGATTGTGGTATCTGGAAACG‑T‑CCTCGCATTTGTTATGCAAAATGCAACAAAGCCA‑GTGAAAT‑CACTGGCTCGCGTCTTCCGAAG‑ATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAag                                                                                                                                                                                            >  1:350135/1‑249 (MQ=255)
                                                                gaTGGATTGTGGTATCTGGAACCG‑T‑CCTCGCATTTGTTATGAAAAATGCAAAAAAGCCA‑GTGAAAC‑CACTGGCTCGCGTCTCCCGAAG‑ATGTTTCAAATCGCACAAGGCCGGGCGGCAATGCCCCGCTTAAGATGATCGGTTCCCATTCTGCCCGGGCAGAAAACTTCCG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAagag                                                                                                                                                                                          <  1:174494/250‑1 (MQ=255)
                                                                 agggtgtgtgTTATCTGGAAGCG‑C‑CCTGGCATTTGATAAGCAAACTGCAACAAATCCA‑GTGAAAT‑CACTGCCTCGTGTCCTCCGAGT‑ATTTTTAAATTCGCACAGGGCCGGGCGCAAATGCCACCCTTAAGAGGTCCGGTTGCCTTTCCCCCCGGCCAGCAAACTCACG‑TAAGTA‑GCCCCGCAAAATCAGCACCCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAagag                                                                                                                                                                                          <  2:38523/244‑1 (MQ=255)
                                                                 aTGGATTGTGGTATCTGGAAACG‑T‑CCTCGCATTTGTTATGCAAAATGCAACAAAGCCA‑GTGAAAT‑CACTGGCTCGCGTCTTCCGAAG‑ATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAAt                                                                                                                                                                                                                                          >  1:482220/1‑202 (MQ=255)
                                                                          ggTATCTGGAAACG‑T‑CCTCGCATTTGTTATGCAAAATGCAACAAAGCCA‑GTGAAAT‑CACTGGCTCGCGTCTTCCGAAG‑ATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAAt                                                                                                                                                                                >  2:52059/1‑251 (MQ=255)
                                                                           gTATCTGGAAACG‑T‑CCTCGTATTTGTTATGCAAAAGGAAAAAAAGCCA‑GTAAACT‑CTCTGGCTCGCGTCTACCGAAG‑ATGTTTCAAATCGCACAAGGCCAGGCGGAAATCCCCCGCTAAGGATGCGCGGTTGCTGTTCTGCCGGGGCACCAAGCT‑ACGGCTTGGA‑GCCGCGCGCAATCAGGAACCCGCCAATTCCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATa                                                                                                                                                                               <  1:52059/251‑1 (MQ=255)
                                                                              tatctgaaacctcccTCGCATTTGTTATGCAAAATGCAACAAAGCCA‑TTAAAAT‑CACTGGCTCGCGTCTTCAGAAG‑ATGTTCCAAATCCCACAGGGCCAGCCGGCAATCCCACGCTTAAGATGATCGGTTCCCTTTCTGCCCGGGCAGAAAACTTCCG‑CGAGTA‑GCCGCGCGAAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATa                                                                                                                                                                               <  1:97902/237‑1 (MQ=255)
                                                                               cTGGAAACG‑T‑CCTCGCATTTGTTATGCAAAATGCAACAAAGCCA‑GTGAAAT‑CACTGGCTCGCGTCTTCCGAAG‑ATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATagcg                                                                                                                                                                            >  2:421390/1‑250 (MQ=255)
                                                                                  gggaaccgcccccgcttttgtaaggCAAAATGAAAAAAGGCCA‑GTGAAAT‑CGCGGGCCCGCGTCTTCCGAAG‑ATTTTTAAAATCGCACAAGGCCGGGCGCCAAGGCCACGCTTAAGATGTGCGGTTGCCATTCGGCCGGGGAGGAAAACTTACG‑CGAGTA‑GCCGCGCGAAATAAGAACCCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGTTGAAATagcgaa                                                                                                                                                                          <  1:393012/227‑1 (MQ=255)
                                                                                       g‑T‑CCTCGCATTTGTTATGAAAAATGCAACAAAGCCA‑GTGAAAT‑CACTGGCTCGCGTCTTCCGACG‑ATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATa                                                                                                                                                                    >  2:478191/1‑250 (MQ=255)
                                                                                       g‑T‑CCTCGCATTTGTTAGGCCAAATGCAACAAAGCCA‑GTGAAAT‑CACTGGCTCGCGTCTTCCGAAG‑ATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGATTAAGATGATCGGGTGCCATTCTGCCCGGGCAGCAAAATTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATa                                                                                                                                                                    >  2:51610/1‑250 (MQ=255)
                                                                                        cgccctctttttttttATGCAAAATCCAA‑AAAGCCAAGGAAAAT‑CTCTGCCGCGCTTCTTCCGATG‑TTTTTTCAAATCGCACAGGGCAAGGGGACAATCCCCCGCTTAAGATGTTCGTGTCCCATTCTCCCGGGGCAGCAAACTCGCG‑CGAGGA‑GCCGCGCGAAATCAGGACCCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATa                                                                                                                                                                    <  2:545467/237‑1 (MQ=255)
                                                                                                    tGTTATGCAAAATGCAACAAAGCCA‑GTGCAAT‑CACTGGCTCGCGTCTTCCGAAG‑ATGTTTCAAATCGCACAAGGCCAGGCGGCAACGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGGA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAAAAGCGAAGCGATAAGGCATAATccc                                                                                                                                                        >  1:190675/1‑251 (MQ=255)
                                                                                                          gCAAAATGCAACAAAGCCA‑GTGAAAT‑CACTGGCTCGCGTCTTCCGAAG‑ATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGAAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAgc                                                                                                                                                  >  2:150352/1‑251 (MQ=255)
                                                                                                            aaaaTGCAACAAAGCCA‑GTGAAAT‑CACTGGCTCGCGTCTTCCGAAG‑ATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAgcgg                                                                                                                                                >  1:19331/1‑251 (MQ=255)
                                                                                                                 gCAACAAAGCCA‑GTGAAAT‑CACTGGCTCGCGTCTTCCGAAG‑ATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGAc                                                                                                                                            >  1:533618/1‑250 (MQ=255)
                                                                                                                                              cgcgTCTTCCAAAT‑TTTTTTAAAATCGCAAAAGCCCGGGCGGCAATGCCACACTAAAGATGACCGGTTGCCATTCGGCCCGGACAGAAACCTTACG‑CGAGTA‑CCCGCGCGCAATAAGAAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCa                                                                                                                <  1:189643/251‑1 (MQ=255)
                                                                                                                                                 gTCTTCCGAAG‑ATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGACGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGc                                                                                                             >  2:535413/1‑251 (MQ=255)
                                                                                                                                                         aaG‑ATTTTTCAAACCGCAAAAGGCCGGGCGCAAATCCCACGTTAAAGATGATCGGTTCCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTG‑GCCGCCCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCg                                                                                                     <  1:77757/251‑1 (MQ=255)
                                                                                                                                                               gTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGAt                                                                                                >  2:357457/1‑251 (MQ=255)
                                                                                                                                                                tttCAAACCGCACAGGGCCAGGCGGCAGTGCCACGTTTAAGGTGGTCGGTTGCCATTCTGCCCGGGCAGCAACCTGACG‑GGAGCA‑GCCGCGCGAAATCAGGACCCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGAt                                                                                                <  2:483259/250‑1 (MQ=255)
                                                                                                                                                                  ccAAACCGCACAAGGCCGGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATCCGGCCCGGGCAACAAACTGACG‑CGAGTA‑GCCGCGCACAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATccc                                                                                             <  1:421390/250‑1 (MQ=255)
                                                                                                                                                                    aaaTCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATcccc                                                                                            >  1:68058/1‑250 (MQ=255)
                                                                                                                                                                     aaTCGAACAAGGCCAGGCGGCAATGCCACGCTTAAGATGAGCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATACCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATcccaa                                                                                           >  2:390949/1‑248 (MQ=255)
                                                                                                                                                                      aTCGCACAGGGCAGGGCGACAATCCCACTCTTAATATGATCGGTTGCCTTTCCCCCCGGGCAGCAAACTTCCG‑CGAGTA‑CCCGCGCGAAATCAGGACCCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATCCCCAAt                                                                                         <  1:162613/251‑1 (MQ=255)
                                                                                                                                                                      aTCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATCCCCAAt                                                                                         >  2:196213/1‑251 (MQ=255)
                                                                                                                                                                           acaAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATCCCCAATTCCAc                                                                                    >  2:295691/1‑251 (MQ=255)
                                                                                                                                                                             aaGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATCCCCAATTCCACtt                                                                                  >  1:119504/1‑251 (MQ=255)
                                                                                                                                                                                    gggggcAATGCCCCGCTTAAGAGGATCGGTTGCCATTCGGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCGGCGCGAAATCAGGACCCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGAt                                                                                                <  1:408998/227‑1 (MQ=255)
                                                                                                                                                                                    ggcggcAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAACGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGAt                                                                                                >  2:408998/1‑230 (MQ=255)
                                                                                                                                                                                        ccAATGCCACGCTTAAGATGATCGGTTGCCTTTCTGCCCGGGAAGCAAACTTACG‑CGAGTC‑GCCGCGCGAAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGc                                                                       <  1:535413/250‑1 (MQ=255)
                                                                                                                                                                                          aaTGCCACACTTAAGAGGATCGTCTGCCATTCGGCCCGGGCAACAACCTTCCG‑CGAGTG‑GCCCCGCGAAATCAGGACCCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATCCCCAATTCCGCTTTTTGCCCCGGCAc                                                                     <  1:93568/251‑1 (MQ=255)
                                                                                                                                                                                              ccACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGAAAACTTACG‑CGAGTA‑CCCGCGCGAAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAAc                                                                 <  1:254820/251‑1 (MQ=255)
                                                                                                                                                                                               naCGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAAc                                                                 >  2:254820/2‑250 (MQ=255)
                                                                                                                                                                                                  gCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCt                                                             >  1:370353/1‑251 (MQ=255)
                                                                                                                                                                                                         atgatCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATcac                                                       >  2:297869/1‑250 (MQ=255)
                                                                                                                                                                                                         atgatCGGTTGCCATTCTGCCCGGGCAGCAAAATTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATcac                                                       >  1:495546/1‑250 (MQ=255)
                                                                                                                                                                                                          ggatgggTTCCCATTCTGCCCGGGCAGAAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGAACCCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAg                                                     <  1:297869/246‑1 (MQ=255)
                                                                                                                                                                                                                 ttGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACg                                               >  2:553266/1‑250 (MQ=255)
                                                                                                                                                                                                                      aTTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAATAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATg                                         >  2:410740/1‑251 (MQ=255)
                                                                                                                                                                                                                      aTTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATg                                         >  1:170085/1‑251 (MQ=255)
                                                                                                                                                                                                                           gCCCGGGAAGAAAACTTACG‑CGAGTA‑GCCGCGCGAAATCAGAAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGt                                       <  1:357457/248‑1 (MQ=255)
                                                                                                                                                                                                                               gggCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTg                                >  1:441574/1‑251 (MQ=255)
                                                                                                                                                                                                                                  cagcaAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGtt                             >  1:133953/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         cgcgcgCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAac        >  2:121289/1‑251 (MQ=255)
                                                                                                                                                                                                                                                               cAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAACACGCgg  >  2:315223/1‑251 (MQ=255)
                                                                                                                                                                                                                                                               cAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTGTGCGGGCCAAGTTTATTTAACACGCgg  >  2:72044/1‑251 (MQ=255)
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                          
GTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACG‑T‑CCTCGCATTTGTTATGCAAAATGCAACAAAGCCA‑GTGAAAT‑CACTGGCTCGCGTCTTCCGAAG‑ATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACG‑CGAGTA‑GCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGGCGACATATAAACCAGTAATGCGGAGCTAAGCAGGCTGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAACACGCGG  >  NC_000913/2707193‑2707691

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