Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	12,357	T→C	T65T (ACT→ACC)	dnaK →	chaperone protein DnaK

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

ATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCAT  >  NC_000913/12112‑12595                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               aTTACAGACTCACAACCACATGATGACCTAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCGATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTTGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAATACACCCTGt                                                                                                                                                                                                                                            <  2:438388/250‑1 (MQ=255)    aCAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGc                                                                                                                                                                                                                                        <  1:412237/251‑1 (MQ=255)            acaACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGATGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAAc                                                                                                                                                                                                                                >  1:438960/1‑251 (MQ=255)                    atgatgACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGAtt                                                                                                                                                                                                                        >  1:18440/1‑251 (MQ=255)                      gatgaCCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTgg                                                                                                                                                                                                                      >  1:522983/1‑251 (MQ=255)                         gacagcatatatAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGAACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACAACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTggccg                                                                                                                                                                                                                   >  2:616001/7‑248 (MQ=255)                         gaCCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTCTCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTggccg                                                                                                                                                                                                                   >  1:561031/1‑248 (MQ=255)                              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)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               ATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCAT  >  NC_000913/12112‑12595

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