Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	664,542	G→A	L217F (CTT→TTT)	rlpA ←	rare lipoprotein RlpA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	664,542	0	G	A	97.7%	136.7 / NA	43	L217F (CTT→TTT)	rlpA	rare lipoprotein RlpA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/24);  minor base C (1/0);  total (19/24)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GGCCCCGACTTGCACCATAAAGTTGCCGCTGGCGCTTTGCGAGACGGCTTGCGGTGTCACCATTGCAGGCGAGGTTGCAGGCGTTGTCGACGGAGCTGTAACAACGGGCTGTGGCGCAGGCGTCGGTTCGCTGCCTTCCAGTACACCAGGCGCTAAGGTCGTTGGTGCGCCGAGGAAACCGCTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAGAATGTCACCCTGCGGGCCAGACACTGAACTTGTTCCCGCGCCACCGCTTAAATCGGGAGGTGCAGGCAGGGCGTAAGTCTGTTTGGCGACTGTGGTACAAGCCATACCAGGACCAGAAAGCGAACCATCCTGGGCAACAATAATCGGATCGATACGAACTTTGGTGTTGTTTGACGTGTTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATC  >  NC_000913/664300‑664785                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    gcccccGACTTGCACCATAAAGTTGCCGCTGGCGCTTTGCGAGACGGCTTGCGGTGTCACCATTGCAGGCGAGGTTGCAGGCGTTGTCGACGGAGCTGTAACAACGGGCTGTGGCGCAGGCGTCGGTTCGCTGCCTTCCAGTACACCAGGCGCTAAGGGCGTTGGTGCGCCGAGGAAACCGCTACTGGTTACCGGCGCGCCGGTCGGAGCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAAAATGTCAc                                                                                                                                                                                                                                             >  2:242516/3‑251 (MQ=255)    cccGACTTGCACCATAAAGTTGCCGCTGGCGCTTTGCGAGACGGCTTGCGGGGTCACCATTGCAGGCGAGGTTGCAGGCGTTGTCGACGGAGCTGTAACAACGGGCTGTGGCGCAGGCGTCGGTTCGCTGCCTTCCAGTACACCAGGCGCTAAGGTCGTTGGTGCGCCGAGGAAACCGCTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAAAATGTCACCCt                                                                                                                                                                                                                                          <  2:279513/251‑1 (MQ=255)           tGCACCATAAAGTTGCCGCTGGCGCTTTGCGAGACGGCTTGCGGTGTCACCATTGCAGGCGAGGTTGCAGGCGTTGTCGACGGAGCTGTAACAACGGGCTGTGGCGCAGGCGTCGGTTCGCTGCCTTCCAGTACACCAGGCGCTAAGGTCGTTGGTGCGCCGAGGAAACCGCTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAAAATGTCACCCTGCGGGcc                                                                                                                                                                                                                                   >  2:325082/1‑251 (MQ=255)            gCACCATAAAGTTGCCGCTGGCGCTTTGCGAGACGGCTTGCGGTGTCACCATTGCAGGCGAGGTTGCAGGCGTTGTCGACGGAGCTGTAACAACGGGCTGTGGCGCAGGCGTCGGTTCGCTGCCTTCCAGTACACCAGGCGCTAAGGTCGTTGGTGCGCCGAGGAAACCGCTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAAAATGTCACCCTGCGGGCCa                                                                                                                                                                                                                                  <  1:76775/251‑1 (MQ=255)                cATAAAGTTGCCGCTGGCGCTTTGCGAGACGGCTTGCGGTGTCACCATTGCAGGCGAGGTTGCAGGCGTTGTCGACGGAGCTGTAACAACGGGCTGTGGCGCAGGCGTCGGTTCGCTGCCTTCCAGTACACCAGGCGCTAAGGTCGTTGGTGCGCCGAGGAAACCGCTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAAAATGTCACCCTGCGGGCCAgcc                                                                                                                                                                                                                               >  2:140720/1‑248 (MQ=255)                 aTAAAGTTGCCGCTGGCGCTTTGCGAGACGGCTTGCGGTGTCACCATTGCAGGCGAGGTTGCAGGCGTTGTCGACGGAGCTGTAACAACGGGCTGTGGCGCAGGCGTCGGTTCGCTGCCTTCCAGTACACCAGGCGCTAAGGTCGTTGGTGCGCCGAGGAAACCGCTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAAAATGTCACCCTGCGGGCCAGacac                                                                                                                                                                                                                             >  1:440634/1‑251 (MQ=255)                           cGCTGGCGTATTGCGAGACGGCTTGCGGTGTCCCCATGGCAGGCGAGGTTGCAGGCGTTGTCGACGGAGCTGTAACAACGGGCTGTGGCGCAGGCGTCGGTTCGCTGCCGTCCAGTACACCAGGCGCTAAGGTCGTTGGTGCGCCGAGGAAACCGCTACGGGTGACCGGCGCGCCGGTCGGATCTTCGCTGTTGAGCGTCGAATTACTGACCGGAAAAATGTCACCCTGCGGGCCAGACGCTGAACTTGtt                                                                                                                                                                                                                   <  2:252905/251‑1 (MQ=255)                                              ggCTTGCGGTGTCACCATTGCAGGCGAGGTTGCAGGCGTTGTCGACGGAGCTGTAACAACGGGCTGTGGCGCAGTCGTCGGTTCGCTGCCTTCCAGTACACCAGGCGCTAAGGTCGTTGGTGCGCCGAGGAAACCGCTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAAAATGTCACCCTTCGGGCCAGATACTGAACTTGTTCCCGCGCCACCGCTTAAAt                                                                                                                                                                                                <  1:381188/251‑1 (MQ=255)                                              ggCTTGCGGTGTCACCATTGCAGGCGAGGATGCAGGCGTTGTCGACGGAGCTGTAACAACGGGCTGTGGCGCAGGCGTCGGTTCGCTGCCTTCCAGTACACCAGGCGCTAAGGTCGTTGGTGCGCCGAGGAAACCGCTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAAAATGTCACCCTGCGGGCCAGACACTGAACTTGTTCCCGCGCCACCGCTTAAAt                                                                                                                                                                                                >  2:136626/1‑251 (MQ=255)                                              ggCTTGCGGTGTCACCATGGCAGGCGAGGTTGCAGGCGTTGTCGACGGAGCTGTAACAACGGGCTGTGGCGCAGGCGTCGGTTCGCTGCCTTCCAGTACACCAGGCGCTAAGGTCGTTGGTGCGCCGAGGAAACCGCTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAAAATGTCACCCTGCGGGCCAGACACTGAACTTGTTCCCGCGCCACCGCTTAAAt                                                                                                                                                                                                <  1:199668/251‑1 (MQ=255)                                                         tCACCATTGCAGGCGAGGTTGCAGGCGTTGTCGACGGAGCTGAAACAACGGGCTGTGGCGCAGGCGTCGGTTCGCTGCCTTCCAGTACACCAGGCGCTGAGGTCGTTGGTGCGCCGAGGAAACCGATACTGGTTACCGGCGCGCAGGTCGGATCTTCGCTTTTTATCGTCGAATTACTGACCGGAAAAATGTCACCCTGCGGGACAGACACTGAACTTGTTCCCGCGCCACCGCTTAAATCGGGCGGTgc                                                                                                                                                                                      >  1:263886/1‑250 (MQ=255)                                                          cACCATTGCAGGCGAGGTTGCAGGCGTTGTCGACGGAGCTGTAACAACGGGCTGTGGCGCAGGCGTCGGTTCGCTGCCTTCCAGTACACCAGGCGCTAAGGTCGTTGGTGCGCCGAGGAAACCGCTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAAAATGTCACCCTGCGGGCCAGACACTGAACTTGTTCCCGCGCCACCGCTTAAATCGGGAGGTgca                                                                                                                                                                                     <  2:93950/250‑1 (MQ=255)                                                                 gCAGGCGAGGTTGCAGGCGTTGTCGACGGAGCTGTAACAACGGGCTGTGGCGCAGGCGTCGGTTCGCTGCCTTCCAGTACACCAGGCGCTAAGGTCGTTGGTGCGCCGAGGAAACCGCTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAAAATGTCACCCTGCGGGCCAGACACTGAACTTGTTCCCGCGCCACCGCTTAAATCGGGAGGTGCAGGCAGGGc                                                                                                                                                                             <  1:406727/251‑1 (MQ=255)                                                                  cAGGCGAGGTTGCAGGCGTTGTCGACGGAGCTGTAACAACGGGCTGTGGCGCAGGCGTCGGTTCGCTGCCTTCCAGTACACCAGGCGCTAAGGTCGTTGGTGCGCCGAGGAAACCGCTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAAAATGTCACCCTGCGGGCCAGACACTGAACTTGTTCCCGCGCCACCGCTTAAATCGGGAGGTGCAGGCAGGGCg                                                                                                                                                                            >  1:382653/1‑251 (MQ=255)                                                                   aGGCGAGGTTACAAGCGTTGTCGCCGGAGCTGTAACAACAGGCTGTGGCGCAGGAGTCGGTTCGCTGCCTTCCAGTACACCAGGCGCAACGGTCGTTGGCGCGAAGAGGAAACCGCTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAAAATGTCACCCTGCGGGCCAGACACTGAACTTGTTACCGCGCAACCGCTTAACTCGGAAGGTgcaagacggncg                                                                                                                                                                            >  2:415742/1‑241 (MQ=255)                                                                         ggTTGCAGGCGTTGTCGACGGAGCTGTAACAACGGGGTGTGGCGAAGGCGTCGGTTCGCTGCCTTCCAGTACACCAGGCGCTAAGGTCGTTGGTGCGCCGAGGAAACCGCTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAAAATGTCACCCTGCGGGCCAGACACTGACATTGTTCCCGCGCCACCGCTTAAATCGGGAGGTGCAGCCAGGGCGTAAGTcc                                                                                                                                                                     >  2:333460/1‑250 (MQ=255)                                                                                         gACGGAGCTGTAACAACGGGCTGTGGCGCAGGCGTCGGTTCGCTGCCTTCCAGTACACCAGGCGCTAAGGTCGTTGGTGCGCCGAGGAAACCGCTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAAAATGTCACCCTGCGGGCCAGACACTGAACTTGTTCCCGCGCCACCGCTTAAATCGGGAGGTGCAGGCAGGGCGTAAGTCTGTTTGGCGACTGTGGt                                                                                                                                                     <  1:192724/251‑1 (MQ=255)                                                                                                cTGTAACAACGGGCTGTGGCGCAGGCGTCGGTTCGCTGCCTTCCAGTACACCAGGCGCTAAGGTCGTTGGTGCGCCGAGGAAACCGCTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAAAATGTCACCCTGCGGGCCAGACACTGAACTTGTTCCCGCGCCACCGCTTAAATCGGGAGGTGCAGGCAGGGCGTAGGTCTGTTTGGCGACTGTGGTACAAGcc                                                                                                                                              >  1:340155/1‑251 (MQ=255)                                                                                                 tgtatcagcaggCTGGGGCGCAGGCTTCGGTTCGCTGACTTACAGTACACCAGGCGCTAAGGTCGTTGGTGCGCCGAGGAAACCGCTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAAGATGTCACCCTGCGGGCCAGACACTGAACTTGTTCCCGCGCCACCGCTTAAATCGGGAGGTGCAGGCAGGGCGACAGTCTGATTAGCGACTGTGGTACCAGcc                                                                                                                                              >  2:445454/11‑250 (MQ=255)                                                                                                      caacGGGTTGGGGCGCAGGCGCCGTTTCGCTGCCTTCAAGTACCCCAGGCGCTAAGGTCGTTGATGCGCCGAGGAACCCGCTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTGGAATTACTGACCGGAAAAATGTCACCCTGCGGGCCAGACCCTGAACTTGTTCCCGCGCCACCGCTAAAATCGGGAGGTGCAGGCGGGGCGTAAGTTTGTTTGGCGACTTTGGTACAAGCCATACCa                                                                                                                                        <  2:373175/251‑1 (MQ=255)                                                                                                        acGGGCTGTGGCGCAGGCGTCGGTTCGCTGCCTTCCAGTACACCAGGCGCTAAGGTCGTTGGTGCGCCGAGGAAACCGCTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAAAATGTCACCCTGCGGGCCAGACACTGAACTTGTTCCCGCGCCACCGCTTAAATCGGGAGGTGCAGGCAGGGCGTAAGTCTGTTTGGCGACTGTGGTACAAGCCATACCAgg                                                                                                                                      >  1:69370/1‑251 (MQ=255)                                                                                                                ggccttgggggtcgttccGCTGCCTTCGATTCCACCAGTCGCTAAGTTCATGGGTCCGCCGAGGAACCCGCTACGGTTTACCGGCGGGCGGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGAAAAAATGTCGCCCTGCGGCCCAGACACTGTACTTGTTCCCGGCCCACCGCTTAAATCGGGAGGTGCAGGCAGGGCGTAAGTCGGTTTGGCGACTGTGGTAGAAGCCATACCAGGACCAGaaa                                                                                                                              <  2:86217/234‑1 (MQ=255)                                                                                                                 ggCGCAGGCGTCGGTTCGCTGCCTTCCAGTACACCAGGCGCTAAGGTCGTTGGTGCGCCGAGGAAACCGCTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTGAGCGTCGAATTACTGACCGGAAAAATGTCACCCTGCGGGCCAGACACTGAACTTGTTCCCGCGCCCCCGCTTAAATCGGGAGGTGCAGGCAGGGCGTAAGTCTGTTTGGCGACTGTGGTACAAGCCATACCAGGACCAGAAAg                                                                                                                             <  1:305174/251‑1 (MQ=255)                                                                                                                       ggCGTCGGTTCGCTGCCTTCCAGTACCCCAGGCGCTAAGGTCGTTGGTGCGCCGAGGAAACCGCTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAAAATGTCACCCTGCGGGCCAGACACTGAACTTGTTCCCGCGCCACCGCTTAAATCGGGAGGTGCAGGCAGGGCGTAAGTCTGTTTGGCGACTGTGGTACAAGCCATACCAGGACCAGAAAGCGAAcc                                                                                                                       <  1:107535/251‑1 (MQ=255)                                                                                                                                            aGTACACCAGGCGCTAAGGTCGTTGGTGCGCCGAGGAAACCGCTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAAAATGTCACCCTGCGGGCCAGACACTGAACTTGTTCCCGCGCCACCGCTTAAATCGGGAGGTGCAGGCAGGGCGTAAGTCTGTTTGGCGACTGTGGTACAAGCCATACCAGGACCAGAAAGCGAACCATCCTGGGCAACAATAATCgg                                                                                                  >  1:398075/1‑251 (MQ=255)                                                                                                                                                       cgcTAAGGTCGTTGGTGCGCCGAGGAAACCGCTACTGGTTACGTGCGCGCAGGTCGTATCTTCGCTTTTTAGCGTCGAATTCCTGACCGGAAAAATGTCGCCCGGCGGGCCAGACACTGAACTTGTTCCCGCGCCACCGCTTAAAGCGGGAGGTGCAGGCAGGGCGTACGTCTTTTTGGCGCGTGTGGTACAAGCCATACCAGGCCCAGAAAGCGAACCATCCTGGGCAACAATAATCGGATCGATACGaa                                                                                       <  2:348739/251‑1 (MQ=255)                                                                                                                                                       cgcTAAGGTCGTTGGTGCGCCGAGGAAACCGCTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAAAATGTCACCCTGCGGGCCAGACACTGAACTTGTTCCCGCGCCACCGCTTAAATCGGGAGGTGCAGGCAGGGCGTAAGTCTTTTTGGCGACTGTGGTACAAGCCATACCAGGACCAGAAAGCGAACCATCCTGGGCAACAATAATCGGATCGATACGaa                                                                                       <  1:231883/251‑1 (MQ=255)                                                                                                                                                         cTAAGGACGTTGGTGCGCCGAGGAAACCGCTACTGGCTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAAAATGTCACCCTGCGGGCCAGACACTGAACTTGTTCCCGCGCCACCGCTTAAATCGGGAGGTGCAGGCAGGGCGTAAGTCTGTTTGGCGACTGTGGTACAAGCCATACCAGGACCAGAAAGCGAACCATCCTGGGCAACAATAATCGGATCGATACGAACt                                                                                     >  1:162714/1‑251 (MQ=255)                                                                                                                                                           aaGGTCGTTGGGGCGCCGAGGAAACCGTTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAAAATGTCCCCCTGCGGGCCAGACACTGAACTTGTTCCCGCGCCACCGCTTAAATCGGGAGGTGCAGGCAGGGCGTAAGTCTGTTTGGCGACTGTGGTACAAGCCATACCAGGACCAGAAAGCGAACCATCCTGGGCAACAATAATCGGATCGATACGAACtt                                                                                    <  1:129814/250‑1 (MQ=255)                                                                                                                                                                        cgcggAGGAAACCGCTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAAAATGTCCCCCTGCGGGCCAGACACTGAACGTGTTCCCGCGCCACCGCTTAAATCGGGAGGTGCAGGCAGGGCGTAAGTCTGTTTGGCGACTGGTGTACAAGCCATACCAGGACCAGAAAGCGAACCATCCTGGGCAACAATAATCGGATCGATACGAACTTTGGTGTTGTTTGAc                                                                      <  2:247183/247‑1 (MQ=255)                                                                                                                                                                         gcCGAGGAAACCGCTATGGGTTACCGGCGCGCGGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGGCCGGAAAAATGTCACCCTGGGGGCCAGACACAGAACTTGTTCCCGCGCCACCGCTTAAATCGGGAGGTGCAGGCAGGGCGTAAGTCTGTTTGGCGACTGTGGTACAAGCCATACCAGGACCAGAAAGCGAACCATCCTGGGCAACAATAATCGGATCGATACGAACTTTGGTGTTGTTTGACg                                                                     <  1:389623/251‑1 (MQ=255)                                                                                                                                                                                aaaCCGCTGCTGGTTCCCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCAAATTTCTGACCGGAAAAATGTCACCCTGCGGGCCAGACACTGAACTTGTTCCCGCGCCACCGCTTAAATCGGGAGGTGCAGGCAGGGCGTAAGTCTGTTTGGCGACTGTGGTACAAGCCATACCAGGACCAGAAAGCGAACCATCCTGGGCAACAATAATCGGATCGATACGAACTTTGGTGTTGTTTGACGTGTTAAg                                                              <  2:354332/251‑1 (MQ=255)                                                                                                                                                                                  cccgcttcgtgtttccGGCGCTCCGGTCGGACTTTCGCTTTTTAGCGCCAAATTACTGACCGGAAAAATGTCCCCCTGCGGGCCAGACACTGACCTTGTTCCCGCTCCACCGCTTAAATCGGGGGGTGCAGGCAGGGCGTACGTCTGTTTGGCGACTGTGGTACAAGCCATACCAGGACCAGAAAGCGAACCATCCTGGGCAACAATAAGCGGATCGATACGAACTTTGGTGTTGTTTGACGTGTTAAGAc                                                            <  1:90524/237‑1 (MQ=255)                                                                                                                                                                                       tacgggTTACGGGCGCGCTGGTTGGTTCTTCGCTTTTTAGCGTCGAATTAGTGACCGGAAAAATGTCACCCTGCGGGCCAGACTCTGATCTTTTTCCCGCTCCACCGTTTAAAGCGGGAGTTGCAGGCAGGGCGTAAGTCTGTTTGTCGCTTGTGTTACAAGCCATACCAGGACCAGAAAGCGAACCATCCTGGGCAACAATAATCGGATCGATACGAACTTTGGTGTTGTTTGACGTGTTAAGACGGTc                                                        <  1:458063/246‑1 (MQ=255)                                                                                                                                                                                            gttgaaggcgcgcaggtggcatttccGCTGTTTAGCGCCGACTTACTGCCCGGAAAAATGTCCCCCGGCGGGCCGGACCTTGAACGTGTCCACGCGCCTCCGTTGAAATGGGGAGGAGCAGGAAGGGAGTAAGTCTGTGTGCCGGCTGCGATACAAGCCATACCCGGAGCAGAAAGCGAACCATCCTGGGCAACAATAAGCGGATCGATACGAACTTTGGTGGTGTTTGACGGGTTAAGACGGTCAGCTGc                                                  <  2:263886/226‑1 (MQ=255)                                                                                                                                                                                                              ggttctgcgttttttAGTGTCGAATTACTGCCCGGAAAAATGTCTCCCTGCGGTCCAGACACTGATCTTGTTCCCGCGCCACCGCTTAAAGCGGGAGGTGCAGGCAGGGCGTAAGTTTTTTTGGCTACTGTGGTACAAGCCATACCAGGACCAGAAAGCGAACCATCCTGGGCAACAATAATCTGATCGATACGAGCTTTGGTGTTGTTTGACGTGTTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAAt                                <  1:260819/241‑1 (MQ=255)                                                                                                                                                                                                                atcttagatttttAGCGTCGAATTCCTGACCGGTAAAATGTCATACTGAGGGCCAGAAACTGAAATTGTTCCCGCGCCACCGCTTAAATAGGGCCGTGCAGGCAGAGCGTAAGTCTGTTTGGCGCCAGAGGTACAAGCCATAACAGGACAAGAAAGCGAAACATCCTGGGCAACCATAATCGGATCGATACGAACTTTGGTGTTGTTTAACGTGTTAAGAAGGTCCGCTGCCGCGCGAGATAGTGAAATa                               >  2:196685/9‑250 (MQ=255)                                                                                                                                                                                                                 tCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAAAATGTCACCCTGCGGGCCAGACACTGAACTTGTTCCCGCGCCACCGCTTAAATCGGGAGGTGCAGGCAGGGCGTAAGTCTGTTTGGCGACTGTGGTACAAGCCATACCAGGACCAGAAAGCGAACCATCCTGGGCAACAATAATCGGATCGATACGAACTTTGGTGTTGTTTGACGTGTTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAAc                             >  2:212363/1‑251 (MQ=255)                                                                                                                                                                                                                         ttttAGCGTCGAATTACTGACCGGAAAAATGTCCCCCTGCGGGCCAGACACTGAACTTGTTCCCGCGCCACCGCTTAAAGCGGGAGGTGCAGGCAGGGCGTAAGTCTGTTTGGCGACTGTGGTACAAGCCATACCAGGACCAGAAAGCGAACCATCCTGGGCAACAATAATCGGATCGATACGAACTTTGGTGTTGTTTGACGTGTTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTc                       <  1:133353/249‑1 (MQ=255)                                                                                                                                                                                                                                   gAATTACTGCCCGGAAAAATGCCACCCGGCGTGCCAGGCACTGAACTTGTTCCCGCGCCACCGCTTAAATCGGGAGGTGCAGGCAGGGCGTAAGTCTGTTTGGCGACTGTGGTACAAGCCATACCAGGACCAGAAAGCGAACCATCCTGGGCAACAATAATCGGATCGATACGAACTTTGGTGTTGTTTGACGTGTTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAgg           <  1:242516/251‑1 (MQ=255)                                                                                                                                                                                                                                    aaTTACTGACCAGATCACTGTCGCCCTGCGGGCCAGACACTGAAATTGTTCCCGCGCCACCGCTTAAATCGGGAGGTGCAGGCCAGGCGCTAGTATGTTTGGCGACTGCGGCCCGCGCCATCCCAGGACCAGAAAGCGAACCCTTCTGGGCAACAATAATCGGATCGAAACGACATTTGGTGATGTATAACGTGTTAAGACGATCAGCTCCCGCACGCCAAAGAGAAATAAcggcggagttgccgtacggg          >  2:56662/1‑233 (MQ=255)                                                                                                                                                                                                                                            actggAAAAATGTCACCCTGCGGGCCAGACACTGAACTTGTTCCCGCGCCACCGCTTAAATCGGGAGGTGCAGGCAGGCCGCAAGTCTGCTTGGCGACTGTGGTACAAGCCATACCAGGACCAGAAAGCGAACCATCCTGGGCAACAATAATCGGATCGATACGAACTTTGGAGGTGATTGACGTGTTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAAAAACGCGGTCGTGGCCGCATGGGCCGCGATc  >  2:311337/4‑251 (MQ=255)                                                                                                                                                                                                                                             ccGGAAAAATGTCACCCTGCGGGCAAGACACTGAACTTGTTCCCGCGCCACCGCTTAAATCGGGAGGTGCAGGCAGGGCGTAAGTCTGTTTGGCGACTGTGGTACAAGCCATACCAGGACCAGAAAGCGAACCATCCTGGGCAACAATAATCGGATCGATACGAACTTTGGTGTTGTTTGACGTGTTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATc  >  2:308591/1‑250 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    GGCCCCGACTTGCACCATAAAGTTGCCGCTGGCGCTTTGCGAGACGGCTTGCGGTGTCACCATTGCAGGCGAGGTTGCAGGCGTTGTCGACGGAGCTGTAACAACGGGCTGTGGCGCAGGCGTCGGTTCGCTGCCTTCCAGTACACCAGGCGCTAAGGTCGTTGGTGCGCCGAGGAAACCGCTACTGGTTACCGGCGCGCCGGTCGGATCTTCGCTTTTTAGCGTCGAATTACTGACCGGAAGAATGTCACCCTGCGGGCCAGACACTGAACTTGTTCCCGCGCCACCGCTTAAATCGGGAGGTGCAGGCAGGGCGTAAGTCTGTTTGGCGACTGTGGTACAAGCCATACCAGGACCAGAAAGCGAACCATCCTGGGCAACAATAATCGGATCGATACGAACTTTGGTGTTGTTTGACGTGTTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATC  >  NC_000913/664300‑664785

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