Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,468,881	A→G	5.1%	T419T (ACA→ACG)	mgtA →	magnesium transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,468,881	0	A	G	5.1%	253.3 / 3.4	79	T419T (ACA→ACG)	mgtA	magnesium transporter
Reads supporting (aligned to +/- strand):  ref base A (36/39);  new base G (2/2);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.91e-01

CACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAG  >  NC_000913/4468737‑4469029                                                                                                                                                 |                                                                                                                                                     cACTGATAAAACCGGTACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCgg                                                                                                                                                  >  1:93539/1‑149 (MQ=255)  aCTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATATCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGt                                                                                                                                                 >  1:317745/1‑149 (MQ=255)  aCTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGt                                                                                                                                                 <  1:11417/149‑1 (MQ=255)  aCTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGt                                                                                                                                                 <  4:532475/149‑1 (MQ=255)     gATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGTGGTGCt                                                                                                                                              >  3:277484/1‑149 (MQ=255)     gATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCt                                                                                                                                              >  3:454637/1‑149 (MQ=255)      aTAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTc                                                                                                                                             >  1:549204/1‑149 (MQ=255)         aaaCCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGaa                                                                                                                                          <  2:104767/149‑1 (MQ=255)          aaCCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAg                                                                                                                                         >  4:275755/1‑149 (MQ=255)          aaCCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAg                                                                                                                                         <  1:181952/149‑1 (MQ=255)                   ccTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATg                                                                                                                                >  2:67525/1‑149 (MQ=255)                     tGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGaa                                                                                                                              <  4:32762/149‑1 (MQ=255)                       aCGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAaga                                                                                                                            <  2:267241/149‑1 (MQ=255)                       aCGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAaga                                                                                                                            <  4:105718/149‑1 (MQ=255)                           aGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCa                                                                                                                        <  2:160715/149‑1 (MQ=255)                             gATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACGGCGGTGCTCGAAGGTACGGATGAAGAGTCAgc                                                                                                                      >  3:504782/1‑149 (MQ=255)                             gATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACGGCGGTGCTCGAAGGTACGGATGAAGAGTCAgc                                                                                                                      >  1:99330/1‑149 (MQ=255)                                      gtgCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGcgctcgct                                                                                                             <  1:120124/149‑1 (MQ=255)                                       tgCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAGAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGa                                                                                                                                  >  4:363177/1‑127 (MQ=255)                                       tgCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGATCGATACAGCGGTGCTCGAAGGTACGGa                                                                                                                                  <  3:363177/127‑1 (MQ=255)                                        gCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTgg                                                                                                           <  4:313783/149‑1 (MQ=255)                                        gCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTgg                                                                                                           >  4:523715/1‑149 (MQ=255)                                        gCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTgg                                                                                                           <  2:89066/149‑1 (MQ=255)                                           ggAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCa                                                                                                        >  3:27572/1‑149 (MQ=255)                                                 tCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGtt                                                                                                  >  3:495967/1‑149 (MQ=255)                                                  cATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTg                                                                                                 <  4:121016/149‑1 (MQ=255)                                                   aTACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTgg                                                                                                >  2:415479/1‑149 (MQ=255)                                                    tACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCGTAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGc                                                                                               >  4:468013/1‑149 (MQ=255)                                                       cGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCaaa                                                                                            >  1:508448/1‑149 (MQ=255)                                                       cGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCaaa                                                                                            >  1:582537/1‑149 (MQ=255)                                                            tctcCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAAtt                                                                                       <  1:415479/149‑1 (MQ=255)                                                            tctcCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAAtt                                                                                       >  2:276005/1‑149 (MQ=255)                                                              tcCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATtga                                                                                     >  1:198671/1‑149 (MQ=255)                                                               ccGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATtgat                                                                                    <  3:522238/149‑1 (MQ=255)                                                                        ccAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCg                                                                           <  2:84029/149‑1 (MQ=255)                                                                                cgcgTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGAttt                                                                   >  4:474144/1‑149 (MQ=255)                                                                                 gcgTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTc                                                                  <  1:67525/149‑1 (MQ=255)                                                                                   gTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGa                                                                <  3:252353/149‑1 (MQ=255)                                                                                   gTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGa                                                                >  2:574826/1‑149 (MQ=255)                                                                                    tgctgcATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATCCAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAg                                                               <  2:416544/149‑1 (MQ=255)                                                                                    tgctgcATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAg                                                               <  3:287827/149‑1 (MQ=255)                                                                                          aTAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGc                                                         <  4:151297/149‑1 (MQ=255)                                                                                            aGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCg                                                       >  3:60063/1‑149 (MQ=255)                                                                                             gcgcGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCgg                                                      <  2:397877/149‑1 (MQ=255)                                                                                               gcgTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGAt                                                    >  4:261835/1‑149 (MQ=255)                                                                                                   ggTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCg                                                <  4:454637/149‑1 (MQ=255)                                                                                                   ggTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCg                                                >  3:79918/1‑149 (MQ=255)                                                                                                      tGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACGGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCggtg                                             <  4:504782/149‑1 (MQ=255)                                                                                                      tGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACGGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCggtg                                             <  2:99330/149‑1 (MQ=255)                                                                                                      tGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCggtg                                             <  1:266129/149‑1 (MQ=255)                                                                                                      tGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCggtg                                             <  2:549204/149‑1 (MQ=255)                                                                                                      tGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCggtg                                             <  4:60063/149‑1 (MQ=255)                                                                                                         aCAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTa                                          >  1:158114/1‑149 (MQ=255)                                                                                                         aCAGCCATTATCAGACCGGACTTAAAAACCTGCTCGATAAAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTAGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGTGCGTCGCCGGATGTCGGTGGTa                                          <  4:353645/149‑1 (MQ=255)                                                                                                            gTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTTGTAGTg                                       >  1:119448/1‑149 (MQ=255)                                                                                                            gTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTg                                       >  1:495314/1‑149 (MQ=255)                                                                                                             tCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCggt                                              <  4:144844/141‑1 (MQ=255)                                                                                                             tCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCggt                                              >  3:144844/1‑141 (MQ=255)                                                                                                             tCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTgg                                      <  2:317745/149‑1 (MQ=255)                                                                                                                 tATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGa                                  <  1:59610/149‑1 (MQ=255)                                                                                                                  aTCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGaa                                 <  4:5125/149‑1 (MQ=255)                                                                                                                    cAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCgg                                               <  1:129196/133‑1 (MQ=255)                                                                                                                    cAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCgg                                               >  2:129196/1‑133 (MQ=255)                                                                                                                      gACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATa                             >  4:304675/1‑149 (MQ=255)                                                                                                                      gACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATa                             >  4:430226/1‑149 (MQ=255)                                                                                                                         cGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCg                          >  3:393738/1‑149 (MQ=255)                                                                                                                          ggACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGa                         >  4:230109/1‑149 (MQ=255)                                                                                                                              ttAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCAc                     >  4:275941/1‑149 (MQ=255)                                                                                                                                aaaaaCCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATa                             <  3:215899/139‑1 (MQ=255)                                                                                                                                aaaaaCCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATa                             >  4:215899/1‑139 (MQ=255)                                                                                                                                aaaaaCCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCa                   >  3:134650/1‑149 (MQ=255)                                                                                                                                aaaaaCCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCa                   >  3:395216/1‑149 (MQ=255)                                                                                                                                 aaaaCCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTGGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCAt                  <  2:508448/149‑1 (MQ=255)                                                                                                                                  aaaCCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCAtg                 <  2:495314/149‑2 (MQ=255)                                                                                                                                   aaCCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCa                <  3:436970/149‑1 (MQ=255)                                                                                                                                   aaCCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCa                <  3:275755/149‑1 (MQ=255)                                                                                                                                       tGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTg            >  3:147476/1‑149 (MQ=255)                                                                                                                                           cGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGttt        <  1:68716/149‑1 (MQ=255)                                                                                                                                                cAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCaaa   <  2:198671/149‑1 (MQ=255)                                                                                                                                                 aGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAg  <  3:411306/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     CACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAG  >  NC_000913/4468737‑4469029

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