Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,346,653	0	C	A	28.9%	58.5 / 14.1	45	L30F (TTG→TTT)	pdeR	cyclic di‑GMP phosphodiesterase PdeR
Reads supporting (aligned to +/- strand):  ref base C (24/8);  new base A (2/11);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 4.97e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAA  >  NC_000913/1346405‑1346902                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                            gCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCa                                                                                                                                                                                                                                                            >  1:257312/1‑250 (MQ=255)     aGGAAAGCCCGGGTGATAACTCACGAGAAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACa                                                                                                                                                                                                                                                       <  1:192618/251‑1 (MQ=255)               gtgtGATAAATCACGAGCAACCGCCGGTGTGTCATGACATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCATACCAGA‑TGCACAGAAACATC‑GCTTTCCTCGACGGTCC‑GACTCAT‑CATCAGGCTGGAGGTGATGGAC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGAACACTGCAATGA‑ACGATCAGGGTCTGTGGTTTCGGTGG‑TAGAAAA‑GCAAAAAACATCCCTGCttt                                                                                                                                                                                                                                             >  2:694046/1‑251 (MQ=255)                tgtgATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTc                                                                                                                                                                                                                                            >  1:499177/1‑251 (MQ=255)                 gtgAAAAACCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTc                                                                                                                                                                                                                                            <  1:545813/250‑1 (MQ=255)                    aTAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTc                                                                                                                                                                                                                                        <  1:83887/251‑1 (MQ=255)                            gaaGCACCCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCa                                                                                                                                                                                                                              <  1:81526/249‑1 (MQ=255)                             ggcccaCCGCCGGTGTGTCATGCCATTCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGTTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAg                                                                                                                                                                                                                             <  2:493242/246‑1 (MQ=255)                             gAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCAATTCCCGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑Gcc                                                                                                                                                                                                                               >  2:264331/1‑249 (MQ=255)                                        ggTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGAAATAT‑AAGATCAGGTTCTGTGGTTTCGGTGG‑TAACCAA‑ACGAAAAACATCGATTCTTTCAGTC‑AATC‑GCCCGTAGGGAGtgtg                                                                                                                                                                                                                  >  2:112197/1‑251 (MQ=255)                                             gtCAGGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGA‑TGCAAAAAAAGAGC‑GCTTTCTCCGAGGGTCA‑AACT‑ATACCTAGGGCGGGAGGTGATGACC‑GTCTTTCCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACGGAAATTA‑ACGATAAGGTTTTGGGTTTTGGGGGG‑TAGAAAA‑GCAAAAAACATCGCTTTTTCCCGCC‑AACC‑GCCAGTATGGATTGTGCGAt                                                                                                                                                                                                              <  1:417014/250‑1 (MQ=255)                                                             gcgTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACa                                                                                                                                                                                              >  2:751109/1‑250 (MQ=255)                                                                  cccacccaaTTCCGGTTTTTTGTTTTTTC‑TCCTACCAAAAGGCACAGAAA‑ATCCGCTTTCATCGACG‑TCAAGACTCAT‑CATCGGGTTGGGGGTTA‑GACCCTTCATTTCCCTG‑AGGCGAGCAGCTTTTTCGG‑CAGAAAACTGAAATTG‑ACGTTCTGTTTTTGGGTTTTGGGTGG‑TAAAAAA‑CCGAAAAACACCGCTGTTTCCCATC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGtt                                                                                                                                                                                         <  1:722993/242‑1 (MQ=255)                                                                             cccGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCt                                                                                                                                                                              >  1:575805/1‑250 (MQ=255)                                                                                   tATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGtt                                                                                                                                                                       >  2:182708/1‑251 (MQ=255)                                                                                    aTTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTACTCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACA‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGAATGGATTGTGCGATCCGAGCACGTTGTACAACGTTGTGGACTACCTAACGGttt                                                                                                                                                                      >  2:326631/1‑251 (MQ=255)                                                                                               gTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGAACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCat                                                                                                                                                           >  1:306902/1‑251 (MQ=255)                                                                                                       aacatgaaaaaaaaGGTC‑GCTTTCGTCGACGTGCA‑GACTCAT‑CACCGGGCGTGAGTGGATGCCC‑TTC‑TTTCCCTGTATGGGAGCGGCCT‑TTCGGGCAGATAAC‑GCAATTTTAGGATGAGGTTGTGTGGTTGGGGGGG‑TAGAAAG‑GCGAAAAACCTCGCTGCTTTCCTTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATAtcctcc                                                                                                                                                    <  2:269068/240‑1 (MQ=255)                                                                                                         a‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGCGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGt                                                                                                                                                 >  1:140756/1‑251 (MQ=255)                                                                                                                           ttttCATCGACGTTCA‑GCCTAAC‑CATCGGGGGGGGGGTGAGGCCC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTGCGG‑AAGATAACGGAAATGT‑CGGAACGGGTTTTGGGGTTT‑GGTGGGAAAAAAA‑GCGAAAAACACCGCTGCTTCCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTcc                                                                                                                                   <  1:590688/248‑1 (MQ=255)                                                                                                                             ttCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCgg                                                                                                                               >  1:660372/1‑251 (MQ=255)                                                                                                                                   gACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaact                                                                                                                         >  2:167234/1‑249 (MQ=255)                                                                                                                                   gACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGAGGTtaaat                                                                                                                         >  1:509205/1‑251 (MQ=255)                                                                                                                                      ccggccagcctcatcagcgggggggAGGTGATGACC‑GTCATTCCCGTG‑ATGCGAGCAGCCTTTTCGG‑CAGATACCTGAAATGT‑ACGTTCAGGTTCTGTGTTTTGGGTGG‑AAAAAAA‑GCGAAAACCACCCCTGTTTCCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaata                                                                                                                        <  1:339316/228‑1 (MQ=255)                                                                                                                                       gTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaataaa                                                                                                                      >  1:465890/1‑250 (MQ=255)                                                                                                                                                 at‑catcaGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaacaac                                                                                                             >  2:615991/1‑250 (MQ=255)                                                                                                                                                        ggggggggggTGATGACC‑GTCTTTTCCCGG‑AGGCGAGCTGCCTTTTCGG‑AGGAAAACTGCAATGT‑ACGCTCGGGTTCTGTGTTTTGGGTGG‑TAAAAAA‑GCGCAAAACTCCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACtt                                                                                                      <  1:615991/243‑1 (MQ=255)                                                                                                                                                                ggTGATGACC‑GTCTTTTCCTTG‑AGGCGACCAGCCTGTTCGG‑AAGATAACTGAATTGT‑CGGATCGGGTTTTGGGGTTTCGGTGG‑TAAAAAA‑GCGCAAAACACCGCTGCTTCCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGaaa                                                                                              <  1:750916/251‑1 (MQ=255)                                                                                                                                                                               cccgatcggaggcgggcagcCTGTTCGG‑CAGATAACTGAAATGT‑ACGATCAGGTTCTGTGGTTTGGGTGG‑TAAAAAA‑GCGAAAAACATCGCTGCTTCCCGTA‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTa                                                                                 <  1:748340/236‑1 (MQ=255)                                                                                                                                                                                  cccTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGCAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCtcaa                                                                               >  2:272560/1‑246 (MQ=255)                                                                                                                                                                                      g‑ATGCGAGCAGCTTTTCCGG‑CAAATAACGGCAATGT‑CCGATCAGGTTCTGTGTTTTCGGTGG‑TAAAAAA‑GCGCAAACCATCGTTTTTTCCCGTC‑ATCC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCa                                                                           <  1:514655/249‑1 (MQ=255)                                                                                                                                                                                      g‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTa                                                                                 >  2:748340/1‑243 (MQ=255)                                                                                                                                                                                        aTGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTg                                                                        >  1:574844/1‑251 (MQ=255)                                                                                                                                                                                                agcCTGTTGGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGGGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCACTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAAATCATAGtttttttt                                                                >  2:513196/1‑251 (MQ=255)                                                                                                                                                                                                 gcCTGTTGGG‑AAGAAACCTGAAATGT‑ACGATCGGGTTCTGTGGTTTGGGGGG‑AAAAAAA‑GCGAAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTta                                                               <  1:182708/251‑1 (MQ=255)                                                                                                                                                                                                     tgtggggcagaaacgttcAATGT‑ACGATCGGTTTTTGTGGTGCGGGTGG‑TAAAAAA‑GCGAAAAACATCGCTTCTTTCCG‑C‑AACCCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTtata                                                             <  1:751109/233‑1 (MQ=255)                                                                                                                                                                                                                aaacctgaattatccgttaaGGTTCTGTGGTTTGGGTGG‑AAAAAAA‑GCGAAAAACAGCGCTGCTTCCCGAC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGATTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGt                                                  <  1:264331/232‑1 (MQ=255)                                                                                                                                                                                                                 gaacatcaatttccaaTCAGTTTTTGTGGTTTCGGTGT‑AAGAAAA‑GAGAAAAACTTCGCT‑CCTGCCCTCTACCC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGt                                                  <  1:530907/235‑1 (MQ=255)                                                                                                                                                                                                                      ccaatttacgctcggtttCTGTGTTTTCGG‑GG‑TAAAAAAAGCAAAAAACATCGCTG‑TTTCCCTT‑AACCCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCa                                            <  1:272560/235‑1 (MQ=255)                                                                                                                                                                                                                                   aGGTTCTGTGGTTTCGGTGG‑TAGAAAC‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCCGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGGTCAACGACTTGTTTTTTTTATTTATTTTTTTGTAATCAATTTTTGGtaaaa                                >  2:483421/1‑250 (MQ=255)                                                                                                                                                                                                                                          gtgGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATAGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAAAGCATTGTTTTTTTTATATATTTTTTTGTAATCCAGTTTTGGAAAAACGCCAGtt                        >  1:598535/1‑251 (MQ=255)                                                                                                                                                                                                                                              tttCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTTTATTTTTATTTTTTATCCAGTTTTGGAAAAATGCCCGTTTTc                     >  2:463732/1‑250 (MQ=255)                                                                                                                                                                                                                                                tCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTTTATATTTATTTGTAATCCAGTTTTTGaaaaaagcacgtttttta                   >  2:247201/1‑237 (MQ=255)                                                                                                                                                                                                                                                           aaaa‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGCGAAAGTTGTCCAACGTTGTGGACTCCCTTACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGACTCCGCGGTTAAATAAAAGGAACAACATGACt                                                                                                       >  1:736961/1‑156 (MQ=255)                                                                                                                                                                                                                                                           aaaa‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAAcggc                                                                                      >  2:194386/1‑173 (MQ=255)                                                                                                                                                                                                                                                                   cAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTTTTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAATCAGTTaaaa  >  1:647266/1‑250 (MQ=255)                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                            GCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAA  >  NC_000913/1346405‑1346902

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