Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,233,375	G→A	5.4%	T448M (ACG→ATG)	nhaB ←	sodium:proton antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,233,375	0	G	A	5.4%	241.6 / 2.2	74	T448M (ACG→ATG)	nhaB	sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base G (38/32);  new base A (2/2);  total (40/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

CGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACA  >  NC_000913/1233230‑1233522                                                                                                                                                  |                                                                                                                                                    cGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAg                                                                                                                                                  >  2:60877/1‑149 (MQ=255)  gTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGc                                                                                                                                                 >  6:139615/1‑149 (MQ=255)      aCTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACg                                                                                                                                             <  5:4178/149‑1 (MQ=255)       cTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCCGCAGGAACAGGAACGCAGCCTGACCGTTCGGAGTAGCGACg                                                                                                                                             >  4:251643/1‑148 (MQ=255)          gACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGa                                                                                                                                           <  3:226463/147‑1 (MQ=255)          gACGCAGAGCAAGCCGACGAGGGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGa                                                                                                                                           <  7:53703/147‑1 (MQ=255)                 agCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGAtt                                                                                                                                  <  6:70937/149‑1 (MQ=255)                  gCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCATAGCGACGGAGGGCAGATTg                                                                                                                                 >  5:299974/1‑149 (MQ=255)                  gCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCATAGCGACGGAGGGCAGATTg                                                                                                                                 >  5:305684/1‑149 (MQ=255)                   cAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTgg                                                                                                                                >  3:60719/1‑149 (MQ=255)                   cAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCCGACCGTTCGGCGTAGCGACGGAGGGCAGATTgg                                                                                                                                >  5:205886/1‑149 (MQ=255)                     aGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTa                                                                                                                              <  1:244475/149‑1 (MQ=255)                      gCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTAc                                                                                                                             >  2:236529/1‑149 (MQ=255)                            gaGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTAt                                                                                                                       <  5:139615/149‑1 (MQ=255)                             aGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTAt                                                                                                                       >  3:364784/1‑148 (MQ=255)                                   aGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATgg                                                                                                                 <  7:57620/148‑1 (MQ=255)                                   aGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGc                                                                                                                <  6:368771/149‑1 (MQ=255)                                   aGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGc                                                                                                                <  1:237666/149‑1 (MQ=255)                                    ggACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGc                                                                                                                >  8:227846/1‑148 (MQ=255)                                      aCGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGAc                                                                                                             <  3:303096/149‑1 (MQ=255)                                         aGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGc                                                                                                          <  1:310348/149‑1 (MQ=255)                                           ggTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCTAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCca                                                                                                        <  5:116799/149‑1 (MQ=255)                                               tAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCcagcag                                                                                                    >  6:79118/1‑149 (MQ=255)                                                aaGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCagcagc                                                                                                   >  7:285023/1‑149 (MQ=255)                                                  ggCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAAt                                                                                                                   <  3:378403/131‑1 (MQ=255)                                                  ggCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAAt                                                                                                                   >  4:378403/1‑131 (MQ=255)                                                    cAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACGGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGt                                                                                               <  5:344591/149‑1 (MQ=255)                                                          cATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGTGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCt                                                                                         >  7:340620/1‑149 (MQ=255)                                                            tCCACACCATGAGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACTGATGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTc                                                                                       <  1:60877/149‑1 (MQ=255)                                                              cacaCCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCaa                                                                                     >  3:105270/1‑149 (MQ=255)                                                               acacCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTa                                                                                                                        <  6:210896/113‑1 (MQ=255)                                                               acacCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTa                                                                                                                        >  5:210896/1‑113 (MQ=255)                                                               acacCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAAc                                                                                    <  7:426474/149‑1 (MQ=255)                                                               acacCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAAc                                                                                    <  8:49713/149‑1 (MQ=255)                                                                 acCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGc                                                                                                                <  6:52022/119‑1 (MQ=255)                                                                 acCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGc                                                                                                                >  5:52022/1‑119 (MQ=255)                                                                   cATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGa                                                                                <  8:285023/149‑1 (MQ=255)                                                                   cATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGa                                                                                >  8:434097/1‑149 (MQ=255)                                                                      gCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTAc                                                                                                                             >  5:155380/1‑101 (MQ=255)                                                                      gCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTAc                                                                                                                             <  6:155380/101‑1 (MQ=255)                                                                          ccATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGACGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAgcgc                                                                          >  7:406467/1‑148 (MQ=255)                                                                              aaGAGAGGCGAATCAATGGCGCGAGTGCTGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATg                                                                                                                  <  8:355736/104‑1 (MQ=255)                                                                                    ggCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCAtt                                                               <  1:236529/149‑1 (MQ=255)                                                                                     gCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCagcagc                                                                                                   >  3:198568/1‑112 (MQ=255)                                                                                     gCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCagcagc                                                                                                   <  4:198568/112‑1 (MQ=255)                                                                                        aaTCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCg                                                           >  3:288875/1‑149 (MQ=255)                                                                                             aTGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCATAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTc                                                      <  6:305684/149‑1 (MQ=255)                                                                                             aTGGCGCGAGTGCAGAGGGCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCATAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTc                                                      <  6:299974/149‑1 (MQ=255)                                                                                               ggCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGCCGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGc                                                    >  7:236168/1‑149 (MQ=255)                                                                                                     aGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGtt                                              <  7:227846/149‑1 (MQ=255)                                                                                                      gTGCAGCGGTCAGCAGGAACAGGACCGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCTGATTGGTACCGGTAGTAATGGCGCCCGCCCGCCGCTCGTAGTGCTTCAACGTGATAGCGCCACTGTCCATTGCCGCTGTCGCTTCGtg                                              >  6:381532/1‑147 (MQ=255)                                                                                                       tGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGa                                            >  1:407799/1‑149 (MQ=255)                                                                                                         gagagTTCAGCAGGAACAGGAACGCATCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGata                                          <  3:251643/148‑1 (MQ=255)                                                                                                          agagGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGAc                                                                                                             <  2:148059/81‑1 (MQ=255)                                                                                                          agagGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGAc                                                                                                             >  1:148059/1‑81 (MQ=255)                                                                                                              gTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAAt                                     >  4:88473/1‑149 (MQ=255)                                                                                                              gTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAAt                                     >  7:88647/1‑149 (MQ=255)                                                                                                               tCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATc                                    >  7:210485/1‑149 (MQ=255)                                                                                                               tCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATc                                    >  1:401538/1‑149 (MQ=255)                                                                                                               tCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATc                                    >  8:283677/1‑149 (MQ=255)                                                                                                                 agcagGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGt                                  <  7:283677/149‑1 (MQ=255)                                                                                                                 agcagGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGt                                  >  4:184709/1‑149 (MQ=255)                                                                                                                  gcagGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTc                                 <  2:401538/149‑1 (MQ=255)                                                                                                                     ggaacaggaacGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTccc                               >  8:382258/1‑148 (MQ=255)                                                                                                                          aggaacGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGaa                         <  4:105270/149‑1 (MQ=255)                                                                                                                              acGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACg                     >  3:217995/1‑149 (MQ=255)                                                                                                                              acGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACg                     >  5:369715/1‑149 (MQ=255)                                                                                                                               cGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCCACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGt                    >  7:139105/1‑149 (MQ=255)                                                                                                                                     cTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATcc               >  3:86629/1‑148 (MQ=255)                                                                                                                                        aCCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGaaa           >  5:14614/1‑149 (MQ=255)                                                                                                                                            ttCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGa       >  1:30337/1‑149 (MQ=255)                                                                                                                                            ttCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGa       <  8:340620/149‑1 (MQ=255)                                                                                                                                             tCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGAt      >  5:91401/1‑149 (MQ=255)                                                                                                                                                 cGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGAca  <  3:184709/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACA  >  NC_000913/1233230‑1233522

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