Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,045,401	A→G	5.2%	D155G (GAT→GGT)	yeeJ →	putative adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,045,401	0	A	G	5.2%	241.2 / 2.1	77	D155G (GAT→GGT)	yeeJ	putative adhesin
Reads supporting (aligned to +/- strand):  ref base A (34/39);  new base G (2/2);  total (36/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 = 9.33e-01

TGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTGGGCGTGGATGAAGATTTTAGCCTGAAGAACTCCCA  >  NC_000913/2045255‑2045545                                                                                                                                                   |                                                                                                                                                 tGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTG‑AAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGata                                                                                                                                                <  3:317525/148‑1 (MQ=255)  gATAATGTCCTCCAGGGTGATGAACTGTATGTCACGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGATCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGata                                                                                                                                                <  7:179852/148‑1 (MQ=255)  gATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGata                                                                                                                                                <  8:332360/148‑1 (MQ=255)   aTAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCATAGCAAATCGGGTCTCTGCTCGCCGAAGatat                                                                                                                                               <  2:341985/148‑1 (MQ=255)       tGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGaa                                                                                                                                            >  8:3847/1‑147 (MQ=255)       tGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACa                                                                                                                                          <  5:53053/149‑1 (MQ=255)       tGTCCGCCAGGGGGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCTAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGata                                                                                                                                                <  1:413270/143‑1 (MQ=255)           cGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGa                                                                                                                                      <  7:264537/149‑1 (MQ=255)           cGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGa                                                                                                                                      <  8:410762/149‑1 (MQ=255)                      agaACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCACGCCGCCGGGTAATAGCAGTTACAACCTCGAGCAACAGATAGCCAGTACTTCACGGCAAATCGGGTCTCTGATCTCCGAAGATATGAACAGCGAGCAAGCGGCa                                                                                                                            >  2:344909/2‑148 (MQ=255)                         aCTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAata                                                                                                                        <  8:242068/149‑1 (MQ=255)                          cTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAatat                                                                                                                       >  3:403719/1‑149 (MQ=255)                             gATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGCCAAGCTCGAGCAACAGATAGCCCGTACTTCTCAGCAAATCGGGTCTCTGCTCGCCGAAGATAAGAACAGCGAGCAAGCGGCAAATATGgc                                                                                                                    >  7:409224/1‑149 (MQ=255)                             gATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACCGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGgc                                                                                                                    >  7:444570/1‑149 (MQ=255)                             gATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATg                                                                                                                      <  3:330188/147‑1 (MQ=255)                             gATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATg                                                                                                                      <  4:131673/147‑1 (MQ=255)                             gATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGgc                                                                                                                    >  5:270565/1‑149 (MQ=255)                             gATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGgc                                                                                                                    <  6:270565/149‑1 (MQ=255)                                gTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGcgcg                                                                                                                 >  1:113314/1‑149 (MQ=255)                                 tCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGGCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGt                                                                                                                <  4:314592/149‑1 (MQ=255)                                   ccGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATg                                                                                                                                              <  2:347067/117‑1 (MQ=255)                                   ccGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATg                                                                                                                                              >  1:347067/1‑117 (MQ=255)                                   ccGGCACAAGTTAGTGAAAAAAAATTAACCCAGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCTGCAAATATGCAGCGTg                                                                                                               >  1:80150/1‑148 (MQ=255)                                    cGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGa                                                                                                                                             <  5:210733/117‑1 (MQ=255)                                    cGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGa                                                                                                                                             >  6:210733/1‑117 (MQ=255)                                       cacaAGTTAGTGAAAAAAAATTAACCCAGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATgg                                                                                                          <  8:271303/149‑1 (MQ=255)                                        acaAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATggg                                                                                                         >  1:385594/1‑149 (MQ=255)                                         caAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGc                                                                                                        <  5:390537/149‑1 (MQ=255)                                            gTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGat                                                                                                                                                 >  8:21809/1‑105 (MQ=255)                                             ttAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACCACCTCGAGCAACAGATCGCCAGTACTTCACAGCAAATCGGGTCGCTGCTCGCCGAAGATATCACCAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCctct                                                                                                    <  6:127121/149‑1 (MQ=255)                                                 tGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGGTATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCtcttct                                                                                                 >  7:298703/1‑148 (MQ=255)                                                 tGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGGTATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCtcttct                                                                                                 >  7:306042/1‑148 (MQ=255)                                                   aaaaaaaaTTAACCCCGCCGCCGGTTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCa                                                                                               <  2:113314/148‑1 (MQ=255)                                                   aaaaaaaaTTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAg                                                                                              <  7:133753/149‑1 (MQ=255)                                                   aaaaaaaaTTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAg                                                                                              >  5:55085/1‑149 (MQ=255)                                                   aaaaaaaaTTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGAAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAg                                                                                              <  4:388455/149‑1 (MQ=255)                                                     aaaaaaTTAACCCCTCCGCCGGGTAATACCAGTGCCAACCTCGATCAACAGATAGCCAGTACTTCACATCTACTCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGc                                                                                            <  8:409224/149‑1 (MQ=255)                                                     aaaaaaTTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGc                                                                                            >  2:194440/1‑149 (MQ=255)                                                     aaaaaaTTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTATCAGGc                                                                                            >  4:254650/1‑149 (MQ=255)                                                        aaaTTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATTGCGCGTGGATGGGCCTCTTCTCAGGCTTc                                                                                         <  5:325680/149‑1 (MQ=255)                                                        aaaTTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTc                                                                                         <  7:403598/149‑1 (MQ=255)                                                                 ccgccgccgGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAAt                                                                                >  7:35291/1‑149 (MQ=255)                                                                     cgccgGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCCGTACTTCCCAGCAACTCGGGTCTCTGCGCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCGCTTCTCAGGCTTCAGGCGCAATGACa                                                                            <  6:375723/149‑1 (MQ=255)                                                                     cgccgGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACa                                                                            <  7:3847/149‑1 (MQ=255)                                                                      gccgGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGa                                                                              >  8:172560/1‑146 (MQ=255)                                                                            tAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGt                                                                     <  6:55085/149‑1 (MQ=255)                                                                               tAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTaa                                                                  <  4:195373/149‑1 (MQ=255)                                                                                    gTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGcc                                                               >  5:131333/1‑147 (MQ=255)                                                                                    gTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGcc                                                               <  6:131333/147‑1 (MQ=255)                                                                                       aCAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTc                                                          <  1:194440/149‑1 (MQ=255)                                                                                         aaCCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCgg                                                        >  4:46385/1‑149 (MQ=255)                                                                                             tCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTAcc                                                    >  2:143864/1‑149 (MQ=255)                                                                                                aGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGc                                                  >  8:248764/1‑148 (MQ=255)                                                                                                            ccAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGgc                                                                                     <  2:219717/101‑1 (MQ=255)                                                                                                            ccAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGgc                                                                                     >  1:219717/1‑101 (MQ=255)                                                                                                            ccAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTg                                     >  7:72876/1‑149 (MQ=255)                                                                                                               gTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTGGGc                                  >  2:362428/1‑149 (MQ=255)                                                                                                                tACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTCAAGCCGCTTCGGTACCGCAAGAATCAc                                         <  8:195901/141‑1 (MQ=255)                                                                                                                tACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTCAAGCCGCTTCGGTACCGCAAGAATCAc                                         >  7:195901/1‑141 (MQ=255)                                                                                                                   ttCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTGGGCGTg                               >  1:169773/1‑148 (MQ=255)                                                                                                                   ttCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTGGGCGTg                               >  3:268111/1‑148 (MQ=255)                                                                                                                   ttCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCAGCTTCGGTACCGCAAGAATCACGCTGGGCGTg                               >  2:363333/1‑148 (MQ=255)                                                                                                                     cacaGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTGGGCGTgg                              <  2:169773/147‑1 (MQ=255)                                                                                                                     cacaGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTGGGCGTGGa                             <  3:46385/148‑1 (MQ=255)                                                                                                                      acaGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTGGGCGTGGAt                            >  7:31537/1‑148 (MQ=255)                                                                                                                      acaGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTGGGCGTGGAt                            >  1:470614/1‑148 (MQ=255)                                                                                                                        aGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTGGGCGTGGATGaa                         <  8:35291/149‑1 (MQ=255)                                                                                                                         gCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCTATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTGGGCGTGGATGAAg                        >  5:37592/1‑149 (MQ=255)                                                                                                                          cAAATCGGGTCTCTGCTCGCCGAAGGTATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTGGGCGTGGATGAAGa                       <  8:298703/149‑1 (MQ=255)                                                                                                                          cAAATCGGGTCTCTGCTCGCCGAAGGTATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTGGGCGTGGATGAAGa                       <  8:306042/149‑1 (MQ=255)                                                                                                                           aaaTCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTGGGCGTGGATGAAGAt                      <  2:385594/149‑1 (MQ=255)                                                                                                                            aaTCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAATCCGCTTCGGTACCGCAAGAATCACGCTGGGCGTGTATGAAGAtt                     >  8:92290/1‑149 (MQ=255)                                                                                                                                 ggTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTGGGCGTGGATGAAGATTTTAGc                <  4:403719/149‑1 (MQ=255)                                                                                                                                  gTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTGGGCGTGGATGAAGATTTTAGcc               >  5:380896/1‑149 (MQ=255)                                                                                                                                      ctGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTGGGCGTGGATGAAGATTTTAGCCTgaa           <  3:254650/149‑1 (MQ=255)                                                                                                                                        gCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTGGGCGTGGATGAAGATTTTAGCCTgaag          >  4:137050/1‑148 (MQ=255)                                                                                                                                         cTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTGGGCGTGGATGAAGATTTTAGCCTgaa           <  1:143864/146‑1 (MQ=255)                                                                                                                                            gCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTGGGCGTGGATGAAGATTTTAGc                <  8:112195/138‑1 (MQ=255)                                                                                                                                            gCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTGGGCGTGGATGAAGATTTTAGc                >  7:112195/1‑138 (MQ=255)                                                                                                                                            gCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTGGGCGTGGATGAAGATTTTAGCCTGAAGAACTc     <  8:444570/149‑1 (MQ=255)                                                                                                                                             ccGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTGGGCGTGGATGAAGATTTTAGCCTGAAGAACTcc    <  7:248764/149‑1 (MQ=255)                                                                                                                                               gAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGTCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTGGGCGTGGATGAAGATTTTAGCCTGAAGAACTCCCa  <  1:344909/149‑1 (MQ=255)                                                                                                                                               gAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTGGGCGTGGATGAAGATTTTAGCCTGAAGAACTCCCa  <  8:72876/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 TGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGCGTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGCAAGAATCACGCTGGGCGTGGATGAAGATTTTAGCCTGAAGAACTCCCA  >  NC_000913/2045255‑2045545

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