Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,629,371	T→A	5.4%	F28I (TTC→ATC)	yfgG →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,629,371	0	T	A	5.4%	264.6 / 4.4	74	F28I (TTC→ATC)	yfgG	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (38/32);  new base A (2/2);  total (40/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

TTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACTTTCACAGAACATAACGGCACT  >  NC_000913/2629224‑2629519                                                                                                                                                    |                                                                                                                                                     ttttGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTAtctt                                                                                                                                                     <  4:13850/149‑1 (MQ=255)     gTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTAtcttcttc                                                                                                                                                 <  2:245809/149‑1 (MQ=255)      tCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTAtcttcttct                                                                                                                                                <  2:230971/149‑1 (MQ=255)      tCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTAtcttc                                                                                                                                                    <  2:218794/145‑1 (MQ=255)      tCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTAtcttc                                                                                                                                                    >  1:218794/1‑145 (MQ=255)          gttgtGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTgg                                                                                                                                            >  1:658850/1‑149 (MQ=255)              tgCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGt                                                                                                                                        <  2:187133/149‑1 (MQ=255)                ccTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGttt                                                                                                                                      >  2:267540/1‑149 (MQ=255)                ccTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGttt                                                                                                                                      >  3:307398/1‑149 (MQ=255)                 cTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGtttt                                                                                                                                     >  3:314738/1‑149 (MQ=255)                 cTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGtttt                                                                                                                                     <  3:398641/149‑1 (MQ=255)                    aaTGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATc                                                                                                                                  >  4:414127/1‑149 (MQ=255)                      tGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTa                                                                                                                                <  3:324427/149‑1 (MQ=255)                       gTAAGTACCGTCCTCAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTAc                                                                                                                               <  4:588176/149‑1 (MQ=255)                       gTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTAc                                                                                                                               >  3:506720/1‑149 (MQ=255)                         aaGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTc                                                                                                                             <  2:378381/149‑1 (MQ=255)                         aaGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTc                                                                                                                             <  2:249612/149‑1 (MQ=255)                             aCCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTcgtcc                                                                                                                         <  1:647073/149‑1 (MQ=255)                               cGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTcgtccgt                                                                                                                       >  4:626504/1‑149 (MQ=255)                                  ccACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCgg                                                                                                                    >  4:608978/1‑149 (MQ=255)                                    acaGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTg                                                                                                                  >  3:400914/1‑149 (MQ=255)                                    acaGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTg                                                                                                                  >  1:397921/1‑149 (MQ=255)                                    acaGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTg                                                                                                                  >  2:426679/1‑149 (MQ=255)                                       gCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCt                                                                                                               <  3:626504/149‑1 (MQ=255)                                           gggACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATc                                                                                                                                  >  3:131883/1‑126 (MQ=255)                                           gggACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATc                                                                                                                                  <  4:131883/126‑1 (MQ=255)                                              aCATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGgcagca                                                                                                        >  1:390719/1‑149 (MQ=255)                                                    ttCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTAtcttc                                                                                                                                                    <  2:150613/99‑1 (MQ=255)                                                    ttCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTAtcttc                                                                                                                                                    >  1:150613/1‑99 (MQ=255)                                                    ttCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCa                                                                                                  >  4:599897/1‑149 (MQ=255)                                                    ttCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCa                                                                                                  >  2:76237/1‑149 (MQ=255)                                                    ttCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCa                                                                                                  >  1:135869/1‑149 (MQ=255)                                                    ttCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCa                                                                                                  >  1:30266/1‑149 (MQ=255)                                                    ttCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCa                                                                                                  >  2:130705/1‑149 (MQ=255)                                                      cAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCaga                                                                                                <  4:400914/149‑1 (MQ=255)                                                       aaGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCagag                                                                                               >  1:22056/1‑149 (MQ=255)                                                        aGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGc                                                                                              >  3:150174/1‑149 (MQ=255)                                                            aCCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCaaaa                                                                                          >  2:612043/1‑149 (MQ=255)                                                              cTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCaaaaaa                                                                                        <  4:314738/149‑1 (MQ=255)                                                              cTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCATCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCt                                                                                                               <  4:301403/126‑1 (MQ=255)                                                              cTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCATCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCt                                                                                                               >  3:301403/1‑126 (MQ=255)                                                                    tgAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCt                                                                                  >  2:618844/1‑149 (MQ=255)                                                                      aGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCa                                                                                <  2:390719/149‑1 (MQ=255)                                                                       gTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAg                                                                               <  1:618844/149‑1 (MQ=255)                                                                              aaCCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCca                                                                        >  3:33821/1‑149 (MQ=255)                                                                                ccAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCcaca                                                                      <  3:401488/149‑1 (MQ=255)                                                                                        cGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGt                                                              <  1:76237/149‑1 (MQ=255)                                                                                        cGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGt                                                              <  1:130705/149‑1 (MQ=255)                                                                                           aaaCGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGa                                                           >  3:385516/1‑149 (MQ=255)                                                                                                acacCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCAc                                                      <  4:430099/149‑1 (MQ=255)                                                                                                acacCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCAc                                                      <  1:612043/149‑1 (MQ=255)                                                                                                    cGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGt                                                  >  1:35934/1‑149 (MQ=255)                                                                                                         tAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAAc                                             <  4:307398/149‑1 (MQ=255)                                                                                                          aaCAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACg                                            >  1:355600/1‑149 (MQ=255)                                                                                                           aCAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGt                                           >  4:358622/1‑149 (MQ=255)                                                                                                           aCAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGt                                           <  4:488661/149‑1 (MQ=255)                                                                                                            cAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGtt                                          <  3:358622/149‑1 (MQ=255)                                                                                                               tCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGc                                       >  3:132216/1‑149 (MQ=255)                                                                                                                 gCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCgg                                     >  2:190330/1‑149 (MQ=255)                                                                                                                   aTGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGtt                                   <  2:658850/149‑1 (MQ=255)                                                                                                                     gACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTAc                                 <  3:599897/149‑1 (MQ=255)                                                                                                                           tATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTtctc                           >  4:245763/1‑149 (MQ=255)                                                                                                                                tACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACtt                      <  4:132216/149‑1 (MQ=255)                                                                                                                                 aCTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACttt                     <  4:150174/149‑1 (MQ=255)                                                                                                                                  cTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACTTTc                    >  1:266473/1‑149 (MQ=255)                                                                                                                                    gCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCTATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGtt                                   >  1:320823/1‑132 (MQ=255)                                                                                                                                    gCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGtt                                   <  2:320823/132‑1 (MQ=255)                                                                                                                                    gCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACTTTcac                  <  2:397921/149‑1 (MQ=255)                                                                                                                                      tcatcaGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTtc                             <  1:158597/136‑1 (MQ=255)                                                                                                                                      tcatcaGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACTTTCACAg                >  3:395610/1‑149 (MQ=255)                                                                                                                                      tcatcaGCTGTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTtc                             >  2:158597/1‑136 (MQ=255)                                                                                                                                        atcaGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACTTTCACAGaa              >  3:555457/1‑149 (MQ=255)                                                                                                                                        atcaGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACTTTCACAGaa              <  1:426679/149‑1 (MQ=255)                                                                                                                                        atcaGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACTTTCACAGaa              <  4:555457/149‑1 (MQ=255)                                                                                                                                         tcaGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACTTTCACAGAAc             >  1:307890/1‑149 (MQ=255)                                                                                                                                         tcaGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACTTTCACAGAAc             <  4:506720/149‑1 (MQ=255)                                                                                                                                             cTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACTTTCACAGAACATaa         >  2:518706/1‑149 (MQ=255)                                                                                                                                             cTTTATCATCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCAc                                                      >  1:151949/1‑104 (MQ=255)                                                                                                                                             cTTTATCATCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCAc                                                      <  2:151949/104‑1 (MQ=255)                                                                                                                                                 aTCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACTTTCACAGAACATAACGGc     >  1:74543/1‑149 (MQ=255)                                                                                                                                                 aTCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACTTTCACAGAACATAACGGc     <  2:74543/149‑1 (MQ=255)                                                                                                                                                  tcttcttctTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGt                                           >  2:533467/1‑110 (MQ=255)                                                                                                                                                  tcttcttctTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGt                                           <  1:533467/110‑1 (MQ=255)                                                                                                                                                  tcttcttctTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACTTTCACAGAACATAACGGCa    <  3:414127/149‑1 (MQ=255)                                                                                                                                                  tcttcttctTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACTTTCACAGAACATAACGGCa    <  1:190330/149‑1 (MQ=255)                                                                                                                                                   cttcttcttTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACTTTCACAGAACATAACGGCAc   <  4:395610/149‑1 (MQ=255)                                                                                                                                                   cttcttcttTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACTTTCACAGAACATAACGGCAc   <  2:22056/149‑1 (MQ=255)                                                                                                                                                    ttcttcttTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACTTTCACAGAACATAACGGCACt  >  2:190502/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACTTTCACAGAACATAACGGCACT  >  NC_000913/2629224‑2629519

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