Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,609,163	T→G	M173R (ATG→AGG)	yjjG →	dUMP phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,609,163	0	T	G	98.6%	267.4 / NA	70	M173R (ATG→AGG)	yjjG	dUMP phosphatase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (40/29);  total (41/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAG  >  NC_000913/4609018‑4609305                                                                                                                                                  |                                                                                                                                               cAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTg                                                                                                                                           <  1:26655/151‑1 (MQ=255)   ggTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGc                                                                                                                                         <  2:221520/151‑1 (MQ=255)   ggTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGc                                                                                                                                         <  1:44976/151‑1 (MQ=255)           tGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGc                                                                                                                                 >  2:417151/1‑151 (MQ=255)            ggAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGcc                                                                                                                                >  2:97816/1‑151 (MQ=255)            ggAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGcc                                                                                                                                <  1:270683/151‑1 (MQ=255)            ggAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGc                                                                                                                                 <  1:1002884/150‑1 (MQ=255)             gAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGcc                                                                                                                                >  2:85420/1‑150 (MQ=255)              aaCGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAACCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCg                                                                                                                               >  1:910144/1‑150 (MQ=255)                       gCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGatat                                                                                                                     <  2:1094635/151‑1 (MQ=255)                        ccTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATAtt                                                                                                                    >  1:806006/1‑151 (MQ=255)                          tGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGAGCTGAGGGTTGGCGACACTGCCGAGTCCGATATtc                                                                                                                   >  2:203802/1‑150 (MQ=255)                           gCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATtct                                                                                                                  <  2:983627/150‑1 (MQ=255)                           gCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATtct                                                                                                                  <  1:9735/150‑1 (MQ=255)                                aTTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTgg                                                                                                            >  2:1106610/1‑151 (MQ=255)                                   aCTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCa                                                                                                          <  2:236905/150‑1 (MQ=255)                                    cTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATc                                                                                                        >  1:678822/1‑151 (MQ=255)                                      tCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATc                                                                                                        <  2:570277/149‑1 (MQ=255)                                      tCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATc                                                                                                        <  1:112842/149‑1 (MQ=255)                                              cTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGc                                                                                              >  2:1073105/1‑151 (MQ=255)                                                ggTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCtt                                                                                            <  2:228206/151‑1 (MQ=255)                                                  tGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCCATACTGAACGGTCAAGCGTGCTGATGGTTGGGCACACTGCCGAGGCCGATATTCTCGGTGGCATCAACGCCGGGCTTg                                                                                           >  2:6911/1‑150 (MQ=255)                                                   gATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTACCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGc                                                                                          <  2:552073/150‑1 (MQ=255)                                                       tCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGAcc                                                                                      >  2:947297/1‑150 (MQ=255)                                                        ccGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACct                                                                                     <  2:34261/150‑1 (MQ=255)                                                        ccGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACct                                                                                     <  2:804830/150‑1 (MQ=255)                                                         cGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACctgc                                                                                   >  2:771216/1‑151 (MQ=255)                                                          gaagaaGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACctgct                                                                                  <  2:799454/151‑1 (MQ=255)                                                          gaagaaGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACctgct                                                                                  >  1:1036258/1‑151 (MQ=255)                                                           aagaagTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACctgctg                                                                                 >  1:1073480/1‑151 (MQ=255)                                                             gaagTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTgctggc                                                                               >  1:947357/1‑151 (MQ=255)                                                                gTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTgctggctg                                                                             >  2:423426/1‑150 (MQ=255)                                                                   ggCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATg                                                                         >  1:420400/1‑151 (MQ=255)                                                                      gTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGcac                                                                      >  1:342315/1‑151 (MQ=255)                                                                        tGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGcaca                                                                     <  1:982143/150‑1 (MQ=255)                                                                        tGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCacac                                                                    >  1:701692/1‑151 (MQ=255)                                                                         gCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACAcc                                                                   >  2:65210/1‑151 (MQ=255)                                                                         gCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACAcc                                                                   >  2:615133/1‑151 (MQ=255)                                                                           cAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGGTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCAt                                                                 >  1:190320/1‑151 (MQ=255)                                                                            aaaCCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTCCTGGCTGAATGCACACCATc                                                                <  2:116327/151‑1 (MQ=255)                                                                              aCCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATcgc                                                              >  1:374331/1‑151 (MQ=255)                                                                               ccGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATcgcg                                                             >  2:210046/1‑151 (MQ=255)                                                                                cGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATcgcg                                                             >  1:608806/1‑150 (MQ=255)                                                                                  aaTAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGc                                                          <  2:667195/151‑1 (MQ=255)                                                                                     aaGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAAc                                                       <  1:421151/151‑1 (MQ=255)                                                                                          aaTTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGa                                                   >  1:639417/1‑150 (MQ=255)                                                                                          aaTTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCCCGAGCAACCAg                                                    >  2:67226/1‑149 (MQ=255)                                                                                           aTTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGa                                                   >  1:764750/1‑149 (MQ=255)                                                                                              ttCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGc                                               >  2:122509/1‑150 (MQ=255)                                                                                              ttCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGc                                               <  2:936942/150‑1 (MQ=255)                                                                                              ttCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGc                                               <  2:243805/150‑1 (MQ=255)                                                                                                  aTTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATcgc                                          <  2:894876/151‑1 (MQ=255)                                                                                                   ttATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATcgcg                                         >  2:957938/1‑151 (MQ=255)                                                                                                    tATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCgcgc                                        >  2:877305/1‑151 (MQ=255)                                                                                                        cgcTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCAc                                    >  2:486857/1‑151 (MQ=255)                                                                                                         gcTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCAcc                                   >  1:741133/1‑151 (MQ=255)                                                                                                            ggAACAGGCGGGCAATCCTGACCGTTCACGCGTGGTGAGGGTTGGCGACGCTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTg                                 <  2:914247/150‑1 (MQ=255)                                                                                                            ggAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTgg                                >  1:1026996/1‑151 (MQ=255)                                                                                                             gAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTgg                                >  1:1089138/1‑150 (MQ=255)                                                                                                               aCAGGCGGGCAATCCTGTCCGTTCGCGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGa                               <  2:29583/149‑1 (MQ=255)                                                                                                                      ggCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTttctt                      >  2:649166/1‑151 (MQ=255)                                                                                                                       gCAATCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTttcttc                     <  2:411765/151‑1 (MQ=255)                                                                                                                           tCCTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTg                 >  1:370217/1‑151 (MQ=255)                                                                                                                            ccTGACCGTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTg                 <  1:730976/150‑1 (MQ=255)                                                                                                                                   gTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACt         <  1:310796/151‑1 (MQ=255)                                                                                                                                   gTTCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACt         >  1:209689/1‑151 (MQ=255)                                                                                                                                     tCACGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTgg       >  2:260009/1‑151 (MQ=255)                                                                                                                                       aCGCGTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTg        >  2:978901/1‑148 (MQ=255)                                                                                                                                         gcgTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGagc    >  2:1115403/1‑150 (MQ=255)                                                                                                                                          cgTGCTGAGGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGagcag  <  1:493144/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                               CAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAG  >  NC_000913/4609018‑4609305

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