Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,811,560	G→A	6.8%	A209A (GCG→GCA)	ydjN →	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,811,560	0	G	A	6.8%	182.4 / 2.9	59	A209A (GCG→GCA)	ydjN	putative transporter
Reads supporting (aligned to +/- strand):  ref base G (22/33);  new base A (2/2);  total (24/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

TCTGAGCGTTCCGCAGCTGGTCTTGTCCTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTGCAGGTTCTAACTTGCAAGACATCATC  >  NC_000913/1811413‑1811707                                                                                                                                                    |                                                                                                                                                    tctgaGCGTTCCGCAGCTGGTCTTGTCCTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATgcgc                                                                                                                                                    >  5:117035/1‑149 (MQ=255)     aGCGTTCCGCAGCTGGTCTTGTCCTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCTTAGCTGCATTGAAACTGCTGAAGGATGATGCGcc                                                                                                                                                   <  1:187862/146‑1 (MQ=255)                    gTCTTGTCCTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCt                                                                                                                                  >  7:124791/1‑148 (MQ=255)                    gTCTTGTCCTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCt                                                                                                                                  <  3:49716/148‑1 (MQ=255)                     tCTTGTCCTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTa                                                                                                                                <  6:117035/149‑1 (MQ=255)                       ttGTCCTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCACCGAAAGGTGAACGCGTCTTAGc                                                                                                                              >  6:133955/1‑149 (MQ=255)                       ttGTCCTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCACCGAAAGGTGAACGCGTCTTAGc                                                                                                                              >  6:137213/1‑149 (MQ=255)                         gTCCTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAg                                                                                                                               <  8:16245/146‑1 (MQ=255)                         gTCCTTTATCCCGAAAAACCCGTTTGACGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAACCTGATGAAGCATGATGCGCCGAAAGGTGAACGCCTATTAGCCg                                                                                                                            >  5:116011/1‑149 (MQ=255)                          tCCTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGc                                                                                                                           >  8:169457/1‑149 (MQ=255)                            cTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTa                                                                                                                         <  8:118535/149‑1 (MQ=255)                            cTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTa                                                                                                                         >  4:58038/1‑149 (MQ=255)                             tttATCCCGAAAACCCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTAGCCGCATTCCTCGGCGTATCTGCATTGAAACTGCTGAAGTATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTAt                                                                                                                        <  7:130511/149‑1 (MQ=255)                                 tCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGa                                                                                                                     >  3:52324/1‑148 (MQ=255)                                   ccGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATa                                                                                                                   <  4:194170/148‑1 (MQ=255)                                   ccGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATAc                                                                                                                  >  5:201826/1‑149 (MQ=255)                                   ccGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATAc                                                                                                                  >  5:112319/1‑149 (MQ=255)                                   ccGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATAc                                                                                                                  <  4:36976/149‑1 (MQ=255)                                   ccGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATAc                                                                                                                  >  3:36976/1‑149 (MQ=255)                                    cGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATAcc                                                                                                                 >  5:64897/1‑149 (MQ=255)                                               tttGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTg                                                                                                      <  4:213768/149‑1 (MQ=255)                                                 ttcagATCTTACCTGAGCCTATGCGACGTCAATTATCAGCGTTGTAATTTTTGCCGCATTCCTCGGCGTAGCTGGATTGAAACTTCTGAAGGATGATGCGCCGAAAGTTGAACGCGTCTTAGCCGCTATCTATACCCTAGAAAGCTgg                                                                                                     <  4:22289/144‑1 (MQ=255)                                                   ccGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGt                                                                                                   <  7:169457/148‑1 (MQ=255)                                                           aCCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGtgat                                                                                                <  1:100267/143‑1 (MQ=255)                                                           aCCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGtga                                                                                                 <  2:37267/142‑1 (MQ=255)                                                           aCCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGaaa                                                                                            <  8:134907/147‑1 (MQ=255)                                                             cGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTgg                                                                                        >  4:126748/1‑149 (MQ=255)                                                                aGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTc                                                                                     <  3:58038/149‑1 (MQ=255)                                                                    aaTCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCt                                                                                 >  7:37394/1‑149 (MQ=255)                                                                      tCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTgg                                                                               >  6:179148/1‑149 (MQ=255)                                                                        cGACGTCAAATATCAGCGGGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAATGGTGAACGCGTCTTAGCCGCTATAGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGGCTGGTc                                                                             <  6:195726/149‑1 (MQ=255)                                                                         gACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCa                                                                            <  2:35592/149‑1 (MQ=255)                                                                               aaTTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGt                                                                      >  6:176181/1‑149 (MQ=255)                                                                                aTTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGtt                                                                     <  7:213846/149‑1 (MQ=255)                                                                                   aTCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGa                                                                   <  5:234141/148‑1 (MQ=255)                                                                                             tAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGATCCCTTACGGc                                                        >  7:228833/1‑149 (MQ=255)                                                                                             tAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGc                                                        >  1:51551/1‑149 (MQ=255)                                                                                                tttttGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGtt                                                     <  6:64897/149‑1 (MQ=255)                                                                                                      ccGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGct                                               <  5:176181/149‑1 (MQ=255)                                                                                                       cGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGctc                                              <  5:165600/149‑1 (MQ=255)                                                                                                          aTTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTaa                                           <  6:112319/149‑1 (MQ=255)                                                                                                           ttCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAAt                                          <  7:232243/149‑1 (MQ=255)                                                                                                              cTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGa                                        <  8:124791/148‑1 (MQ=255)                                                                                                              cTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGGTAACGCGTCTTAGCCGCGATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGTTCTAATGAc                                       <  5:35768/149‑1 (MQ=255)                                                                                                                           cATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTg                            <  5:179148/147‑1 (MQ=255)                                                                                                                           cATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTGc                           <  4:52324/148‑1 (MQ=255)                                                                                                                            aTTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTGCAg                         >  2:33213/1‑149 (MQ=255)                                                                                                                            aTTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCAGGCTCTAATGACCAAAGTCGTTGCAg                         >  2:91703/1‑149 (MQ=255)                                                                                                                              tGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTaa                                           >  4:194648/1‑129 (MQ=255)                                                                                                                              tGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTaa                                           <  3:194648/129‑1 (MQ=255)                                                                                                                                       tgAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTGCAGGTTCTAACTTg              >  6:177838/1‑149 (MQ=255)                                                                                                                                         aaGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGtt                             <  7:43441/132‑1 (MQ=255)                                                                                                                                         aaGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGtt                             >  8:43441/1‑132 (MQ=255)                                                                                                                                          aGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGATCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTGCAGGTTCTAACTTGCaa           <  8:228833/149‑1 (MQ=255)                                                                                                                                          aGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTGCAGGTTCTAACTTGCaa           <  6:201826/149‑1 (MQ=255)                                                                                                                                          aGGATGATGCACCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTGCAGGTTCTAACTTGCaa           <  5:133955/149‑1 (MQ=255)                                                                                                                                          aGGATGATGCACCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTGCAGGTTCTAACTTGCaa           <  5:137213/149‑1 (MQ=255)                                                                                                                                            gatgatGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTGCAGGTTCTAACTTGCAAGa         >  7:4687/1‑149 (MQ=255)                                                                                                                                                 tgCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTg                            <  6:2924/125‑1 (MQ=255)                                                                                                                                                 tgCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTg                            >  5:2924/1‑125 (MQ=255)                                                                                                                                                    gcCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTGCAGGTTCTAACTTGCAAGACatcatc  >  4:102116/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TCTGAGCGTTCCGCAGCTGGTCTTGTCCTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTGCAGGTTCTAACTTGCAAGACATCATC  >  NC_000913/1811413‑1811707

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