Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,000,304	Δ1 bp	8.0%	coding (753/765 nt)	yigE ←	DUF2233 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,000,304	0	C	.	8.0%	297.2 / 19.1	88	coding (753/765 nt)	yigE	DUF2233 family protein
Reads supporting (aligned to +/- strand):  ref base C (40/41);  new base . (2/5);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

GGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATT  >  NC_000913/4000159‑4000450                                                                                                                                                  |                                                                                                                                                    ggTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATACAATAAATTGCAACGTCATGTTGGCCGGATAGGGCGGTTACGCCGCATCCGGCATCTGCGCCATCTTTAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCg                                                                                                                                                  >  3:243913/1‑149 (MQ=255)         gcACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCG‑TCGAGCGGaaag                                                                                                                                         <  4:196928/149‑2 (MQ=255)         gcACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCG‑TCGAGCGGaaag                                                                                                                                         <  2:478611/149‑2 (MQ=255)          cACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCG‑TCGAGCGGaaagg                                                                                                                                        <  2:69647/149‑3 (MQ=255)            ccTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCG‑TCGAGCGGaaaggt                                                                                                                                       >  3:564873/1‑145 (MQ=255)            ccTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCG‑TCGAGCGGaaaggt                                                                                                                                       <  4:389684/148‑4 (MQ=255)            ccTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCG‑TCGAGCGGaaaggt                                                                                                                                       >  3:389684/1‑145 (MQ=255)            ccTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCG‑TCGAGCGGaaaggt                                                                                                                                       <  4:564873/148‑4 (MQ=255)             cTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATg                                                                                                                                     >  4:120402/1‑149 (MQ=255)             cTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATg                                                                                                                                     >  1:511791/1‑149 (MQ=255)             cTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATg                                                                                                                                     >  3:264134/1‑149 (MQ=255)               tATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGt                                                                                                                                   <  2:99290/149‑1 (MQ=255)               tATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGt                                                                                                                                   <  1:74992/149‑1 (MQ=255)                 tttGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATAATGGTGa                                                                                                                                 <  2:313945/149‑1 (MQ=255)                  ttGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGAc                                                                                                                                >  2:29245/1‑149 (MQ=255)                   tGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACa                                                                                                                               >  3:289121/1‑149 (MQ=255)                     gCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACaaa                                                                                                                             >  1:310005/1‑149 (MQ=255)                      ccTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACaaaa                                                                                                                            >  1:521552/1‑149 (MQ=255)                       cTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAg                                                                                                                           <  3:143320/149‑1 (MQ=255)                          aaaaaTCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGAt                                                                                                                        >  2:40154/1‑149 (MQ=255)                          aaaaaTCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGAt                                                                                                                        >  4:344891/1‑149 (MQ=255)                                  ttCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGc                                                                                                                >  2:133586/1‑149 (MQ=255)                                  ttCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGc                                                                                                                <  1:53446/149‑1 (MQ=255)                                   tCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGcc                                                                                                               <  2:261807/149‑1 (MQ=255)                                       aTTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACgg                                                                                                           <  4:434231/149‑1 (MQ=255)                                          cAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGAt                                                                                                        >  2:120514/1‑149 (MQ=255)                                            aTAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTg                                                                                                      >  1:90551/1‑149 (MQ=255)                                             tAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGc                                                                                                     >  1:265797/1‑149 (MQ=255)                                               aaTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGccc                                                                                                   <  4:467905/149‑1 (MQ=255)                                               aaTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGccc                                                                                                   <  4:197127/149‑1 (MQ=255)                                                aTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGcccc                                                                                                  <  3:47795/149‑1 (MQ=255)                                                aTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGcccc                                                                                                  <  1:368117/149‑1 (MQ=255)                                                    cAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCt                                                                                              <  4:264134/149‑1 (MQ=255)                                                       cGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTc                                                                                           <  3:545987/149‑1 (MQ=255)                                                        gTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCa                                                                                          <  1:120514/149‑1 (MQ=255)                                                          cATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCata                                                                                        <  1:179455/149‑1 (MQ=255)                                                                 gCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACatat                                                                                 <  2:343244/149‑1 (MQ=255)                                                                 gCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACatat                                                                                 <  2:427908/149‑1 (MQ=255)                                                                 gCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACatat                                                                                 <  3:495760/149‑1 (MQ=255)                                                                 gCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACatat                                                                                 <  1:285542/149‑1 (MQ=255)                                                                   cGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATAtgt                                                                               <  3:344891/149‑1 (MQ=255)                                                                     gATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGa                                                                             >  1:510676/1‑149 (MQ=255)                                                                     gATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGa                                                                             >  4:479540/1‑149 (MQ=255)                                                                         gggCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATc                                                                         <  2:521552/149‑1 (MQ=255)                                                                          ggCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCg                                                                        >  1:523565/1‑149 (MQ=255)                                                                           gCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGt                                                                       >  4:266094/1‑149 (MQ=255)                                                                             gTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGc                                                                     <  2:255240/149‑1 (MQ=255)                                                                                      gcATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGAt                                                            >  2:571462/1‑149 (MQ=255)                                                                                       cATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATa                                                           <  2:129851/149‑1 (MQ=255)                                                                                            ggCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAagca                                                      >  4:514280/1‑149 (MQ=255)                                                                                             gCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATaaa                                                         <  3:433396/145‑1 (MQ=255)                                                                                             gCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATaaa                                                         >  4:433396/1‑145 (MQ=255)                                                                                                  tGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAgctgct                                                >  4:66811/1‑149 (MQ=255)                                                                                                        aTCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCg                                                                        <  3:115301/119‑1 (MQ=255)                                                                                                        aTCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCg                                                                        >  4:177722/1‑119 (MQ=255)                                                                                                        aTCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCg                                                                        <  3:177722/119‑1 (MQ=255)                                                                                                        aTCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCg                                                                        >  4:115301/1‑119 (MQ=255)                                                                                                          cTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTc                                        <  2:101193/149‑1 (MQ=255)                                                                                                           ttCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCa                                       >  3:17583/1‑149 (MQ=255)                                                                                                           ttCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCa                                       <  1:245198/149‑1 (MQ=255)                                                                                                           ttCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCa                                       >  1:392675/1‑149 (MQ=255)                                                                                                              aGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATCCATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGtt                                    <  4:243913/149‑1 (MQ=255)                                                                                                              aGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGtt                                    <  4:17583/149‑1 (MQ=255)                                                                                                               gTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGttt                                   >  1:363609/1‑149 (MQ=255)                                                                                                               gTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGttt                                   <  1:45736/149‑1 (MQ=255)                                                                                                                tATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTc                                  <  1:395009/149‑1 (MQ=255)                                                                                                                 aTCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCg                                 >  1:106530/1‑149 (MQ=255)                                                                                                                    tGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCtt                              >  3:123221/1‑149 (MQ=255)                                                                                                                     gACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCttt                             <  1:127838/149‑1 (MQ=255)                                                                                                                      acacAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTg                            >  1:63287/1‑149 (MQ=255)                                                                                                                          aaaaCTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCAt                        <  3:529566/149‑1 (MQ=255)                                                                                                                          aaaaCTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCAt                        >  3:155080/1‑149 (MQ=255)                                                                                                                          aaaaCTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCAt                        >  1:282709/1‑149 (MQ=255)                                                                                                                          aaaaCTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCAt                        >  3:534292/1‑149 (MQ=255)                                                                                                                           aaaCTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATa                       >  3:92697/1‑149 (MQ=255)                                                                                                                            aaCTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATaa                      <  2:523565/149‑1 (MQ=255)                                                                                                                              cTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACa                    <  1:492898/149‑1 (MQ=255)                                                                                                                                aTCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAgg                  >  4:18900/1‑149 (MQ=255)                                                                                                                                     tttAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCaaaa             <  2:510676/149‑1 (MQ=255)                                                                                                                                     tttAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCaaaa             <  2:106530/149‑1 (MQ=255)                                                                                                                                      ttAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAAt            >  1:126732/1‑149 (MQ=255)                                                                                                                                          ccTTTCCGCTCGA‑CGGATATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATa        >  1:24436/1‑149 (MQ=255)                                                                                                                                          ccTGTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATa        >  4:299912/1‑149 (MQ=255)                                                                                                                                           cTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTc                                  >  4:354680/1‑122 (MQ=255)                                                                                                                                           cTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTc                                  <  3:354680/122‑1 (MQ=255)                                                                                                                                            tttCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATaaa      >  1:401997/1‑149 (MQ=255)                                                                                                                                             ttCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATaaaa     <  3:479540/149‑1 (MQ=255)                                                                                                                                                cGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAAtt  <  2:63287/149‑1 (MQ=255)                                                                                                                                                  cTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAAc                                            <  1:291310/105‑1 (MQ=255)                                                                                                                                                  cTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAAc                                            >  2:291310/1‑105 (MQ=255)                                                                                                                                                  |                                                                                                                                                    GGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATT  >  NC_000913/4000159‑4000450

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: ‑