Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	460,060	C→T	I391I (ATC→ATT)	lon →	Lon protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	460,060	0	C	T	92.7%	95.5 / ‑1.2	43	I391I (ATC→ATT)	lon	Lon protease
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (23/15);  total (26/17)
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.23e-01

AGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCG‑ATCCTCTGCC‑TGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAA  >  NC_000913/459820‑460305                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       aGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCAGTATGGGCTGGAGCGCGTGAAAGAGCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCG‑ATCCTCTGCC‑TGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAATCGGAAATTCGTGGTCAcc                                                                                                                                                                                                                                              <  2:46739/251‑1 (MQ=255)         gTCAGGCGCAGGAAATCCTTGATAACGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCT‑ATCCTATGCC‑TGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCCTTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGAGTGCGTGATGACGCGGAAATTCGTGGTCACCGACGTACt                                                                                                                                                                                                                                      >  2:372224/1‑251 (MQ=255)                   ggAAATCCTGGATTCTGACCTTTATGGTCTGGACGGCGTAAAAGACCATATCCTTGTGTAACTGGGGGTGCAAACCCGTG‑CAAAAAAAATAAGGGGCCGGATCCCCTGCC‑TGGTAGGGCCGCCGGGGTGAGGTAAAACCTCTGTGGGTCAGTCCTTTGCCAAACCCACCGGGGGTAAATATGTCCGTATGGCGCTGGGCGGCGTGGGTGATGAAGCGGAAATTCGTGGTCCCCGCCGTACTTACCTCGGtt                                                                                                                                                                                                                            <  2:175188/251‑1 (MQ=255)                              aTACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCG‑ATCCTCTGCC‑TGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGt                                                                                                                                                                                                                 >  1:15218/1‑251 (MQ=255)                              aTACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGAATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCG‑ATCCTCTGCC‑TGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATTCGTGGTCACCGCCGTACTTACCTCGGTTCTATGCCgggg                                                                                                                                                                                                                 >  2:247191/1‑250 (MQ=255)                                 ccGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCG‑ATCCTCTGCC‑TGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATTCGTGGTCACCGCCGTACTTACATCGGTCCTATGCCGGGTaa                                                                                                                                                                                                               >  2:374637/1‑250 (MQ=255)                                                     cgcgTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCG‑ATCCTCTGCC‑TGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGa                                                                                                                                                                                           >  1:165416/1‑250 (MQ=255)                                                      gcgTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCG‑ATCCTCTGCC‑TGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGaa                                                                                                                                                                                          >  1:188592/1‑250 (MQ=255)                                                          gAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCG‑ATCCTCTGCC‑TGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAATTATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGACAATGGCGAAAGTgg                                                                                                                                                                                     >  2:101629/1‑251 (MQ=255)                                                                          gAGTATCTTGCGGTTCAAAGCCGTGTCAACCAAATCAAGGGACCG‑ATCCTCTGCC‑TGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGACCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATTCGTGGTCACCGCCGTCCTTACATCGGTTCTATGCCGGGCAAACTGATCCAGAAAATGGCGAACGTGGGCGTGCAAAAAccgat                                                                                                                                                                     >  2:399734/1‑248 (MQ=255)                                                                                 ttGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCG‑ATCCTCTGCC‑TGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTg                                                                                                                                                             >  1:49707/1‑251 (MQ=255)                                                                                    cGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCG‑ATCCTCTGCC‑TGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTc                                                                                                                                                          >  1:390448/1‑251 (MQ=255)                                                                                     ggTTCAAAGCCGTGTCAACAAAATAAAGGGACCG‑ATCCTCTGCC‑TGGTAGGGCAGCCGGGGATAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATCTGTCCGTAGGGCGCTGGGCGGCGTGCGTGATGACGCGGAAATTAGTGGTCCCAGACGTACTTACATCGGTTCTCTGCCGGGATCACTGAGCCAGAAAACGCCGAAACGGGGCGTG‑CAAAACCCGCTGCTCCTGCTCg                                                                                                                                                         >  2:159852/1‑251 (MQ=255)                                                                                        tCAAAGCAGTGTCAACAAAATCAAGGGACCG‑ATCC‑GTGCCCTGGAAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGTTCAGTCCATTGCCAAATCCACCGGGCTAAAAGATGTCCGTATGGCGCTGGGCGGCGTGCGTGATAAAGCGGAAATTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGTGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGAt                                                                                                                                                       <  1:179163/250‑1 (MQ=255)                                                                                         cAAAGCCGTGTCAACAAAATCAAGGGACCG‑ATCCTCTGCC‑TGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATg                                                                                                                                                      >  1:66699/1‑250 (MQ=255)                                                                                          aatgCCGGGTCAACAAAATCAAGGGACCG‑ACCCTCGGCC‑TGTTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAGTTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATgag                                                                                                                                                    <  1:76826/248‑1 (MQ=255)                                                                                                   tCAACAAAATCATGGGAGCG‑ATCCTCTGCC‑TGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAAAATGTCCGTATGGCGCAGGGCGGCCGGCGTGATGAAGCGCAAATACGTGGTCACCGACGTACTTACGTCGGTTCTATGCCGGGTACACTGATCCTGAAAATGGCGAAAATGGGCGGG‑ACGAACCCACTGTTCCTGCTCGTTGAGATCCACaa                                                                                                                                            >  2:149948/1‑250 (MQ=255)                                                                                                     aacaaAATCAAGGGACCG‑ATCCTCTGCC‑TGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAAt                                                                                                                                         >  1:5845/1‑251 (MQ=255)                                                                                                     aacaaAATCAAGGGACCG‑ATCCTCCGCC‑TGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAAt                                                                                                                                         >  1:117671/1‑251 (MQ=255)                                                                                                                     cG‑ATCCTCTGCC‑TGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGt                                                                                                                         >  1:287913/1‑251 (MQ=255)                                                                                                                                          gccgcgGGGGGTAGGTAAAACCTCTCTTGGGCAGCCCATTGCCAAAGCCCCCGGGGGTAAATATGTCCGTATGGCGCGGGGGGGCGTGCGTGATGAAGCGGAAATTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCAc                                                                                                      <  1:247191/251‑1 (MQ=255)                                                                                                                                             gccgGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAACTATGTCCGTATGGCGCTGGGCGGCGTGCGAGATGAAGCGGAAATTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGAAATGCGTGGCGATCCGGCCTCTGCACTGc                                                                                                   >  1:418501/1‑251 (MQ=255)                                                                                                                                                          aaaaCTTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATAGGTCCGTATGGCGCGGGGCGGCGTGCGTGAGGAAGCGGAAATTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGa                                                                                      <  1:27645/251‑1 (MQ=255)                                                                                                                                                             cccTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATcc                                                                                   <  1:213530/250‑1 (MQ=255)                                                                                                                                                                tctctTGATCAGTCCATTACCAAAGCCACCGGGCGTAATTATGTCCGGCTGGCGCTGTGCGTCTTGCGTGATTCAGCGAAAATTCATTATCCCCGCCGTACTTGCATCGGATATATTCCGATTAAAATGATCCAGCAATTGGCGACTTTGGACGTG‑AAAAACCCGCTGTTCCTGCTCGATGCAATCGACAAACTGACATCTGACATGCATGTAGATCCAGCctctttccttcttgaatagttggctccat                                                                                 >  2:261844/1‑225 (MQ=255)                                                                                                                                                                   tttggtgggTCCATGGCAAAGGCCACCGGGCGTAATTATGCCCGTATGGCGCTGGTCGGCGTGCGTGAGGAAGCGGAAATCCGTGGTCTCCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGc                                                                              <  1:66491/242‑1 (MQ=255)                                                                                                                                                                        tCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACg                                                                        >  1:364002/1‑251 (MQ=255)                                                                                                                                                                        tCAGTCCATTGCCAAAGACACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATTCGTGTTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGAGGGAGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGCTCGACAAAATGTCTTCTGCCATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAAAACg                                                                        >  2:248892/1‑251 (MQ=255)                                                                                                                                                                         cAGTCAATTGCCAAAGCCACCGGGCGTAAATATGTTCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAACTTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGAAATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACgg                                                                       >  1:144272/1‑250 (MQ=255)                                                                                                                                                                                                 ggTAAATATGTCCGTAGGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGa                                               <  1:102622/250‑1 (MQ=255)                                                                                                                                                                                                     aaTATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTg                                           >  1:111715/1‑251 (MQ=255)                                                                                                                                                                                                      atatGCCCGTAGGGCGCTGGGGGGCGTGCGTGATGAAGCGGAAATTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCTGAACGTAGCGTTCAGCGACCACTACCTGGAAGTgg                                          <  1:257875/251‑1 (MQ=255)                                                                                                                                                                                                          gTCCGTCTGGCGCTGGGCGGCGTGCGTGATGTAGCTTAAATTCGGGGTCTCCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGAtt                                       <  2:165416/250‑1 (MQ=255)                                                                                                                                                                                                                ggggCGCTGGGGGGGGTGCGGGAGGAAGTGGAAATTCGTTGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGAtct                                <  2:188592/249‑1 (MQ=255)                                                                                                                                                                                                                    cgcTGGGCGGCGTGCGTGATGAAGCGGAAATTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGc                            >  1:316486/1‑251 (MQ=255)                                                                                                                                                                                                                    cgcTGGGAGGCGTGGTTGATGAAGCGGAAATTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCTACCACTACCTGGAAGTGGATTACGATCTCAGc                            <  1:118290/251‑1 (MQ=255)                                                                                                                                                                                                                     ccTGGGGGGCGTGTGTGATGAGGCGGAAATTCGTGGTCGCCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCg                           <  1:426334/250‑1 (MQ=255)                                                                                                                                                                                                                       tGGGCGGCGTGCGTGATGAAGCGGAAATTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACa                                                                                                                              >  1:3926/1‑150 (MQ=255)                                                                                                                                                                                                                           cggcgTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGTGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTg                      >  1:133651/1‑250 (MQ=255)                                                                                                                                                                                                                               gtgcgtgATGAAGCGGAAATTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGGTCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGCCCACTACCTGGCAGTGGATTACGATCTCAGCGACGTGATGtt                 >  1:44877/1‑251 (MQ=255)                                                                                                                                                                                                                                 gggtgtggaAGGGGAAATTTGTGGACACCGCCGTTCTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTT‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTTACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTg               <  1:276564/244‑1 (MQ=255)                                                                                                                                                                                                                                  cgtgATGAAGCGGAAATTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGt              <  1:70651/251‑1 (MQ=255)                                                                                                                                                                                                                                              gAAATGCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGaa  <  1:277796/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       AGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCG‑ATCCTCTGCC‑TGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTG‑AAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAA  >  NC_000913/459820‑460305

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