Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,362,581	C→A	37.3%	noncoding (46/76 nt)	pheU ←	tRNA‑Phe

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,362,581	0	C	A	37.3%	43.9 / 51.6	69	noncoding (46/76 nt)	pheU	tRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base C (25/18);  new base A (17/9);  total (42/27)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCTGTGT  >  NC_000913/4362436‑4362727                                                                                                                                                  |                                                                                                                                                   gTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACg                                                                                                                                                 <  1:250769/149‑1 (MQ=255)          ggATAATCTGGTTGAGTTCGGTAGCATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGACCTGAAACGAAAAACGGACCTCCGTGGAGGGCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTc                                                                                                                                        <  2:226381/149‑1 (MQ=255)          ggAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAAAGAAAAACGGACATCCGTGGAGGTCCTTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGTAAACTGGGAtttta                                                                                                                                        >  1:190438/1‑148 (MQ=255)           gagaATCTGGGTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGGCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCa                                                                                                                                       <  2:983579/149‑1 (MQ=255)            atattgtGGTTGAGTTCGGTAGTATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAATCGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGATCCAAGGAAACGGGGATTTTCaa                                                                                                                                      <  1:386188/143‑1 (MQ=39)            agaATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGATCCAAGGACACGGGGATTTTCaa                                                                                                                                      >  1:49624/1‑149 (MQ=255)              aaTCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATc                                                                                                                                    <  1:510541/149‑1 (MQ=255)               aTCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGTAATCGAACCAAGGACACGGGGATTTTCAATcc                                                                                                                                   >  1:955132/1‑149 (MQ=255)                tCTGGTTGAGTACGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATccc                                                                                                                                  >  1:751019/1‑149 (MQ=255)                   ggTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGACCGTTTATATGAATTTGGTGCCCGGACACGGAATCGACCCAAGGAAACGGAGCTTTTCAATCCCCTg                                                                                                                               >  1:933014/1‑149 (MQ=255)                     ttGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTTTATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGAAACGGGGATTTTCAATCCCCAGct                                                                                                                             >  2:65958/1‑149 (MQ=255)                     ttGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGct                                                                                                                             >  2:482426/1‑149 (MQ=255)                          ttCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGAAACGGGGATTTTCAATCCCCTGCTCTAcc                                                                                                                        >  2:961180/1‑149 (MQ=255)                           tCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCg                                                                                                                       >  2:1004231/1‑149 (MQ=255)                             ggTAGAATTGATTTGGAGGCAGAACGCTTAACTCGTGGCGACCCGAAACCAAACCCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACTGGGATTTTCAATCCCCTGCTCTACCGAc                                                                                                                     <  1:236399/149‑1 (MQ=37)                                    ttGAATTGGAGGCAGAACGCATAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGACCGTTTATAGGAATTTGGTGCCAGGACTCGGAATCCAACCAACTAAAACGGCATTTTCAATCCCCTGCTCAACCGAATGaccaa                                                                                                              >  2:953267/1‑145 (MQ=255)                                        tttGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACCGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTCCCGACTGAGCTATCCg                                                                                                          >  2:767968/1‑149 (MQ=255)                                              ggCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGAAACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAAc                                                                                                    >  1:123038/1‑149 (MQ=255)                                               gCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCACCTGCTCTACCGACTGAGCTATCAGGTCaagt                                                                                                   >  1:637652/1‑147 (MQ=255)                                                  gAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACggg                                                                                                 <  2:377836/148‑1 (MQ=255)                                                  gAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACggg                                                                                                 >  1:377836/1‑148 (MQ=255)                                                            aTCGTGGCGTCCTGAAACGAAAAACGGACCTCCGGGGAGGTCCGTTTATTTGAATTTGGTGCCCGGACTCGGAATCGACCCAAGGCCACGGGGCTTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAAc                                                                                      >  1:82145/1‑149 (MQ=21)                                                             tCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGAAACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAAcc                                                                                     >  2:805928/1‑149 (MQ=255)                                                                gggCGTCCGGAAACGACACACGGTCCTCCGGGGCGGGCCGATTATATGAATTTGGTGCCCGGACGCGGAAGCGAACCAAGGACACGGCGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTa                                                                                  <  2:714273/148‑1 (MQ=17)                                                                 ggCGTCCTGAAACGAAAAACGGACCTCCGTGTAGGTCCGTTTATATGAATTTGGTGCCCTGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTaa                                                                                 <  2:65177/149‑1 (MQ=255)                                                                     tCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTg                                                                             >  1:278247/1‑149 (MQ=255)                                                                      ccTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGAAACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGc                                                                            <  1:716854/149‑1 (MQ=255)                                                                        tGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGAAACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGcac                                                                          >  2:459650/1‑149 (MQ=255)                                                                          aaaCGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGCGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACAt                                                                        >  2:541795/1‑149 (MQ=255)                                                                              gAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCg                                                                    >  2:650937/1‑149 (MQ=255)                                                                              gAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGAAACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACAtggcg                                                                    >  1:415569/1‑145 (MQ=255)                                                                              gAAAAACGGACATCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTACAACGTAAGATGCACATCTCg                                                                    >  2:760224/1‑149 (MQ=18)                                                                                   aCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGAAACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAAc                                                               >  2:429468/1‑149 (MQ=255)                                                                                   aCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGAAACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAAc                                                               >  1:153883/1‑149 (MQ=255)                                                                                    cGGCCCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGAAACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAATCCGTAATCTGCACATCTCGTCAAcc                                                              <  2:851180/149‑1 (MQ=255)                                                                                    cGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCGAGTCAAcc                                                              >  2:615254/1‑149 (MQ=255)                                                                                         cTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAAtt                                                         >  1:368809/1‑149 (MQ=255)                                                                                         cTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGAAACGGGTATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGCGCGCATTAAACCGTAATCTGCACATATCGTCAACCTAAtt                                                         >  1:993439/1‑149 (MQ=255)                                                                                          tCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAAttt                                                        <  2:3394/149‑1 (MQ=255)                                                                                          tCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAAttt                                                        <  2:633629/149‑1 (MQ=255)                                                                                            cGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGAAACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCa                                                      <  1:961180/149‑1 (MQ=255)                                                                                               ggAGGTCCGTTTATATGAATTAGGTGCTCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGATCTATCCGGGCAACGGTTCGCATTAATTCTTAATCTGCACATCTCGTCAACCTAATTTCAGTa                                                   <  1:79151/149‑1 (MQ=18)                                                                                                gAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGAAACGGGGAATTTCAATCCCCTGCTCTACCGACTGAGCGATCCGGGCAACGGGGCGCATTAATCCGTAATCTGCCCCTCTCGTGCACCTAAATTCAGTaa                                                  >  2:176777/1‑149 (MQ=18)                                                                                                   gTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGAAACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCTTTAAACCGTAATCTCCACATCTCGTCAACCTAATTTCTGTAAAAg                                               <  2:366830/149‑1 (MQ=255)                                                                                                   gTCCGTTTATATGAATTTGGTGCCCGGACACGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGATCTATCCGGGGAACTGGGCGCATTAAACCGTAATCTGGACATCGCGGCATCCTAATTTCAGTAAAAg                                               >  2:563097/1‑149 (MQ=18)                                                                                                    tCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGc                                              <  2:279658/149‑1 (MQ=255)                                                                                                        tttATATGAATTTGGTGCCCGGACTCTGAATCGAACCAAGGACACGGGGATTTTAAATCCCCTGATCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTATACCGTAATCTGCACATCTAGTCAACCTAATTTCAGTAAAACCGAtt                                          >  1:197279/1‑149 (MQ=21)                                                                                                            tatGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGAAACGGGGATTTTCAATCCCCTGCTCTACCGACTTAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTATAAGCGATTCAAc                                      <  1:628083/149‑1 (MQ=255)                                                                                                            tatGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGAAACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTTACCCGTAATCTGCACATGTCGTCAACCTAAATTCAGTAAAAGCGATTAAAc                                      >  2:763258/1‑149 (MQ=18)                                                                                                             atGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCAGATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAAcg                                     >  1:292507/1‑148 (MQ=255)                                                                                                             atGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGAAACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCTTTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACt                                     >  1:1008515/1‑149 (MQ=255)                                                                                                              tGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCACTCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTg                                    >  2:537967/1‑149 (MQ=255)                                                                                                                aaTTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACAACTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTTCt                                  >  1:608922/1‑149 (MQ=255)                                                                                                                 aTTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCACCCTAATTTCAGTAAAAGCGATTCAACTGCtt                                 >  1:381471/1‑149 (MQ=255)                                                                                                                 aTTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCtt                                 >  2:542984/1‑149 (MQ=255)                                                                                                                  tttGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTa                                <  1:1004231/149‑1 (MQ=255)                                                                                                                      gTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACGGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGAt                            >  2:788389/1‑149 (MQ=255)                                                                                                                       tGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAAtt                                                         <  2:165408/119‑1 (MQ=255)                                                                                                                       tGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAAtt                                                         >  1:165408/1‑119 (MQ=255)                                                                                                                             gACTCGGAATCGAACCAAGGAAACGGGGATTTTCAATCCCCTGCTCTCCCGACTGAGCTAACCGGGCAACGGGGCGCATTACAACGTAATCTGCACATCTCGTCAACCTAATTTCAGAAAAAGCGATTCACCTGCATAAGATTGCGGCa                     >  2:213534/1‑149 (MQ=21)                                                                                                                                tCGTAATCGAAACTAGGTAACGGGGATTTTCAATCCACTGATCTACCTACTGAGCTATCCGGGCAAAGTGTCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTTCTTAAGATTGCGGCAAAt                  <  2:933014/149‑1 (MQ=255)                                                                                                                                 cGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATc                 <  2:137045/149‑1 (MQ=255)                                                                                                                                   gAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGc               <  2:278247/149‑1 (MQ=255)                                                                                                                                      tCGAACCAAGGACACGGGGACTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGcttc            <  1:445514/149‑1 (MQ=255)                                                                                                                                       cGAACCAAGGAAACGGGGATTTTCAATCCCCTGCTCTACCTACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGcttct           <  1:760224/149‑1 (MQ=255)                                                                                                                                         aaCCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCttt         >  2:592606/1‑149 (MQ=255)                                                                                                                                         aaCCAAGGAAACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCttt         >  1:759780/1‑149 (MQ=255)                                                                                                                                            cAAGGAAACTGGGATTTTCAATCCCCTGCTCTACCGACTGTGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGGCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCt      <  1:593581/149‑1 (MQ=255)                                                                                                                                                gACACGGGGATTTTCAATGCCCTGCTCTACCGACTGAGCTATCCGGGGAACGGGGCGCATTAAATCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCTgtgt  <  2:751019/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCTGTGT  >  NC_000913/4362436‑4362727

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