Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,486,869	T→A	L218M (TTG→ATG)	lptF →	lipopolysaccharide transport system protein LptF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,486,869	0	T	A	98.4%	202.1 / ‑3.3	61	L218M (TTG→ATG)	lptF	lipopolysaccharide transport system protein LptF
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (29/31);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.47e-01

GAAAGCGAACCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGCGTATCACGTTGGTATTCACCGTGTTTATGATGGCACTTATGGTCGTACCGCTGAGCGTGGTTAACCCACGTCAGGGACGCGTACTGTCGATGC  >  NC_000913/4486622‑4487115                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       gaaagcgaaCCCTGGCATGGCGGCGCTGGCGCCAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGGTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGCAGGCACTGCAATGt                                                                                                                                                                                                                                                     >  2:268258/1‑251 (MQ=255)   aagcgaaCCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGtt                                                                                                                                                                                                                                                    >  1:47679/1‑250 (MQ=255)       gaaCCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCCGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTg                                                                                                                                                                                                                                               >  2:142718/1‑251 (MQ=255)       gaaCCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGGCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTg                                                                                                                                                                                                                                               <  2:100078/251‑1 (MQ=255)          cccTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCTACTAATGGCAGCGCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTTGCCGATTCCGGTCATTTAACCCAGCTTCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGAtt                                                                                                                                                                                                                                            <  2:209579/251‑1 (MQ=255)                      ggcgCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACgg                                                                                                                                                                                                                                >  1:128647/1‑251 (MQ=255)                         gcTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACt                                                                                                                                                                                                                             >  2:150295/1‑251 (MQ=255)                          cTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGACAGCGTTGACGGCAGCGATTTCCAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACtt                                                                                                                                                                                                                            >  1:187064/1‑251 (MQ=255)                                cAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCCGTGCTATTCATCGCAAGCGTTGAAGGCCGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATAAACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCAAGCTGCGCGACGGCTCCCAGGTCGTCACTCTAAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAgg                                                                                                                                                                                                                       >  1:181018/1‑250 (MQ=255)                                           ccAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCg                                                                                                                                                                                                            >  1:107792/1‑250 (MQ=255)                                                 aGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTg                                                                                                                                                                                                     <  1:250031/251‑1 (MQ=255)                                                  gCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTgg                                                                                                                                                                                                    <  1:16824/251‑1 (MQ=255)                                                              aGCTCGGTGCTGTTCATCGAAAGCGCTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCACATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGt                                                                                                                                                                                        >  2:137205/1‑251 (MQ=255)                                                                          ttCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGAccc                                                                                                                                                                            >  2:242161/1‑251 (MQ=255)                                                                                  aaGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGAt                                                                                                                                                                     >  1:7955/1‑250 (MQ=255)                                                                                   aGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATAc                                                                                                                                                                   >  2:228957/1‑251 (MQ=255)                                                                                          aCGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGCCATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGAcc                                                                                                                                                              >  1:140757/1‑249 (MQ=255)                                                                                             gCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACTATACCGACCAGATg                                                                                                                                                         <  1:163845/251‑1 (MQ=255)                                                                                                   aTTTCAAAGATGTGTTCCTCGAGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATg                                                                                                                                                   >  1:67480/1‑251 (MQ=255)                                                                                                            aTGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGca                                                                                                                                               >  1:138216/1‑246 (MQ=255)                                                                                                            aTGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGca                                                                                                                                               <  2:138216/246‑1 (MQ=255)                                                                                                            aTGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATtg                                                                                                                                          >  2:295936/1‑251 (MQ=255)                                                                                                              gtgtTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATtgtg                                                                                                                                        <  2:228127/251‑1 (MQ=255)                                                                                                              gtgtTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATtgt                                                                                                                                         <  1:117008/250‑1 (MQ=255)                                                                                                                 ttCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGa                                                                                                                                      >  1:55249/1‑250 (MQ=255)                                                                                                                       gcgcAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGGCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGa                                                                                                                               <  2:226104/251‑1 (MQ=255)                                                                                                                          cAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGacac                                                                                                                            <  1:137205/251‑1 (MQ=255)                                                                                                                                cGACCAAAAGGGAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAAGGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGGCCCGAACGAGACCGACCAGGTGGACATGCGCTCATTGTGGAACACTGACACCGATc                                                                                                                       <  1:292501/250‑1 (MQ=255)                                                                                                                                      aaaGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACCCCGATCGTGCTcg                                                                                                                <  2:10989/251‑1 (MQ=255)                                                                                                                                          gTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTTGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCgcgc                                                                                                             >  2:85258/1‑250 (MQ=255)                                                                                                                                            aaTGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCgcgccg                                                                                                           >  2:145288/1‑248 (MQ=255)                                                                                                                                               gCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAgaac                                                                                                        <  2:73321/250‑1 (MQ=255)                                                                                                                                                        tCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTggt                                                                                               <  2:224489/250‑2 (MQ=255)                                                                                                                                                             ggtggtggCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATc                                                                                          >  2:125884/1‑250 (MQ=255)                                                                                                                                                              gtggtggCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCAc                                                                                        >  2:114468/1‑251 (MQ=255)                                                                                                                                                                ggtggCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGt                                                                                      <  1:84225/251‑1 (MQ=255)                                                                                                                                                                   ggCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGTGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTg                                                                                    <  1:35394/250‑1 (MQ=255)                                                                                                                                                                    gCCGATTCCGGACATGTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGGGGAACACTGACACCGAGCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGt                                                                                  <  2:255801/251‑1 (MQ=255)                                                                                                                                                                     ccGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGt                                                                                  >  2:277355/1‑250 (MQ=255)                                                                                                                                                                     ccGATGCCGGACAATTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTAGCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGCCCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTa                                                                                 <  2:162011/251‑1 (MQ=255)                                                                                                                                                                       gagtCCGGACATTTAACCCAGCTGCGCGACGGCGCCCAGGGCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCAGTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGTGCACATTGTGGAAAACTTACACCGATCGTGCTCGCGCAGAACTGAACGGGTGTATctcggttgtatt                                                                               <  2:83525/248‑12 (MQ=255)                                                                                                                                                                        aTTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGacac                                                                                                                            <  2:81634/205‑1 (MQ=255)                                                                                                                                                                        aTTCCGGACATTTAACCCAGCTGAGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTc                                                                              >  2:6790/1‑251 (MQ=255)                                                                                                                                                                          tCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACGCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATAGCGACCAGATGGACATGCGCACATTGTGGGACACGGACTCCGATCGTGCCCGCTCAGAACTGAACTGGTGCATCGCTTTGGTATTCa                                                                             <  2:284871/250‑1 (MQ=255)                                                                                                                                                                               aCATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCgtgt                                                                       >  1:120906/1‑251 (MQ=255)                                                                                                                                                                                     aaCCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTatga                                                                 <  2:174511/251‑1 (MQ=255)                                                                                                                                                                                       cccAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTatgatg                                                               <  1:171198/251‑1 (MQ=255)                                                                                                                                                                                       cccAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTatgatg                                                               <  1:232414/251‑1 (MQ=255)                                                                                                                                                                                              gcgcGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTATGATGGCACtt                                                         <  1:145288/250‑1 (MQ=255)                                                                                                                                                                                                  gACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTATGATGGCACTTATGGt                                                    <  1:150295/251‑1 (MQ=255)                                                                                                                                                                                                         cccAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCCCGTTGGAATTCACCGTGTTTATGATGGCACTTATGGTCGTAcc                                              <  1:232758/250‑1 (MQ=255)                                                                                                                                                                                                            aGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTATGATGGCACTTATGGTCGTACCGCTg                                          <  2:128647/251‑1 (MQ=255)                                                                                                                                                                                                                     ctctCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTATGATGGCACTTATGGTCGTACCGCTGAGCGTGGtt                                 >  1:205904/1‑251 (MQ=255)                                                                                                                                                                                                                          aaCCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTATGATGGCACTTATGGTCGTACCGCTGAGCGTGGTTAAccc                            <  2:58265/251‑1 (MQ=255)                                                                                                                                                                                                                             cAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTATGATGGCACTTATGGTCGTACCGCTGAGCGTGGTTAACCCACg                         >  1:168264/1‑251 (MQ=255)                                                                                                                                                                                                                              aGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTATGATGGCACTTATGGTCGTACCGCTGAGCGTGGTTAACCCACGt                        <  2:67480/251‑1 (MQ=255)                                                                                                                                                                                                                                           cGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTATGATGGCACTTATGGTCGTACCGCTGAGCGTGGTTAACCCACGTCAGGGACGCGTAc           <  1:142718/251‑1 (MQ=255)                                                                                                                                                                                                                                            gAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTATGATGGCACTTATGGTCGTACCGCTGAGCGTGGTTAACCCACGTCAGGGACGCGTACt          >  2:59065/1‑251 (MQ=255)                                                                                                                                                                                                                                              aGGCACTGCAATGTTACGTGATTTACGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACTTTGGGATTCGCCGTGTTTATGATGGCACTTATGGTCGTACCTCTGAGCGTGGTTAACCCACGTCAGGGACGCGTACTGt        <  2:129167/251‑1 (MQ=255)                                                                                                                                                                                                                                                    tGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTATGATGGCACTTATGGTCGTACCGCTGAGCGTGGTTAACCCACGTCAGGGACGCGTACTGTCGAtgc  >  2:229550/1‑251 (MQ=255)                                                                                                                                                                                                                                                    tGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTATGATGGCACTTATGGTCGTACCGCTGAGCGTGGTTAACCCACGTCAGGGACGCGTACTGTCGAtgc  >  2:164652/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       GAAAGCGAACCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGCGTATCACGTTGGTATTCACCGTGTTTATGATGGCACTTATGGTCGTACCGCTGAGCGTGGTTAACCCACGTCAGGGACGCGTACTGTCGATGC  >  NC_000913/4486622‑4487115

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