Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	12,384	T→C	G74G (GGT→GGC)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	12,384	0	T	C	98.0%	158.3 / NA	50	G74G (GGT→GGC)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (23/26);  minor base A (1/0);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGAC  >  NC_000913/12141‑12630                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       aaTATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCt                                                                                                                                                                                                                                                 <  1:45101/251‑1 (MQ=255)  atatatAGTGGAGACGTATAGATGGGTAAAATAATTGATATCGACCTGGGTACAACCAACTCTTGTGTAACGATTATGGATGGCACCACTCCTAGTGTGCTGGAGAACGCCGAAGGAGATCGCACCACGCCTTCTATCATTGCCTATACCCAGAATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGGGACGAACCCGCAAAACACACTGTTTGCGATTAAACACCTGATTAGACGCCGCtt                                                                                                                                                                                                                                                >  2:363655/1‑251 (MQ=255)       gaGTGGAGACGGTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACGCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAgg                                                                                                                                                                                                                                           <  2:341268/250‑1 (MQ=255)                   ttAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTa                                                                                                                                                                                                                                <  1:509727/250‑1 (MQ=255)                           gTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCAGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTAAAGCGTGAt                                                                                                                                                                                                                       >  2:331672/1‑251 (MQ=255)                              aaaTAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGGACAGCGTGATGtt                                                                                                                                                                                                                    <  1:169645/251‑1 (MQ=255)                               aataatTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGCAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTg                                                                                                                                                                                                                         >  2:271817/1‑245 (MQ=255)                               aataatTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTg                                                                                                                                                                                                                         <  1:271817/245‑1 (MQ=255)                                      ggTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCcatca                                                                                                                                                                                                             <  2:145307/250‑1 (MQ=255)                                                        acCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGTACTAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAattat                                                                                                                                                                                             <  2:96967/248‑1 (MQ=255)                                                           aaCTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATGGCACCACGCCTTCTAGCATTGCCTATTCCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGGGATGTTTCCATCATGCCGTTCAAAATTCTtgctgc                                                                                                                                                                                       <  2:277247/251‑1 (MQ=255)                                                           aaCTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATtgctgc                                                                                                                                                                                       <  1:319342/251‑1 (MQ=255)                                                             ctctTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGctgctg                                                                                                                                                                                     >  1:627541/1‑251 (MQ=255)                                                                ttGTGTAGCGAATATGGATGGAACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACAACGCCTTCTATAATTGCCTATACCCAGGATGGAGAACCTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGACCCCGCAAAACACCCTGTTTGCGCTTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGAGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGAt                                                                                                                                                                                   >  2:618043/1‑250 (MQ=255)                                                                  gtgtAGCGATTATTGATGGCACCACTCCTCGTGTGCTGTAGAACGCCGAAGGCGATCGCATCACGCCTTCTATCATTGCCTAGACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGGCAGGCAGTTACGAACCCGCAAAACACCCTGTTTGCGATTAAACGACTGATTGGCCGCCGCTGCCAGGACGAAGAATTACAGCATGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATaa                                                                                                                                                                                 <  2:343239/250‑1 (MQ=255)                                                                    gtAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACg                                                                                                                                                                               <  1:577177/250‑1 (MQ=255)                                                                      aGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGTCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCg                                                                                                                                                                            <  2:522983/251‑1 (MQ=255)                                                                          aTTATGGATGGCACCCCTCCTCGTGTGCTGGAGAACGCCGAAAGAGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCCTGCCGTTCCAAATTATTGCTGCTGATAACGGCTACGc                                                                                                                                                                        >  2:548879/1‑251 (MQ=255)                                                                             atggatggCACCACTCCGCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATg                                                                                                                                                                     <  2:307155/251‑1 (MQ=255)                                                                                gatggCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGt                                                                                                                                                                  <  1:433122/251‑1 (MQ=255)                                                                                 atggCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGAGTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCAGGCCGTGCAAAAGTATTGCTGCTGAGAACGGCGACGGATGGGTc                                                                                                                                                                 <  2:264914/251‑1 (MQ=255)                                                                                   ggCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGACCCCGCAAAACACCCTGTTTGCGATTAACCGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCgt                                                                                                                                                               >  2:266906/1‑250 (MQ=255)                                                                                   ggCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCGTCCAGGACTAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCg                                                                                                                                                                <  1:610700/250‑1 (MQ=255)                                                                                      accacTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGt                                                                                                                                                            >  1:146957/1‑251 (MQ=255)                                                                                               cgggTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGc                                                                                                                                                     <  2:534782/246‑1 (MQ=255)                                                                                                gtgTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCATATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGACGTTAAAGGca                                                                                                                                                    >  2:612287/3‑248 (MQ=255)                                                                                                                  aaGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACcgccgc                                                                                                                                  >  1:162937/1‑249 (MQ=255)                                                                                                                        aTCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATctct                                                                                                                          >  1:336696/1‑251 (MQ=255)                                                                                                                         tCGCACCACGCCTACTATCATTGCCTATACCCCGGCTGGTGAAACTCTGGTTGGTCAGCCGGCTAAAAGTCAGGCAGTGAAGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGCCGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCCAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCAAGAAAATGGCACCGCCGCAGATCTctg                                                                                                                         >  2:36059/1‑251 (MQ=255)                                                                                                                          cgcccaacGCCTTCTACTATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGGCAGCCGGCTAAACGTCAGGCAGTGACGACCCCGCAAACCACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGGGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATACCGGCGCCGCAGGGTTCGAAGTAAAAGGCCAGAAAAGGTCTCCGCCGCAGATCTctgc                                                                                                                        <  2:402940/245‑1 (MQ=255)                                                                                                                                               gCCTATACCCCGGGTGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCGGGACGACGAAGTACAGCGTGATGTTTCCATCATGCCGTTCACAATTATTGCTGCTGATAACGGCGACGCATGGGTCGGAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATga                                                                                                    >  2:84621/1‑250 (MQ=255)                                                                                                                                                       ccAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGAGACGCACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAACTTATTGCTGCTGATACCGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCg                                                                                           >  1:312717/1‑251 (MQ=255)                                                                                                                                                          ggATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGCAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGACGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTg                                                                                        >  2:231090/1‑251 (MQ=255)                                                                                                                                                                tGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGAtt                                                                                  >  1:501997/1‑251 (MQ=255)                                                                                                                                                                  aaaCTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACCCCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTAc                                                                                <  1:525592/251‑1 (MQ=255)                                                                                                                                                                           ttGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTACTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGa                                                                        >  1:417011/1‑250 (MQ=255)                                                                                                                                                                            tGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAAc                                                                      >  1:263312/1‑251 (MQ=255)                                                                                                                                                                             ggccAGCCGGCTAAACGTCAGGATGTGACGAACCCCAAAAACCCCCTGTTTGCGATAAAACTCCTGATGGGCCGCCGCTTCCAGGACAAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGCCGCATGGGTCGAAGTTAAAGGCCAAAAAATGGCACCGCCGCAGAGCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAAc                                                                      <  2:395852/247‑1 (MQ=255)                                                                                                                                                                             ggTCATCCGGCTAACCGTAAGGCAGTGACGACCCCGCAAAACACCCGGTTGGCGATTAACCGCCTGATTGTCCGCCGCTCCCAGGACGAAAAAGTACAGCGTGATGTTTCAATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAAGGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAAcc                                                                     <  1:616001/251‑1 (MQ=255)                                                                                                                                                                             ggTCAGCCGGCTAAACGTCAGACAGTGACGGAACCGCAAAACCCCATGTTTGCGCTTAAACGACTGATTGGCCGCCGCTTCCAGGACGAAGAAGAACCGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGACAAAAATGAAGAAAAACGCTGAAGATTACATGGGTGAAcc                                                                     >  2:591604/1‑251 (MQ=255)                                                                                                                                                                                cAGCCGGCTAAACGCCAGGCAGTGCCGAACCCGCAAAACCCCCTGTTTGCGATTAAACGCCGGATTGGCCGCCGTTTCCAGGACGAAGAAGTACAGCGGGATGTTTCCATCATGCCGTTCAAAATTATTGCGGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCAAGAAAAGGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGt                                                                  <  2:501997/251‑1 (MQ=255)                                                                                                                                                                                            cGTCAGGAAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGACTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGAACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGt                                                      >  2:393122/1‑251 (MQ=255)                                                                                                                                                                                                          gAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAACATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCgg                                        >  1:281324/1‑251 (MQ=255)                                                                                                                                                                                                            ccccGCAAACCCCCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTGCCAGGACAAAGAAGTACAGCGTGATGTTTCCATCATGCCTTTCAAAATTATTGCGGCTGATACCGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCa                                      <  1:581148/250‑1 (MQ=255)                                                                                                                                                                                                              ccGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCGAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCAt                                     >  1:487127/1‑250 (MQ=255)                                                                                                                                                                                                                         cctGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGaaa                                                                                                                                               <  2:385/132‑1 (MQ=255)                                                                                                                                                                                                                         cctGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGaaa                                                                                                                                               >  1:385/2‑133 (MQ=255)                                                                                                                                                                                                                              tggCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGTTGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAAGGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCa                     <  1:470419/248‑1 (MQ=255)                                                                                                                                                                                                                                   aTTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGGGTGATGTTCCCACCAGGCCTTTCAAAATTATTGCGGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAAGGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCa               <  1:266906/251‑1 (MQ=255)                                                                                                                                                                                                                                                ttGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCTTGCCGTTCCAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGCAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGAc  >  2:47119/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       AATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGAC  >  NC_000913/12141‑12630

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