Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	732,765	0	C	T	23.8%	57.4 / 20.9	43	T1061T (ACC→ACT)	rhsC	rhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base C (19/13);  new base T (3/7);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 1.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.23e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACA  >  NC_000913/732515‑732994                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                      acacGCGGACACAATATGCAGAGCCGCGGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTAGGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACAGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAAcc                                                                                                                                                                                                                                       <  2:215296/251‑1 (MQ=21)  cacGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAAAGGGAACTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACCGAACCCGCATCCAGACGATTTATCAACCGGGGAGCTTCACGCCACTCATCAGGGCTGAAACCg                                                                                                                                                                                                                                      >  1:132690/1‑251 (MQ=255)          cacaATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCggg                                                                                                                                                                                                                              >  1:143622/1‑250 (MQ=21)           acaATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCCTCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCggcg                                                                                                                                                                                                                             >  1:126358/1‑249 (MQ=21)               tatGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCt                                                                                                                                                                                                                         >  1:8173/1‑251 (MQ=21)                 tGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCCTCCAGACGATTTATCCGCCGGGGAGCTTCACGCCACTCATCCGGGTTGGAACCGCCACCGGGGAGCTgg                                                                                                                                                                                                                       >  2:48212/1‑251 (MQ=255)                        ccggtGGTCGAAAGTCGCTATCTTTACGCCCCGCGGGGCCGCGGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAAc                                                                                                                                                                                                                <  2:210796/247‑1 (MQ=255)                              gTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAgcg                                                                                                                                                                                                          >  1:236193/1‑251 (MQ=34)                                cgagaggcGCTATCTTGACGCCCCGGGGGGCCGCAGGGCGGCAAAACGGGTCTGGCGGCGTGAACGGGACCTGACGGGCTGGATGGCGCGGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCGGACCACGATACAGAACGACAAAACCCGCATCCAGACGAGTTACCAGCCGGGGAGCTTCACGCCACTCAGCCGGGTGGAAACCGCCACCGGTGAACTGGCGAAAACGCAGcgcc                                                                                                                                                                                                        <  2:134484/244‑1 (MQ=11)                                cGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGCGGCCAACAGGGCATGGCGCCGTGAACGGGACCCGACGGGCTGGAAGTAGCCGTCAAGGGCACCACAAGTGACCTGGTACGGCTGGGACCGAGACCGCCTGACCAAGATACAGAAAGACAAGACCCGCAACCCGACGATTTATCCACCGGGGAGCTTCACGCCACTCATCAAGGTTGCAACCGCCACCAGTGAACCGGCGAAAACGCCGcgcc                                                                                                                                                                                                        >  1:48262/1‑251 (MQ=2)                                  aaaGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGcgccgc                                                                                                                                                                                                      >  1:303490/1‑251 (MQ=21)                                    aGTCTCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAg                                                                                                                                                                                                    <  1:256746/251‑1 (MQ=18)                                           aTCTTTACGACCCGCGGGGCCGCGGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAACCCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGc                                                                                                                                                                                              <  2:289249/250‑1 (MQ=16)                                                 aCGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATAcc                                                                                                                                                                                       <  1:184151/251‑1 (MQ=33)                                                 aCGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATAc                                                                                                                                                                                        <  1:145355/250‑1 (MQ=33)                                                       cGCTGGGCCGCAGGGGGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTcag                                                                                                                                                                                 <  2:277500/251‑1 (MQ=18)                                                         cTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGGCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTcaga                                                                                                                                                                                >  1:254912/1‑249 (MQ=18)                                                                 cAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCCCGCTACAGAACGACAAAACCCGCATCCCGCCGAATTATCAGCCCGGGAGCTTCACGCCACCCCTCAGGGTTGAAAACGCCCCCCGTGAGCAGGCGAAAACGCACCACCGCAGCCTGGCGGATAcccctgcgccgtccggtg                                                                                                                                                                       >  2:252684/1‑236 (MQ=11)                                                                   gggTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGc                                                                                                                                                                     <  1:153479/251‑1 (MQ=34)                                                                         cAAAACGGGTATGGCGGCGTGACCGTGACCTGACGGGCTGGATGTCTCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGAc                                                                                                                                                               <  1:257220/251‑1 (MQ=255)                                                                           aaaCGGGTATGGCGGCGTGAACGTGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACgg                                                                                                                                                             <  1:19753/251‑1 (MQ=9)                                                                            aaCGGGTATGGCGGCGTGACCGGGACCGGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGGGCCCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGt                                                                                                                                                            <  2:183390/251‑1 (MQ=255)                                                                              cGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTgg                                                                                                                                                          <  2:8173/251‑1 (MQ=14)                                                                                 gTATGGCGGCGTAAACGGGCCCGGACGGGCTGGAGGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAg                                                                                                                                                       <  2:291005/251‑1 (MQ=11)                                                                                       cggcgTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTgtggt                                                                                                                                                 <  1:117291/251‑1 (MQ=32)                                                                                        ggcgTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTgtggtg                                                                                                                                                >  1:126291/1‑251 (MQ=32)                                                                                              aaCGGGACCTGACTGGCTGGCTGTCGCTGTCACGGAAACCGCAAGCGACCTGGTACGGCTGGGGCGGCGACCGCCTGACCACGATACAGCACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCCTCAGGGGTGAAACCGCCACACGTGAGCTGGCGAAAACGCCGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGGCCGTGGCAGGGTGGGGTTcccc                                                                                                                                          >  2:56134/1‑250 (MQ=11)                                                                                                                    gTCGCTGCCACGGAACCCGCAAGTGCCGTGGTACGGCGGGGAGGGCGACCGCCTGACACCGATAAAGAACGACAGACCCCGCATCCAGCCGATTTCTCAGCGGGGGACCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAATACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCg                                                                                                                    <  1:269838/251‑1 (MQ=11)                                                                                                                            cACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTg                                                                                                            >  2:102505/1‑251 (MQ=31)                                                                                                                                  acccGCAAGTGCCCGGGTACGGCGGGGACGGCGACCGCCTGACCACGATACAGAACACCAGACCCCACACCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGGGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGt                                                                                                      <  1:242050/249‑1 (MQ=11)                                                                                                                                    ccgaaaGTGCCCGGGTACGGCGGGGACGGCGCCCGCCGGACCACGATACAGACCGAAAGACCCCGCACCCAGACGATTTATCGGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGa                                                                                                    <  1:96944/247‑1 (MQ=11)                                                                                                                                         aGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGGGaaaacc                                                                                               >  2:189231/1‑248 (MQ=14)                                                                                                                                                                        cTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGccgccg                                                                >  1:195925/1‑251 (MQ=31)                                                                                                                                                                          gacaccGATACAGACCGACAGAACCCGCACCCAGCCGATTAACCAGCCGGGGAGCTCCACGCCACTCATCAGGGTTGAACCCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGccgccg                                                                <  1:173812/244‑1 (MQ=11)                                                                                                                                                                               cGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCgctggctg                                                         >  1:22307/1‑251 (MQ=31)                                                                                                                                                                                             aGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCt                                           >  1:98365/1‑251 (MQ=31)                                                                                                                                                                                                                        cgggggcgcatCACTCCACTCATCGGGTTTGAAACCGCCCCCGGTGCGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATCCCCTTCAGCGGCCCGGTGGCGAAGACGGTGCCATTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGAt                <  1:198325/241‑1 (MQ=11)                                                                                                                                                                                                                        ccGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGcc                    >  1:61448/1‑247 (MQ=34)                                                                                                                                                                                                                          ggggAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGCGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATgg              >  1:182062/1‑251 (MQ=21)                                                                                                                                                                                                                           gggAGATTCACGCCACTCATCAGGGTTGAACCCGCCACCGGTGAGCTGGCGATAACGCAGCGCCGCAGCCTGGAGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTAGCCCGGCTGGAAAGGGAAATCCTGGCTGACAGGGTGAGTGAGGAAAGCAGCCGCTGGCTAGCATCGTGCGGCCTGAAGGTGGCGCAGATGCAAAGCCAGATggg             >  2:150779/1‑250 (MQ=11)                                                                                                                                                                                                                                    cgcccccccatcgGGGTTGAAACGCCCACGGGTGAGCGGGCAAAACCGCAGCGCCGCAGCGTGGGGGATACCCTGCAGCAGCCCGGTGGCAAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGtata    <  1:266864/243‑1 (MQ=11)                                                                                                                                                                                                                                      gCCACTCATCAGGGTTGAACCCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGGGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATaca  <  1:66245/251‑1 (MQ=37)                                                                                                                                                                                                                                      gCCACTCATCAGGGTTGAAACCGCCACCGGTTAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGAATac   >  2:139257/1‑250 (MQ=25)                                                                                                                                                                                                                                               caGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCa                   >  1:96920/1‑225 (MQ=35)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                      ACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACA  >  NC_000913/732515‑732994

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