Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,911,366	+T	coding (160/1041 nt)	pstS ←	phosphate ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,911,364	1	.	T	100.0%	264.2 / NA	68	N54K (AAT→AAA)	pstS	phosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (39/29);  total (39/29)

CACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTA‑TTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCACTCATCGATAAGGTCGCGGCGACAACAGTTGCGACGGTGG  >  NC_000913/3911214‑3911509                                                                                                                                                        |                                                                                                                                                  cacGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTtta                                                                                                                                                    <  1:620762/151‑1 (MQ=255)  acGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTtta‑t                                                                                                                                                  >  2:184256/1‑151 (MQ=255)  acGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTtta‑t                                                                                                                                                  >  1:499176/1‑151 (MQ=255)   cgccgcCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTtta‑tt                                                                                                                                                 >  2:1043314/1‑151 (MQ=255)       gccAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTtta‑tttacc                                                                                                                                             >  1:1123555/1‑147 (MQ=255)        ccAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTtta‑tttaccg                                                                                                                                            <  1:1137200/151‑6 (MQ=255)         cAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCgg                                                                                                                                            >  1:101940/1‑151 (MQ=255)             cACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGtttc                                                                                                                                        <  1:626879/151‑1 (MQ=255)               cGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGtttctt                                                                                                                                      <  2:46176/151‑1 (MQ=255)               cGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGtttctt                                                                                                                                      <  1:882133/151‑1 (MQ=255)               cGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGtttct                                                                                                                                       <  1:1099358/150‑1 (MQ=255)                  tCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTttctttct                                                                                                                                   >  1:653205/1‑151 (MQ=255)                   cGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTg                                                                                                                                  >  1:946750/1‑151 (MQ=255)                      gAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGt                                                                                                                                >  2:535418/1‑150 (MQ=255)                          tGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAgtg                                                                                                                           >  2:992640/1‑151 (MQ=255)                          tGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAgtg                                                                                                                           >  2:852897/1‑151 (MQ=255)                           ggAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAgtg                                                                                                                           <  2:1062437/150‑1 (MQ=255)                              aCAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGGAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTc                                                                                                                         >  2:627934/1‑149 (MQ=255)                                  aCCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGccc                                                                                                                    >  1:370784/1‑150 (MQ=255)                                  aCCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCa                                                                                                                   <  2:16018/151‑1 (MQ=255)                                   ccTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCAt                                                                                                                  <  2:954121/151‑1 (MQ=255)                                          gAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGc                                                                                                             <  1:57602/149‑1 (MQ=255)                                          gAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCAt                                                                                                           <  1:258647/151‑1 (MQ=255)                                           aGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCAt                                                                                                           >  1:346700/1‑150 (MQ=255)                                               aGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACAcc                                                                                                      >  1:205200/1‑151 (MQ=255)                                                  ttttCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCgg                                                                                                    <  2:75927/150‑1 (MQ=255)                                                     tCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGgcgc                                                                                                 <  1:705153/150‑1 (MQ=255)                                                        tCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAggg                                                                                             >  1:685222/1‑151 (MQ=255)                                                          aGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGa                                                                                            >  1:878743/1‑150 (MQ=255)                                                            aCAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAgg                                                                                         >  2:41469/1‑151 (MQ=255)                                                                  gcgcGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTtgcac                                                                                   >  2:484149/1‑151 (MQ=255)                                                                     cgTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTtgcacctg                                                                                >  2:169684/1‑151 (MQ=255)                                                                        cAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTtgcacctgcac                                                                             >  2:169591/1‑151 (MQ=255)                                                                         agagGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTtgcacctgcacc                                                                            <  1:111188/151‑1 (MQ=255)                                                                         agagGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTtgcacctgcacc                                                                            >  2:276319/1‑151 (MQ=255)                                                                          gagGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTgcacctgcacct                                                                           <  1:551208/151‑1 (MQ=255)                                                                             gCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTc                                                                        >  2:479584/1‑151 (MQ=255)                                                                              cACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTc                                                                        >  1:226484/1‑150 (MQ=255)                                                                               aCCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAg                                                                      <  1:340047/151‑1 (MQ=255)                                                                                  aaaaTCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGgct                                                                   <  2:387244/151‑1 (MQ=255)                                                                                  aaaaTCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGgct                                                                   >  2:843465/1‑151 (MQ=255)                                                                                   aaaTCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGgctt                                                                  >  1:115794/1‑151 (MQ=255)                                                                                                tAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAac                                                      >  1:462314/1‑150 (MQ=255)                                                                                                 aGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAacac                                                    <  2:81778/151‑1 (MQ=255)                                                                                                 aGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAacac                                                    >  2:960710/1‑151 (MQ=255)                                                                                                  gCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAacac                                                    >  1:560456/1‑150 (MQ=255)                                                                                                    gATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACaga                                                 <  2:749241/151‑1 (MQ=255)                                                                                                     aTAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGATTCTGCAAACACagag                                                <  2:19699/151‑1 (MQ=255)                                                                                                       aaTCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAgaga                                               >  1:120155/1‑150 (MQ=255)                                                                                                            gTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCAc                                         >  2:214815/1‑151 (MQ=255)                                                                                                               ttACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCACTc                                       >  1:869326/1‑150 (MQ=255)                                                                                                               ttACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCACTc                                       >  2:750867/1‑150 (MQ=255)                                                                                                                  cGCCAGAGGAAGAACCGATAAGCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTCGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCACTCATCg                                   >  1:157990/1‑151 (MQ=255)                                                                                                                  cGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCACTCATCg                                   <  2:1067371/151‑1 (MQ=255)                                                                                                                      aCCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCACTCATCGATaa                               <  2:249205/151‑1 (MQ=255)                                                                                                                      aCCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCACTCATCGATa                                >  1:887240/1‑150 (MQ=255)                                                                                                                      aCCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCACTCATCGATa                                >  2:203875/1‑150 (MQ=255)                                                                                                                      aCCGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCACTCATCGATa                                >  2:774036/1‑150 (MQ=255)                                                                                                                       ccGGAAGAACCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCGTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGTAAGGTGGCACCTGCACCTGTCAGGCTTGCTTTTGCAAACACAGAGAAAGCACTCATCGATaa                               <  1:749698/150‑1 (MQ=255)                                                                                                                          gaagaaCCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCACTCATCGATAAGGTc                           <  2:213216/151‑1 (MQ=255)                                                                                                                            agaaCCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCACTCATCGATAAGGTcgc                         >  1:809925/1‑151 (MQ=255)                                                                                                                              aaCCGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCACTCATCGATAAGGTCgcgg                       <  1:61212/151‑1 (MQ=255)                                                                                                                                ccGATACCCTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCACTCATCGATAAGGTCgcggcg                     <  1:191907/151‑1 (MQ=255)                                                                                                                                      cccTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCACTCATCGATAAGGTCGCGGCGacaaca               <  2:877163/151‑1 (MQ=255)                                                                                                                                      cccTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCACTCATCGATAAGGTCGCGGCGacaac                >  2:812299/1‑150 (MQ=255)                                                                                                                                       ccTGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCACTCATCGATAAGGTCGCGGCGacaac                >  1:295231/1‑149 (MQ=255)                                                                                                                                         tGGTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCACTCATCGATAAGGTCGCGGCGACAACAGtt            >  1:655836/1‑151 (MQ=255)                                                                                                                                           gTAGTTAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCACTCATCGATAAGGTCGCGGCGACAACAGTTGc          <  1:843327/151‑1 (MQ=255)                                                                                                                                                tAACTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCACTCATCGATAAGGTCGCGGCGACAACAGTTGCGACg      <  2:719142/150‑1 (MQ=255)                                                                                                                                                  aCTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCACTCATCGATAAGGTCGCGGCGACAACAGTTGCGACggtg   >  2:768220/1‑151 (MQ=255)                                                                                                                                                  aCTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCACTCATCGATAAGGTCGCGGCGACAACAGTTGCGACggtg   <  1:528373/151‑1 (MQ=255)                                                                                                                                                  aCTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCACTCATCGATAAGGTCGCGGCGACAACAGTTGCGACggtg   <  1:35648/151‑1 (MQ=255)                                                                                                                                                  aCTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCACTCATCGATAAGGTCGCGGCGACAACAGTTGCGACggtg   >  1:307122/1‑151 (MQ=255)                                                                                                                                                   cTTTATTTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCACTCATCGATAAGGTCGCGGCGACAACAGTTGCGACggtgg  >  2:13337/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                  CACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTA‑TTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTGCAAACACAGAGAAAGCACTCATCGATAAGGTCGCGGCGACAACAGTTGCGACGGTGG  >  NC_000913/3911214‑3911509

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