Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,567,505	C→G	5.4%	P201P (CCG→CCC)	eutH ←	ethanolamine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,567,505	0	C	G	5.4%	296.0 / 5.2	89	P201P (CCG→CCC)	eutH	ethanolamine transporter
Reads supporting (aligned to +/- strand):  ref base C (40/44);  new base G (3/2);  total (43/46)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAAC  >  NC_000913/2567357‑2567653                                                                                                                                                     |                                                                                                                                                     cccAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCAcc                                                                                                                                                      >  1:92607/1‑149 (MQ=255)  ccAGCCTAGCTGGATTTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCg                                                                                                                                                     <  2:1529382/149‑1 (MQ=255)  ccAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCg                                                                                                                                                     >  1:1205005/1‑149 (MQ=255)     gCCAAGCATGAATTTCACTACCGCAGCGGCAAGACCGAGCGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGTGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGa                                                                                                                                                  <  1:269947/149‑1 (MQ=255)         aGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCACCGAGGACTTTGGCGAAGATCTGGAAGCCGTTGATAATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCAt                                                                                                                                              <  1:111870/149‑1 (MQ=255)          gCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTTGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATAACCGGGATCATg                                                                                                                                             <  2:1114597/149‑1 (MQ=255)           cAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGt                                                                                                                                            <  1:615724/149‑1 (MQ=255)            aGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGtt                                                                                                                                           >  2:450777/1‑149 (MQ=255)            aGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGtt                                                                                                                                           >  2:1002295/1‑149 (MQ=255)              gAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGTGCCACCAGAATCGCAACAATGATCACCGGGATCATGTtca                                                                                                                                         >  2:1425970/1‑149 (MQ=255)                 tttCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTtcatca                                                                                                                                      <  1:1184366/149‑1 (MQ=255)                 tttCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTtcatca                                                                                                                                      <  1:1181033/149‑1 (MQ=255)                    cACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCCCGGTGATCAttgttg                                                                                                                                         >  2:1545106/1‑138 (MQ=255)                    cACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCCCGGTGATCAttgttg                                                                                                                                         <  1:1545106/143‑6 (MQ=255)                    cACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCatcagga                                                                                                                                   <  2:1243593/149‑1 (MQ=255)                       tACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCCCCAGAATCGCAACAATGATCACCGGGATCATGTTCatcaggatca                                                                                                                                >  2:591199/1‑149 (MQ=255)                       tACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCatcaggatca                                                                                                                                >  2:960473/1‑149 (MQ=255)                         ccGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCCGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCatcagtatcagg                                                                                                                              >  2:809201/1‑149 (MQ=255)                         ccGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCatcaggatcagg                                                                                                                              <  1:1341400/149‑1 (MQ=255)                         ccGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCatcaggatcagg                                                                                                                              <  1:369474/149‑1 (MQ=255)                               cGGCAAGACCGAGGGTGATCAATGCACCGAGGAATTTGGCGAAGATCTGGAAGCCTTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGaaa                                                                                                                        <  1:446695/149‑1 (MQ=255)                               cGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGaaa                                                                                                                        <  1:1193760/149‑1 (MQ=255)                                    aGACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGaa                                                                                                                   <  2:906778/149‑1 (MQ=255)                                     gACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGATAGTGAAt                                                                                                                  <  2:543293/149‑1 (MQ=255)                                     gACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAAt                                                                                                                  >  2:119311/1‑149 (MQ=255)                                         gAGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCTCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCa                                                                                                              >  2:1279240/1‑149 (MQ=255)                                         gAGGGTGATCAATGCAACGAGGAATTTGGCGAAGACCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCa                                                                                                              >  1:50843/1‑149 (MQ=255)                                          aGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCAc                                                                                                             >  1:530343/1‑149 (MQ=255)                                          aGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCAc                                                                                                             <  1:15112/149‑1 (MQ=255)                                            ggTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCg                                                                                                           >  2:1597498/1‑149 (MQ=255)                                            ggTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAAGGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCg                                                                                                           >  2:1147511/1‑149 (MQ=255)                                             gTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCgg                                                                                                          <  2:896421/149‑1 (MQ=255)                                             gTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCgg                                                                                                          <  2:1443277/149‑1 (MQ=255)                                             gTGATCAATGCAACGAGGAATTTGGCGAAGAGCTGGAAGCCGTTTATCATGTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCgg                                                                                                          <  2:1132926/149‑1 (MQ=255)                                                    aTGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCg                                                                                                   <  2:1205005/149‑1 (MQ=255)                                                     tGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGt                                                                                                  >  2:954464/1‑149 (MQ=255)                                                        aaCGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCTTTTTTTCCGGGATGAATTTAAGCACCAGCGCCACAAGATTCGCAACACTGATCACGGGGTTCATGGTCATCAGGATCAGGTCGAATGTGAATTCGACCGGCTGTCCTTTGa                                                                                               <  1:228422/149‑1 (MQ=255)                                                         aCGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGAt                                                                                              <  2:1496685/149‑1 (MQ=255)                                                         aCGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGAt                                                                                              <  1:210568/149‑1 (MQ=255)                                                          cGAGGACTGTGGCGAAGAGCTGGAAGCCGTTGATCATTTTTACCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATc                                                                                             <  2:1276885/149‑1 (MQ=255)                                                             ggAATTTGGCGAAGACCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGc                                                                                          <  2:50843/149‑1 (MQ=255)                                                                aTTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCCCCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGcacc                                                                                       >  1:184200/1‑148 (MQ=255)                                                                 tttGGCGAAGATCTGGAAGCAGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCacac                                                                                      <  2:1572312/149‑1 (MQ=255)                                                                   tGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCg                                                                                    <  1:8898/149‑1 (MQ=255)                                                                   tGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCg                                                                                    <  2:92607/149‑1 (MQ=255)                                                                      cGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAg                                                                                 <  1:119311/149‑1 (MQ=255)                                                                       gAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGt                                                                                <  1:1287895/149‑1 (MQ=255)                                                                       gAAGATCTGGAAGCCGTTGATCATTTATTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAAATCCACCGGCTGGCCGTTGATCTGCACACCGGAGt                                                                                >  2:332725/1‑149 (MQ=255)                                                                       gAAGATCTGGAAGCCGTTGATCA‑TTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTa                                                                               >  2:769965/1‑149 (MQ=255)                                                                          gttCTGGAAGCCGTTGATCATTTTTTCCGGGATTAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATTTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGTAGTACa                                                                             <  2:306085/147‑1 (MQ=255)                                                                           aTCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACAt                                                                            <  1:1002295/149‑1 (MQ=255)                                                                           aTCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACAt                                                                            <  1:960473/149‑1 (MQ=255)                                                                            tCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATa                                                                           <  2:1426044/149‑1 (MQ=255)                                                                             cTGGAATCCGGTGATCATTTTTTCCGGGATGAATATCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGTATCATGTTCATCAGGATCAGGCCGATAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAg                                                                          <  2:205179/149‑1 (MQ=255)                                                                              tGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGc                                                                         >  2:1219990/1‑149 (MQ=255)                                                                                gAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGAt                                                                                              >  2:1154432/1‑126 (MQ=255)                                                                                gAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGAt                                                                                              <  1:1154432/126‑1 (MQ=255)                                                                                      gTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGac                                                                 <  2:111611/149‑1 (MQ=255)                                                                                      gTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGac                                                                 >  1:676754/1‑149 (MQ=255)                                                                                      gTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGac                                                                 <  1:954464/149‑1 (MQ=255)                                                                                      gTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGACCTGCACACCGGAGTACATAGCAACCAGac                                                                 >  2:15378/1‑149 (MQ=255)                                                                                          aTCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGaccacc                                                             >  1:497838/1‑149 (MQ=255)                                                                                             aTTTTTTCCGGGATGATTTTCAGCCGCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACGCCGGAGTACGTAGAAACCAGACCACCAGc                                                          >  1:234482/1‑149 (MQ=255)                                                                                             aTTTTTTCCGGGATGAATTTCAGCCCCAGTGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGc                                                          <  1:1425970/149‑1 (MQ=255)                                                                                                  ttCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATAc                                                     <  1:1371735/149‑1 (MQ=255)                                                                                                     cGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACGAc                                                  <  2:1247753/149‑1 (MQ=255)                                                                                                        gATGAATTTCAGCCCCAGCGCCACCAGATTCGCAACAATGATCACCGGGATCATCTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTGCATAGCATCCAGCCCACAATCGATACAACCaa                                               >  1:869669/1‑149 (MQ=255)                                                                                                         aTGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAAt                                              >  2:1570489/1‑149 (MQ=255)                                                                                                           gAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCCCCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCCGCGATACAAGCAAGCg                                            >  1:302461/1‑149 (MQ=255)                                                                                                             aTTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCAGACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGa                                          >  1:428730/1‑149 (MQ=255)                                                                                                               ttCAGCCCCAGCGCCACCAGAATCGCAACAATGATCAGCGGGATCATGTTCCGCAGGATCAGGGCGAAATTGAATTCAACCGGCTGGCCGTTTGTGTGCACACCGGAGTAAATAGAAACCAGACCAACAGCGATAAAACCaaaacgaat                                        >  2:112669/1‑142 (MQ=255)                                                                                                                   gCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTc                                    <  2:926544/149‑1 (MQ=255)                                                                                                                    ccccAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTca                                   <  2:428730/149‑1 (MQ=255)                                                                                                                      ccAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTcaca                                 >  2:1588784/1‑149 (MQ=255)                                                                                                                       cAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGACCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACaa                                <  1:15378/149‑1 (MQ=255)                                                                                                                        aGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAAt                               >  2:1040811/1‑149 (MQ=255)                                                                                                                         gcgcCACCAGAATCGCAACAATGATCACCGGGATCATGGTCATCAGGATCAGGGCGAAAGGGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATg                              <  2:785429/149‑1 (MQ=255)                                                                                                                             caccaGAATCGCAACAATGATCACCGGGCTCAGGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCCTTGATCTGCACACCGGAGTACATAGCAACCAGACCACAAGCGCTACAACGACTCGGACTGGTCACCATCCCCg                          >  2:103880/1‑149 (MQ=255)                                                                                                                               ccaGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGcc                        >  2:1187048/1‑149 (MQ=255)                                                                                                                                  gAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACACCCAATCGGAATGGTCACAATGCccgcccgc                     >  1:139457/1‑149 (MQ=255)                                                                                                                                   aaTCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGcc                    >  1:372309/1‑148 (MQ=255)                                                                                                                                     tCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAAGCGGAATGGTCACAATGCCCGCCAGCACg                  >  1:40365/1‑149 (MQ=255)                                                                                                                                      gatgaacaTGATCCCGGTGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACAAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGc                 >  1:339011/8‑149 (MQ=255)                                                                                                                                        caacaaTGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCg               >  1:24748/1‑149 (MQ=255)                                                                                                                                          acaaTGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAg             >  2:162443/1‑149 (MQ=255)                                                                                                                                          aacaTGATCCCGGTGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTCGCCGTTGATCTGCCCACCGGAGTCCATAGCAACCAGACCACAAGCGATACAACCAATCGGAATGGTCACAATGACCGCCAGCACGCCGAg             >  1:547230/4‑149 (MQ=255)                                                                                                                                            aaTGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTAGATATGCACACAGGAGTACATAGCAACCAGACAACCAGCGATACAACCAATCGGAATGGTCACAATGACCGCCAGAACGCagagag           >  2:1195976/1‑147 (MQ=255)                                                                                                                                                 tCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGa      >  2:1488873/1‑149 (MQ=255)                                                                                                                                                   aCCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATa    >  2:1560837/1‑149 (MQ=255)                                                                                                                                                   aCCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATa    <  1:1597498/149‑1 (MQ=255)                                                                                                                                                     cGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAac  >  2:719659/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGATCTGGAAGCCGTTGATCATTTTTTCCGGGATGAATTTCAGCCCCAGCGCCACCAGAATCGCAACAATGATCACCGGGATCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAAC  >  NC_000913/2567357‑2567653

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