Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,332,823	0	A	G	34.2%	54.3 / 17.9	38	pseudogene (1580/4605 nt)	yfaS	pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base A (1/24);  new base G (13/0);  total (14/24)
Fisher's exact test for biased strand distribution p-value = 2.59e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTTCGAGACGCAACCATTCATACGTAACAGGAACCTGTTTTGAAGATTCCAGCGCGGCATAACGGAACACAACCGACTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCG‑CTGGCGTC‑GGACA‑CCGTTTCG‑CTG‑CC‑TTCCAGCG‑ACAC‑GGGGAAACGTCCGGCATAACG‑CAAATCG‑TTACCGACCATTGATAATTGCTGAGCG‑CGCAAACTTAACTGCACGCGGGCATTTTTTACCGGCTCGCCATCCGGGTAGAGCAGTTGCAGTTTGCCGCTGACCGCTTCGCCAGTTTTGAACTCTTTTTTGGCGAGAGCTAAACCA  >  NC_000913/2332593‑2333053                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  cTTCGAGACGCACCCATTCATCCGTAACAGGAACCTTTTTTGAAGATTCCAGCGCGGCATAACGGAACACAACCGACTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTaa                                                                                                                                                                                                                                            <  2:472848/238‑1 (MQ=255)  ttCGAGACGCAACCATTCATACGTAACAGGAACCTGTTTTGAAGATTCCAGCGCGGCATAACGGAACACAACCGACTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACg                                                                                                                                                                                                                                          <  1:946403/239‑1 (MQ=255)                                aaCCTGTTTTGAAGATTCCAGCGCGGCATAACGGAACACAACCGACTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCATCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCG‑CTGGCGTC‑GGac                                                                                                                                                                                                                 <  1:469590/232‑1 (MQ=255)                                 aCCTGTTTTGAAGATTCCAGCGCGGCATAACGGAACACAACCGACTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATCGGCGGCGGGGGGGAGTTAAGGCC‑CATGAACG‑CTGGGGTC‑GGTCA‑cccgttc                                                                                                                                                                                                        >  2:655374/1‑234 (MQ=255)                                                      gcgGCATAACGGAACACAACCGACTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCG‑CTGGCGTC‑GGACA‑CCGTTTCG‑CTG‑CC‑TTCCAGCG‑aca                                                                                                                                                                                   <  1:683439/235‑1 (MQ=255)                                                      gcgGAAAAACGGAACACAACCGACGCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCG‑CTGGCGTC‑GGACA‑CCGTTTCG‑CTG‑CC‑TTCCAGCG‑aca                                                                                                                                                                                   <  2:4126/235‑1 (MQ=255)                                                              aCGGAACACAACCGACTCGCCGCTATTACTGTATTGGGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCG‑CTGGCGTC‑GGACA‑CCGTTTCG‑CTG‑CC‑TTCCAGCG‑ACAC‑GGGGAAACGTc                                                                                                                                                                      <  1:194888/239‑1 (MQ=255)                                                                             cTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCG‑CTGGCGTC‑GGACA‑CCGTTTCG‑CTG‑CC‑TTCCAGCG‑ACAC‑GGGGAAACGTCCGGCATAACG‑Caaa                                                                                                                                                       <  1:388402/238‑1 (MQ=255)                                                                                   gcTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCG‑CTGGCGTC‑GGACA‑CCGTTTCG‑CTG‑CC‑TTCCAGCG‑ACAC‑GGGGAAACGTCCGGCATAACG‑CAAATCG‑TTAc                                                                                                                                               <  1:253934/239‑1 (MQ=255)                                                                                       ttACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCG‑CTGGCGTC‑GGACA‑CCGTTTCG‑CTG‑CC‑TTCCAGCG‑ACAC‑GGGGAAACGTCCGGCAt                                                                                                                                                                <  1:259790/220‑1 (MQ=255)                                                                                         aCTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCG‑CTGGCGTC‑GGACA‑CCGTTTCG‑CTG‑CC‑TTCCAGCG‑ACAC‑GGGGAAACGTCCGGCATAACG‑CAAATCG‑TTACCGACCa                                                                                                                                         <  1:45460/239‑1 (MQ=255)                                                                                           tGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCG‑CTGGCGTC‑GGACA‑CCGTTTCG‑CTG‑CC‑TTCCAGCG‑ACAC‑GGGGAAACGTCCGGCATAACG‑CAAATCG‑TTACCGACCAtt                                                                                                                                       <  1:547743/239‑1 (MQ=255)                                                                                               ttGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTTTTACCCCCCCAAACCCC‑CGGGTGTG‑GGACC‑CCGTTTCC‑CTG‑CC‑TCCCCCCG‑CCAC‑GGGGAAACGCCGGCCAAAACC‑AAAACCG‑TT‑CCGACCATTGTTa                                                                                                                                   >  2:360301/1‑239 (MQ=255)                                                                                                  tgCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGGTTTTAACCCCCATGATCG‑CCGGGGGC‑GGAGA‑CCGGTTCG‑CTG‑CCTTTCCTCGG‑CCAC‑GGGGAACGGTCCGCCAAAACC‑CAAATGT‑TTACC‑ACCCATGATaaat                                                                                                                                >  2:600963/1‑237 (MQ=255)                                                                                                         gTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGATTTTACGCCACATGACCC‑CCGGCGTC‑GGACA‑CCCTTTTG‑CTG‑CC‑TTCCAGCG‑CCAC‑GGGGAAACGTCCGGCCTAACG‑CCAATCG‑TTACCCACCATTGATACTTGCTGAg                                                                                                                          >  1:337615/1‑238 (MQ=255)                                                                                                            cTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCG‑CTGGCGTC‑GGACA‑CCGTTTCG‑CTG‑CC‑TTCCAGCG‑ACAC‑GGGGAAACGTCCGGCATAACG‑CAAATCG‑TTACCGACCATTGATAATTGCTGAgcg‑cg                                                                                                                     <  1:360301/239‑1 (MQ=255)                                                                                                                   aGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCG‑CTGGCGTC‑GGACA‑CCGTTTCG‑CTG‑CC‑TTCCAGCG‑ACAC‑GGGGAAACGTCCGGCATAACG‑CAAATCG‑TTACCGACCATTGATAATTGCTGAGCG‑CGCAAAc                                                                                                                <  1:732067/237‑1 (MQ=255)                                                                                                                    gTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGGATTTACGCCCCCATACCG‑CCGGCGGCGGGGCA‑CCGTTT‑G‑CTG‑CC‑CTCCCGCGCACAC‑GGGGAAACGCCCGGC‑TAACGCCAAACTC‑TTCCCGCCCATTTTTAATTG‑TGTGCGCCGCAAACtt                                                                                                              >  2:298519/1‑239 (MQ=255)                                                                                                                       aTGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTTACGCCACATGACCG‑CTGGCGTC‑GGACA‑CCGTTTCG‑CTG‑CC‑TTCCAGCG‑ACAC‑GGGGAAACGTCCGGCATAACG‑CAAATCG‑TTACCGACCATTGATAATTTCTGAGCG‑CGCAAACTTa                                                                                                             >  1:704941/1‑236 (MQ=255)                                                                                                                           gcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCG‑CTGGCGTC‑GGACA‑CCGTTTCG‑CTG‑CC‑TTCCAGCG‑ACAC‑GGGGAAACGTCCGGCATAACG‑CAAATCG‑TTACCGACCATTGATAATTGCTGAGCG‑CGCAAACTTAACTGCAc                                                                                                      <  1:600963/239‑1 (MQ=255)                                                                                                                             cAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGATTTTACGCCACATGACCG‑CTGGCGGC‑GGGCACCCGTTTC‑‑CCG‑CC‑CTCCAGCC‑ACAC‑GGGGGAACGTCCCGCGTAAAG‑CAAAACG‑TTACCCACCCTTTATAATTGCTGCGCG‑CGCGAACCTAAATCCAca                                                                                                     >  2:732067/1‑237 (MQ=255)                                                                                                                                       tCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGATTTAACGCCACATGACCG‑CCGGCGTC‑GGACA‑CCGTTTCG‑CTG‑CC‑TTCCCGCG‑ACAC‑GGGGAAACGGCCGGCATAACG‑CAAATCG‑TTACCGACCCTTGATAATTGCTGAGCG‑CGCAAACTTAACTGCACGCGGGCAttttt                                                                                          >  1:759197/1‑239 (MQ=255)                                                                                                                                            gAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGG‑GGGGGGGGTTTAAGCCCCAAGAACG‑CCGGGGGT‑GGGCACCCGTTTGC‑GGG‑CC‑TTCCA‑CG‑GCACGGGGGAAACGCCGGCCCAAACG‑CAAACGC‑TTCCACCCCATTGTAAATGGCGGAGCG‑CCAAAACTAAACCGCACGCGGGCAtttttt                                                                                         >  2:746197/1‑236 (MQ=255)                                                                                                                                               gATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCG‑CTGGCGTC‑GGACA‑CCGTTTCG‑CTG‑CC‑TTCCAGCG‑ACAC‑GGGGAAACGTCCGGCATAACG‑CAAATCG‑TTAc                                                                                                                                               <  1:639357/179‑1 (MQ=255)                                                                                                                                                  ctctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCG‑CTGGCGTC‑GGACA‑CCGTTTCG‑CTG‑CT‑TTCCAGCG‑ACAC‑GGGGAAACGTCCGGCATAACG‑CAAATCG‑TTACCGACCATTGATAATTACTGAGCG‑CGCa                                                                                                                   <  1:979642/203‑1 (MQ=255)                                                                                                                                                   tctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGATTTTACGCCCCATGACCG‑CTGGCGTC‑GGGCA‑CCGGTTCG‑CTG‑CC‑TTCCCGCG‑GCCA‑GGGGGAAAGCCCGGCCTAACG‑CAAAACG‑TTAACCACCCTTGATAATTGCTGAGCG‑CGCACACTTAAATGAACGCGGGGATTTTTTTCCCGCTCGc                                                                               >  2:389455/1‑238 (MQ=255)                                                                                                                                                         gtgggggtgACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATCGGCGGC‑GGGGGGAGTTAACATCAAATGAACC‑CGGGGGGC‑GGACA‑CCGGTTCGTCTG‑CC‑‑TCTATCG‑ACAC‑GGGGGAAAGTTCGGCACAACC‑TCAAACG‑TTACCGACCACTgt                                                                                                                                     >  2:639357/1‑178 (MQ=255)                                                                                                                                                            gtggtgACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGGATTTACGCCCCCTTACCG‑CTGGGCTC‑GGGCA‑CCGCTTTG‑CTG‑CC‑TTTCCGCG‑ACCC‑GGGGGAACGTCCGGCATAACG‑CAAATCG‑TTTCCCACCATTTATAAATTCTCAGTG‑CGCAAACtttaatc                                                                                                         >  2:979642/1‑198 (MQ=255)                                                                                                                                                             tggtgACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCCCGGTTTATCGGCGGC‑GGGGGGATTTAAAGCCCCAATACCG‑CTGGGGGC‑CGAAA‑CCCGTTCG‑CTG‑CC‑TTCCCGCG‑ACCC‑GGGGAAACGGCCCGCATAACG‑CCAAACG‑TTACCCAGCATTGATAATTTCTGAGCG‑CGCAAAATTAACTGCCCGCGGGCCATTTTTTCCGGCTCGCCATCCcgggtg                                                                    >  2:803371/1‑236 (MQ=255)                                                                                                                                                              ggtgACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCG‑CTGGCGTC‑GGACA‑CCGTTTCG‑CTG‑CC‑TTCCAGCG‑ACAC‑GGGGAAACGTCCGGCATAACG‑CAAATCG‑TTACCGACCATTGATAATTGCTGAGCG‑CGCAAACTTAACTGCACGCGGGCATTTTTTACCGGCTCGCCATCCGGGTag                                                                    <  1:82505/238‑1 (MQ=255)                                                                                                                                                                       aTACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATCGGCGGC‑GGGGGGGATTTACGCCCCCTTACCGCCGGGCGCT‑GGCCC‑CCCTTT‑G‑CTGCCC‑TTCCCCGA‑ACCG‑GGGGAACCGTCC‑GCCTAACG‑CAAAACGTTTACCCTCCATTATAAATTGCTG‑GCG‑CCCCAACCTAACAGCACCCGGGGATTTTTTTCCTGCGCGCCCTCCCGGGTGCgcagtttca                                                          >  2:637976/1‑239 (MQ=255)                                                                                                                                                                        tACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCG‑CTGGCGTC‑GGACA‑CCGTTTCG‑CTG‑CC‑TTCCAGCG‑ACAC‑GGGGAAACGTCCGGCATAACG‑CAAATCG‑TTACCGACCATTGATAATTGCTGAGCG‑CGCAAACTTAACTGCACGCGGGCATTTTTTACCGGCTCGCCATCCGGGTAGAgcagttgca                                                          <  1:238748/238‑1 (MQ=255)                                                                                                                                                                                  cgTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCG‑CTGGCGTC‑GGACA‑CCGTTTCG‑CTG‑CC‑TTCCAGCG‑ACAC‑GGGGAAACGTCCGGCATAACG‑CAAATCG‑TTACCGACCATTGATAATTGCTGAGCG‑CGCAAACTTAACTGCACGCGGGCATTTTTTACCGGCTCGCCATCCGGGTAGAGCAGTTGCAGTTTGCCGCt                                                <  2:704941/238‑1 (MQ=255)                                                                                                                                                                                                    ttAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCG‑CTGGCGTC‑GGACA‑CCGTTTCG‑CTG‑CC‑TTCCAGCG‑ACAC‑GGGGAAACGTCCGGCATAACG‑CAAATCG‑TTACCGACCATTGATAATTGCTGAGCG‑CGCAAACTTAACTGCACGCGGGCATTTTTTACCGGCTCGCCAt                                                                            <  1:777106/192‑1 (MQ=255)                                                                                                                                                                                                     tAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCG‑CTGGCGTC‑GGACA‑CCGTTTCG‑CTG‑CC‑TTCCAGCG‑ACAC‑GGGGAAACGTCCGGCATAACG‑CAAATCG‑TTACCGACCATTGATAATTGCTGAGCG‑CGCAAACTTAACTGCACGCGGGCATTTTTTACCGGCTCGCCATCCGGGTAGAGCAGTTGCAGTTTGCCGCTGACCGCTTCGCCAGtt                                <  1:631537/235‑1 (MQ=255)                                                                                                                                                                                                                   tCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCG‑CTGGCGTC‑GGACA‑CCGTTTCG‑CTG‑CC‑TTCCAGCG‑ACAC‑GGGGAAACGTCCGGCATAACG‑CAAATCG‑TTACCGACCATTGATAATTGCTGAGCG‑CGCAAACTTAACTGCACGCGGGCATTTTTTACCGGCTCGCCATCCGGGTAGAGCAGTTGCAGTTTGCCGCTGACCGCTTCGCCAGTTTTGAACTCTTTTTTGGCg              <  1:746197/239‑1 (MQ=255)                                                                                                                                                                                                                   tCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCG‑CTGGCGTC‑GGACA‑CCGTTTCG‑CTG‑CC‑TTCCAGCG‑ACAC‑GGGGAAACGTCCGGCATAACG‑CAAATCG‑TTACCGACCATTGATAATTGCTGAGCG‑CGCAAACTTAACTGCACGCGGGCATTTTTTACCGGCTCGCCATCCGGGTAGAGCAGTTGCAGTTTGCCGCTGACCGCTTCGCCAGTTTTGAACTCTTTTTTGGCg              <  2:950268/239‑1 (MQ=255)                                                                                                                                                                                                                                ggc‑ggGGAGATTTAACGCCACATGACCG‑CTGGCGTC‑GGACA‑CCGTTTCG‑CTG‑CC‑TTCCAGCG‑ACAC‑GGGGAAACGTCCGGCATAACG‑CAAATCG‑TTACCGACCATTGATAATTGCTGAGCG‑CGCAAACTTAACTGCACGCGGGCATTTTTTACCGGCTCGCCATCCGGGTAGAGCAGTTGCAGTTTGCCGCTGACCGCTTCGCCAGTTTTGAACTCTTTTTTGGCGAGAGCTAAACCa  <  1:660587/238‑1 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  CTTCGAGACGCAACCATTCATACGTAACAGGAACCTGTTTTGAAGATTCCAGCGCGGCATAACGGAACACAACCGACTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCG‑CTGGCGTC‑GGACA‑CCGTTTCG‑CTG‑CC‑TTCCAGCG‑ACAC‑GGGGAAACGTCCGGCATAACG‑CAAATCG‑TTACCGACCATTGATAATTGCTGAGCG‑CGCAAACTTAACTGCACGCGGGCATTTTTTACCGGCTCGCCATCCGGGTAGAGCAGTTGCAGTTTGCCGCTGACCGCTTCGCCAGTTTTGAACTCTTTTTTGGCGAGAGCTAAACCA  >  NC_000913/2332593‑2333053

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