Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	730,748	0	T	C	32.4%	60.9 / 10.0	37	I389T (ATC→ACC)	rhsC	rhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base T (19/6);  new base C (0/12);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 1.00e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGG‑CCGG‑ATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGT‑TACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGG‑CAGGCTTAAGTTACAC‑GTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCAC  >  NC_000913/730506‑730995                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        cgcgccgcTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACCCAGGTGCGCAGCTTTACTTACGATGATAAAGACCGGGG‑CCGG‑ATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGT‑TACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGG‑CGGGCTTAAGTTACAC‑GTATCAGTATGAGAAGGACCGCATCACCATCACCGACa                                                                                                                                                                                                                                                 <  1:307061/251‑1 (MQ=255)    gccgcTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGTCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGG‑CCGG‑ATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGT‑GACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGG‑CAGGCTTAAGTTACTC‑GTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGcc                                                                                                                                                                                                                                              <  2:228969/251‑1 (MQ=255)    gccgcTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGG‑CCGG‑ATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGT‑TACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGG‑CAGGCTTAAGTTACAC‑GTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGcc                                                                                                                                                                                                                                              <  2:168225/251‑1 (MQ=255)                            cccccccGGCGAAGTGGCGGCGGTGTAGGCCCGCAGCAAAACACAGGTGCGCATCTTTACTTACGATGA‑AAAAACCGGGGGCCGG‑ATGGTGGCGCCCCGCCCCACGGGCCGCCCGAAAA‑TCTGTTCCCGT‑TACGACAGCGACGGGCGGGTGAAAGAACAGCAAAACCCGG‑CAGGCTTAAGTTACAC‑GTAGCAGTATGAGAAGGACCGCATCCCCACCACCGACAGCCTGAACCGCCGTGAAGTCCTGcaca                                                                                                                                                                                                                      <  1:212163/250‑1 (MQ=255)                                   ggcgAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGG‑CCGG‑ATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGT‑TACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGG‑CAGGCTTAAGTTACAC‑GTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAggg                                                                                                                                                                                                               >  2:180143/1‑251 (MQ=255)                                        aggggcgggGGTGTATGCCCGCACAAATCCACGGGTGCGCAGTTTTCTTTACGATGATAAACACGGGGG‑CCGG‑AGGTTGGCGCCCCGTCACAGGGGCCGACCGAAAA‑TCTGTTACCGT‑TCCGACAGCGAGGGGGGGTTGACAAAAAAGCAAACCCCGG‑CGGGTTTAAGTAACAC‑GTACCAGTATGAGAGGGCCCGCACCCCCACCACCGCCAGCCTAAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAg                                                                                                                                                                                                          <  1:98369/248‑1 (MQ=255)                                                  gtgtATGACCGCAGCAATTCACAGGTGCGCAGCTTTCTTTGCGATGAAAAACCCCGGGG‑CCGG‑AGGTTGGCGCCCCCCCACACGGCCCG‑CCGGAAAATCTTTTCCCTT‑TACGCCGGCGACGGGGGGGTGAAAGACCACCAAAACCCGG‑TAGCCTTAAGTTCCCC‑GTATCAGTATAAGAAGGACCCAATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGaa                                                                                                                                                                                                <  1:62218/251‑1 (MQ=255)                                                    gtATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGG‑CCGG‑ATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGT‑TACGCCAGAGACGGGCGGGTGACAGAACAGCTAAACCCGG‑CAGGCTTAAGTTACAC‑GTATCAGGATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAgg                                                                                                                                                                                              >  2:216338/1‑250 (MQ=255)                                                       tGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGG‑CCGG‑ATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGT‑TACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGG‑CAGGCTTAAGTTACAC‑GTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCggg                                                                                                                                                                                           >  1:51562/1‑251 (MQ=255)                                                            acagcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGG‑CCGG‑ATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGT‑TACGACAGCGCCGGGCGGGGGACAGAACAGCTAAACCCGG‑CAGGCTTAAGTTCAAC‑GTATCAGTATGAGAAGGACCGCATCACCATCACCGACAACCTTAAACGCCGTGAAGTCCTGCACACGCAGGGTGCAGGAGGGATGACGCGGgtggtg                                                                                                                                                                                      >  2:187876/2‑251 (MQ=17)                                                                caATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGG‑CCGG‑ATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGT‑TACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGG‑CAGGCTTAAGTTACAC‑GTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGaaa                                                                                                                                                                                   >  2:257897/1‑250 (MQ=255)                                                                     cacgggTGCGCACTTTTTCTTACGATAAAAAATACGGGGG‑CGGG‑AGGGTGGCGCCCCGTCCCAGGGGCCGACGGAAAA‑TTGGTTACCTT‑TACGACAGCGAGGGGGGGGGGACAGAACAGCTAACCCCGG‑AAGGCTAAAGTTACCC‑GTATATGTTTGAGAAGGCCCGCATCCCCACCACCGCCAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGaa                                                                                                                                                                             <  1:180143/247‑1 (MQ=255)                                                                       caGGTGCGCAGCTTTACTTACGATGATAAATACCGGGG‑CCGG‑ATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGT‑TACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGG‑CAGGCTTAAGTTACAC‑GTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAac                                                                                                                                                                            >  1:166472/1‑250 (MQ=255)                                                                        gggTGGGCAGATTTACTTACGAGGAAAAAGAC‑GGGG‑CCGGGAGGGTGGCGCACTTCCACAGGGGCCGAGCGGAAA‑TC‑GTTTCCCTGTACGGCAGCGACGGGGGGTTAACAAAACCGCTAAGCCCGA‑CAGGCTTAATTTA‑ACATTATCAGTATGAGAGGGCCCGCATCACCACCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAacac                                                                                                                                                                          <  1:57907/250‑1 (MQ=9)                                                                         ggTGCGCAGCTTTACTTACGATGATAAATACCGGGG‑CCGG‑ATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGT‑TACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGG‑CAGGCTTAAGTTACAC‑GTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACAcg                                                                                                                                                                         >  1:228902/1‑251 (MQ=255)                                                                           ggggcatcTTTCCTTACGAAGAAAAATCCGGGGG‑CCGA‑ATGGTGGCGCCCCGCCACCCGGGCCGACGGAAAA‑TCTGTTCCCGT‑TACGCCAGCGACGGGCGGGTGACAGAACAGCAAACCCCGG‑CGGGCTAAATTTACCC‑GTATCAGTATGAGAAGGCCCGCATCCCCATCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACAcgcg                                                                                                                                                                       <  1:219128/244‑1 (MQ=14)                                                                                              tgatAAATACCGGGG‑CCGG‑ATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGT‑TACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGG‑CAGGCTTAAGTTACAC‑GTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGc                                                                                                                                                       >  2:321706/1‑248 (MQ=255)                                                                                                  aaatccggggg‑CCGG‑GTGGGGGCGCACCGTCACACGGCCCGACCGGAAA‑TTTGTCACCTT‑TACGACGGCGGGGGGGGGGTAACAAAACAGC‑AAACCCGGGCAGGCTAAAGTTACAC‑GTATCAGAAGGGGAAGGACCGAATCACCACCACCGACACCCGGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCa                                                                                                                                                <  1:124026/244‑1 (MQ=9)                                                                                                   aatacagagg‑CCAG‑ATGGTGGCGAACCGTCACACGGGCCGACCGGCAA‑TCTGTTACCGT‑TACGATAGCGACGGGCGGGTGACAGAACAGCTAAACCCGG‑CAGGCTTACGTTACAC‑GTATCAGTATGAGAAGGACCGACTCACCATCACCGACAGCCTGAACCGCCGTGAAGTCATGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGGTAGAAAGGACCACGCGGACGGCAGCGGCACGCagaggc                                                                                                                                                 >  2:187439/9‑247 (MQ=14)                                                                                                                                         gccggccgAAAA‑TCTGTTCCCGT‑TACGACAGCGACGGGCGGGTGAAAAAACACAAAAACCCGG‑CAGGCTAAAGTAACCC‑GTAGCAGTATGAGAAGGACCGCACCCCCACCCCCGACAGCCTAAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGAt                                                                                                           <  1:33417/251‑1 (MQ=14)                                                                                                                                                            aCCGT‑TACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGG‑CAGGCTTAAGTTACAC‑GTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGacaac                                                                                          >  1:293699/1‑250 (MQ=39)                                                                                                                                                                  tACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGG‑CAGGCTTAAGTTACAC‑GTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAAt                                                                                     >  1:225483/1‑250 (MQ=255)                                                                                                                                                                           tcggatcgggggaacaaaacagaaaaCCCGGGCAGGCATAGTTTACAC‑GTAACAGTATGGAAAGGCCCGCAACCCAACCCCCGAAAGCCGAAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCg                                                                             <  2:103318/227‑1 (MQ=14)                                                                                                                                                                                  gggTGACAGAACAGCTAAACCCGG‑CAGGCTTAAGTTACAC‑GTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCg                                                                             >  1:103318/1‑242 (MQ=255)                                                                                                                                                                                       gaaaaacaagcaaaccccgacgggCTTAAGTCACCA‑GTATCAGTATGAGAGGGCCCCAACCACCACCCCCGACAGCCGGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACgg                                                                <  2:121907/229‑1 (MQ=14)                                                                                                                                                                                         aGAACAGCTAAACCCGG‑CAGGCTTAAGTTACAC‑GTATAAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCTCAGGCAGGACAACAGAATAAAGTCCGGATGTGGTGACGGGcc                                                             >  1:221106/1‑251 (MQ=255)                                                                                                                                                                                            aacatcaaaccccggcatgttcaatttACAC‑GTATCAGTAGAAAAAGCCCCGCACCACCACCCCCGACAGCTTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGTGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCtc                                                           <  1:326622/226‑1 (MQ=14)                                                                                                                                                                                              aGCTAAACCCGG‑CAGGTTTAATTTACAC‑GTATCAGTAGGAGAAGGACCCCACCACCATCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCtcatc                                                        <  1:69279/251‑1 (MQ=255)                                                                                                                                                                                              aGCTAAACCCGG‑CAGGCTTAAGTTACAC‑GTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCtcatc                                                        >  1:162069/1‑251 (MQ=255)                                                                                                                                                                                               gCTAAACCCGG‑CAGCCTTAATTGACAC‑GTATCAGTATGAGAAGGCCCGCATCCCCACCCCGGCCAGCCGGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCtcatca                                                       <  1:39421/251‑1 (MQ=14)                                                                                                                                                                                                                ttAATTTACAC‑GTATCAGTATGAGAAGGACCGCATCCCCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCCCGcc                                       >  2:186838/1‑251 (MQ=255)                                                                                                                                                                                                                    gttccac‑GTATCAGTAGGAGAAGGACCGCATCACCACCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGAt                                   <  1:257897/247‑1 (MQ=255)                                                                                                                                                                                                                      tACAC‑GTATCAGTATGAGCAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGACAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATgg                                 >  2:52791/1‑251 (MQ=255)                                                                                                                                                                                                                                   tATGAGAAGGACCGCATCACCATCACCGACAGCCTGATCCGCCGTGAAGTCCTCCACACGCAGGGTGAAGGCGCGCTGCAGCCCGTGATGAAAAAGGAACACGCGGCCGGCAGCGTCACGCCGAGTCCGCTTGACGCGGTGGGCAGGCTCAGGGCACCGACGGATGCCGCCGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCGTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGc                     >  1:44250/1‑251 (MQ=255)                                                                                                                                                                                                                                    aggagaAGGCCCGCATCACCACCACGGACGGCTTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATcaccaagaaggagggaagggaagcgacg                    <  1:241013/249‑24 (MQ=14)                                                                                                                                                                                                                                      gagaAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGtt                  >  1:153363/1‑251 (MQ=255)                                                                                                                                                                                                                                          aGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTAc              >  2:29355/1‑251 (MQ=255)                                                                                                                                                                                                                                                 catcaccatcacCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAacc       >  1:158724/1‑251 (MQ=255)                                                                                                                                                                                                                                                      ccatcaccGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACcaccac  >  1:270886/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        CGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGG‑CCGG‑ATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGT‑TACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGG‑CAGGCTTAAGTTACAC‑GTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCAC  >  NC_000913/730506‑730995

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