Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,619,803	A→G	F120L (TTT→CTT)	uraA ←	uracil permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,619,803	0	A	G	98.4%	235.5 / ‑4.2	62	F120L (TTT→CTT)	uraA	uracil permease
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (38/23);  total (39/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

CCAGGGTGGTAATAGAGATGATGATGGTTTTGGAGTCTGGCGTTTGCCCTTCAGCCGGGAGTAAACCCGCCATACCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAACAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTCATACCCTAACGGCAACAGTAACAATACCGGTGAAATAAAGGCAAAGCTGGA  >  NC_000913/2619662‑2619950                                                                                                                                              |                                                                                                                                                    ccAGGGTGGTAATAGAGATGATGATGGTTTTGGAGTCTGGCGTTTGCCCTTCAGCCGGGAGTAAACCCGCCATACCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTc                                                                                                                                            >  2:973605/1‑151 (MQ=255)     ggtggtAATAGAGATGATGATGGTTTTGGAGTCTGGCGTTTGCCCTTCAGCCGGGAGTAAACCCGCCATACCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTccagc                                                                                                                                        >  1:898666/1‑151 (MQ=255)       tggtAATAGAGATGATGATGGTTTTGGAGTCTGGCGTTTGCCCTTCAGCCGGGAGTAAACCCGCCATACCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTccagcc                                                                                                                                       >  2:461396/1‑150 (MQ=255)       tggtAATAGAGATGATGATGGTTTTGGAGTCTGGCGTTTGCCCTTCAGCCGGGAGTAAACCCGCCATACCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGAAGCACGTccagcc                                                                                                                                       >  1:595258/1‑150 (MQ=255)         gtAATAGAGATGATGATGGTTTTGGAGTCTGGCGTTTGCCCTTCAGCCGGGAGTAAACCCGCCATACCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCcagccagc                                                                                                                                    >  1:901462/1‑151 (MQ=255)           aaTAGAGATGATGATGGTTTTGGAGTCTGGCGTTTGCCCTTCAGCCGGGAGTAAACCCGCCATACCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCg                                                                                                                                  >  2:184991/1‑151 (MQ=255)            aTAGAGATGATGATGGTTTTGGAGTCTGGCGTTTGCCCTTCAGCCGGGAGTAAACCCGCCATACCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAACAGCACGTCCAGCCAGCCgg                                                                                                                                 >  2:91700/1‑151 (MQ=255)              agagATGATGATGGTTTTGGAGTCTGGCGTTTGCCCTTCAGCCGGGAGTAAACCCGCCATACCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTc                                                                                                                               >  2:536360/1‑151 (MQ=255)                agaTGATGATGGTTTTGGAGTCTGGCGTTTGCCCTTCAGCCGGGAGTAAACCCGCCATACCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTccc                                                                                                                             >  1:508067/1‑151 (MQ=255)                    gatgatgGTTTTGGAGTCTGGCGTTTGCCCTTCAGCCGGGAGTAAACCCGCCATACCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCt                                                                                                                         >  1:1009180/1‑151 (MQ=255)                       gatgGTTTTGGAGTCTGGCGTTTGCCCTTCAGCCGGGAGTAAACCCGCCATACCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCttt                                                                                                                       >  1:28154/1‑150 (MQ=255)                               tGGAGTCTGGCGTTTGCCCTTCAGCCGGGAGTAAACCCGCCATACCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGAt                                                                                                              >  2:932043/1‑151 (MQ=255)                                     cTGGCGTTTGCCCTTCAGCCGGGAGTAAACCCGCCATACCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATaaa                                                                                                           >  2:446458/1‑148 (MQ=255)                                         cGTTTGCCCTTCAGCCGGGAGTAAACCCGCCATACCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAAcc                                                                                                    >  1:42747/1‑151 (MQ=255)                                          gTTTGCCCTTCAGCCGGGAGTAAACCCGCCATACCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAAcc                                                                                                    >  1:866001/1‑150 (MQ=255)                                           tttGCCCTTCAGCCGGGAGTAAACCCGCCATACCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAAcc                                                                                                    >  1:1005432/1‑149 (MQ=255)                                           tttGCCCTTCAGCCGGGAGTAAACCCGCCATACCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAAcc                                                                                                    >  1:563576/1‑149 (MQ=255)                                                 cTTCAGCCGGGAGTAAACCCGCCATACCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGa                                                                                             >  2:1125253/1‑150 (MQ=255)                                                 cTTCAGCCGGGAGTAAACCCGCCATACCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGa                                                                                             >  1:433074/1‑150 (MQ=255)                                                      gCCGGGAGTAAACCCGCCATACCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGc                                                                                        >  1:393855/1‑150 (MQ=255)                                                       gggagaagacccccGCCATACCGGCGGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGaa                                                                                                       <  1:219398/123‑1 (MQ=255)                                                       ccGGGAGTAAACCCGCCATACCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGCTAAAAGAAACCAGGCAGAACAGCAc                                                                                      >  2:890680/1‑151 (MQ=255)                                                         gggAGTAAACCCGCCATACCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCAcgc                                                                                    >  1:1047460/1‑151 (MQ=255)                                                            aGTAAACCCGCCATACCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCAcgccgc                                                                                 >  1:1060928/1‑151 (MQ=255)                                                                cgcccGCCATCCCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCAcgcc                                                                                   <  2:403251/143‑1 (MQ=255)                                                                 ccccGCCATCCCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACata                                                                            <  1:373408/150‑1 (MQ=255)                                                                           ccGCCAGCTCCCCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACataata                                                                         <  1:270562/144‑1 (MQ=255)                                                                            cGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAgccgcc                                                                 >  1:607124/1‑151 (MQ=255)                                                                             ggCAGCTACGCTCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACCTCCAGCCAGCCGGTCCCCGCTTTCTTCAAGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAgccgaac                                                                >  1:859736/1‑148 (MQ=255)                                                                                 gCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAgcg                                                            >  2:1125003/1‑151 (MQ=255)                                                                                       cccgccAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACt                                                       <  1:703016/150‑1 (MQ=255)                                                                                       cccgccAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACt                                                       <  1:887714/150‑1 (MQ=255)                                                                                        ccgccAGCTCCAGCCCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGa                                                         <  1:640288/147‑1 (MQ=255)                                                                                        ccgccAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGcgcg                                                          <  1:646511/146‑1 (MQ=255)                                                                                        ccgccAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTc                                                     <  1:460902/151‑1 (MQ=255)                                                                                        ccgccAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTc                                                     <  1:922710/151‑1 (MQ=255)                                                                                         cgccAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTc                                                     >  1:302429/1‑150 (MQ=255)                                                                                          gccAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTCAt                                                   <  1:915685/151‑1 (MQ=255)                                                                                             aGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTCATAcc                                                <  2:459760/151‑1 (MQ=255)                                                                                                tCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTTCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGCCTTCATACCCt                                              >  1:1159073/1‑150 (MQ=255)                                                                                                  cAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTCATACCCTAAc                                           >  2:249876/1‑151 (MQ=255)                                                                                                      ccGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTCATACCCTAACGGCa                                       <  1:244206/151‑1 (MQ=255)                                                                                                        gATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACTAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTCATACCCTAACGGCAAc                                     >  2:412449/1‑151 (MQ=255)                                                                                                           gACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTCATACCCTAACGGCAACAGt                                  >  2:560571/1‑151 (MQ=255)                                                                                                            acatcAACGATTGCGCCCATTGCCGCTGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAGCCAGGCAGAACAGCACGCCGCGCATAATAAAGCCGCCCAGCGCGACTTCATACCCTAACGGCAACAGt                                  <  2:550338/146‑1 (MQ=255)                                                                                                               gCAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTCATACCCTAACGGCAACAGTAAc                               >  2:595795/1‑150 (MQ=255)                                                                                                                cAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTCATACCCTAACGGCAACAGTAACaa                             <  1:973025/151‑1 (MQ=255)                                                                                                                cAACGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTCATACCCTAACGGCAACAGTAACa                              >  2:887236/1‑150 (MQ=255)                                                                                                                   cGATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTCATACCCTAACGGCAACAGTAACAATAc                          >  2:1259523/1‑151 (MQ=255)                                                                                                                    gATTGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTCATACCCTAACGGCAACAGTAACAATAcc                         <  2:482286/151‑1 (MQ=255)                                                                                                                       tGCGCCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTCATACCCTAACGGCAACAGTAACAATACCgg                       >  1:1063443/1‑150 (MQ=255)                                                                                                                        gcgcCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTCATACCCTAACGGCAACAGTAACAATACCGGTg                     <  1:486401/151‑1 (MQ=255)                                                                                                                         cgcCCATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTCATACCCTAACGGCAACAGTAACAATACCGGTGa                    <  1:402083/151‑1 (MQ=255)                                                                                                                           cccATTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTCATACCCTAACGGCAACAGTAACAATACCGGTGaaa                  <  1:793037/151‑1 (MQ=255)                                                                                                                              aTTGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTCATACCCTAACGGCAACAGTAACAATACCGGTGAAATa                >  2:916209/1‑150 (MQ=255)                                                                                                                                tGCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTCATACCCTAACGGCAACAGTAACAATACCGGTGAAATaaa              <  1:492434/150‑1 (MQ=255)                                                                                                                                 gCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTCATACCCTAACGGCAACAGTAACAATACCGGTGAAATAAAgg            >  2:909801/1‑151 (MQ=255)                                                                                                                                 gCCGCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTCATACCCTAACGGCAACAGTAACAATACCGGTGAAATAAAg             <  2:535585/150‑1 (MQ=255)                                                                                                                                    gCAGGTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTCATACCCTAACGGCAACAGTAACAATACCGGTGAAATAAAGGCaa         <  1:13616/151‑1 (MQ=255)                                                                                                                                        gTGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTCATACCCTAACGGCAACAGTAACAATACCGGTGAAATAAAGGCAAAGCt     >  2:70495/1‑151 (MQ=255)                                                                                                                                         tGGAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTCATACCCTAACGGCAACAGTAACAATACCGGTGAAATAAAGGCAAAGCTg    >  2:625947/1‑151 (MQ=255)                                                                                                                                           gAAGCAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTCATACCCTAACGGCAACAGTAACAATACCGGTGAAATAAAGGCAAAGCTGGa  <  1:569946/151‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                    CCAGGGTGGTAATAGAGATGATGATGGTTTTGGAGTCTGGCGTTTGCCCTTCAGCCGGGAGTAAACCCGCCATACCGGCAGCTACGCCCGCCAGCTCCAGACCGATGACGGCAACGATTGCGCCCATTGCCGCAGGTGGAAACAGCACGTCCAGCCAGCCGGTCCCCGCTTTCTTCACGATAAAAGAAACCAGGCAGAACAGCACGCCGCACATAATAAAGCCGCCCAGCGCGACTTCATACCCTAACGGCAACAGTAACAATACCGGTGAAATAAAGGCAAAGCTGGA  >  NC_000913/2619662‑2619950

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