Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,614,957	A→C	S431R (AGT→CGT)	nikA →	nickel/heme ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,614,957	0	A	C	100.0%	338.7 / NA	85	S431R (AGT→CGT)	nikA	nickel/heme ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (43/42);  total (43/42)

CATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTGACCCGTC  >  NC_000913/3614809‑3615104                                                                                                                                                     |                                                                                                                                                    cATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGt                                                                                                                                                   >  1:1108152/1‑151 (MQ=255)   ttCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTc                                                                                                                                                 >  1:735383/1‑151 (MQ=255)    tCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTc                                                                                                                                                 <  1:1193756/150‑1 (MQ=255)    tCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTc                                                                                                                                                 >  2:1213944/1‑150 (MQ=255)      aGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAAt                                                                                                                                              <  1:14533/151‑1 (MQ=255)       ggCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATg                                                                                                                                             <  2:577861/151‑1 (MQ=255)         cTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATgcg                                                                                                                                           >  1:546678/1‑151 (MQ=255)            atatGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGcgcg                                                                                                                                         >  1:719428/1‑150 (MQ=255)                  gcCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTc                                                                                                                                  <  2:1066149/151‑1 (MQ=255)                   ccAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTca                                                                                                                                 <  2:1099967/151‑1 (MQ=255)                    cAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTcac                                                                                                                                <  2:795545/151‑1 (MQ=255)                          ggCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCGCCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCt                                                                                                                           <  2:217501/150‑1 (MQ=255)                            cgcAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGAc                                                                                                                        >  1:682490/1‑151 (MQ=255)                               aGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTc                                                                                                                     >  2:1274933/1‑151 (MQ=255)                                 aTGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTcc                                                                                                                    >  2:506521/1‑150 (MQ=255)                                  tGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTcc                                                                                                                    >  2:639108/1‑149 (MQ=255)                                  tGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTc                                                                                                                     >  1:701939/1‑148 (MQ=255)                                      tcGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGcac                                                                                                              >  2:1167066/1‑151 (MQ=255)                                       cGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGcac                                                                                                              <  1:1284498/150‑1 (MQ=255)                                       cGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGca                                                                                                               <  2:939931/149‑1 (MQ=255)                                       cGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGc                                                                                                                <  2:480512/148‑1 (MQ=255)                                        gCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAg                                                                                                            <  2:120036/151‑1 (MQ=255)                                         cTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAg                                                                                                            <  1:757523/150‑1 (MQ=255)                                         cTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGc                                                                                                           >  2:905421/1‑151 (MQ=255)                                         cTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGc                                                                                                           >  1:660024/1‑151 (MQ=255)                                           gATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCaa                                                                                                         <  1:987196/151‑1 (MQ=255)                                            aTTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCa                                                                                                          >  2:292947/1‑149 (MQ=255)                                             ttGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAgg                                                                                                       >  1:647488/1‑151 (MQ=255)                                               ggCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGa                                                                                                      <  2:522314/150‑1 (MQ=255)                                                  gaagaaGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAg                                                                                                  <  1:1214558/151‑1 (MQ=255)                                                    agaagaGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGcc                                                                                                >  1:1230947/1‑151 (MQ=255)                                                    agaagaGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGcc                                                                                                >  2:668133/1‑151 (MQ=255)                                                      aagaGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCg                                                                                               <  2:485238/150‑1 (MQ=255)                                                       agagagCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCGCCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGAc                                                                                             >  2:767564/1‑151 (MQ=255)                                                       agagagCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGAc                                                                                             <  2:1067174/151‑1 (MQ=255)                                                       agagagCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGAc                                                                                             <  1:927587/151‑1 (MQ=255)                                                        gagagCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGAc                                                                                             >  2:1046025/1‑150 (MQ=255)                                                          gagCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTTCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACaa                                                                                           >  1:692097/1‑150 (MQ=255)                                                            gCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAAcc                                                                                        <  2:276052/151‑1 (MQ=255)                                                            gCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAAcc                                                                                        <  2:123887/151‑1 (MQ=255)                                                             cAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCg                                                                                       <  1:387282/151‑1 (MQ=255)                                                                  tCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCtgat                                                                                  >  2:657501/1‑151 (MQ=255)                                                                   cTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCtgatt                                                                                 >  1:293167/1‑151 (MQ=255)                                                                    tATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCtgattg                                                                                >  2:1269331/1‑151 (MQ=255)                                                                     aTGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCtgat                                                                                  <  2:1020088/148‑1 (MQ=255)                                                                       gCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCtgattgat                                                                              <  2:1199557/150‑1 (MQ=255)                                                                          cGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATaaa                                                                           <  2:781500/150‑1 (MQ=255)                                                                          cGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATaa                                                                            >  2:696889/1‑149 (MQ=255)                                                                             cAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAagag                                                                        <  1:1274022/150‑1 (MQ=255)                                                                             cAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAgaga                                                                       <  2:24355/151‑1 (MQ=255)                                                                              aGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGAt                                                                      <  2:560184/151‑1 (MQ=255)                                                                                cgcgACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCg                                                                    <  2:380334/151‑1 (MQ=255)                                                                                 gcgACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCgg                                                                   <  2:33809/151‑1 (MQ=255)                                                                                  cgACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGc                                                                  >  1:565227/1‑151 (MQ=255)                                                                                   gACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCg                                                                 >  2:387960/1‑151 (MQ=255)                                                                                      ggTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGa                                                                >  2:846818/1‑149 (MQ=255)                                                                                      ggTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAg                                                              <  1:854510/151‑1 (MQ=255)                                                                                         cGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGc                                                           >  1:847910/1‑151 (MQ=255)                                                                                              tGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGc                                                       <  2:546693/150‑1 (MQ=255)                                                                                                 cATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGAcc                                                   <  2:601066/151‑1 (MQ=255)                                                                                                  aTGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGActc                                                  >  1:409276/1‑149 (MQ=255)                                                                                                       tttCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGAc                                             >  1:1010671/1‑151 (MQ=255)                                                                                                        ttCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACg                                            >  1:335161/1‑151 (MQ=255)                                                                                                         tCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACg                                            >  1:540625/1‑150 (MQ=255)                                                                                                           cACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGaa                                          >  1:326267/1‑150 (MQ=255)                                                                                                            aCCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAAc                                        <  1:779935/151‑1 (MQ=255)                                                                                                                 aCCTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCaa                                    >  2:637399/1‑150 (MQ=255)                                                                                                                  ccTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACg                                  >  1:494788/1‑151 (MQ=255)                                                                                                                  ccTGGGGCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACg                                  >  2:469680/1‑151 (MQ=255)                                                                                                                     ggggCGCGCCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTc                                <  2:992096/150‑1 (MQ=255)                                                                                                                         cgcgcCATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGgcg                           >  2:528521/1‑151 (MQ=255)                                                                                                                            gcCATATGATCCACACGCCTTCCGCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGGAACGTAAGGCGCTg                        <  1:366127/151‑1 (MQ=255)                                                                                                                              cATATGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTa                      <  2:157674/151‑1 (MQ=255)                                                                                                                                 atGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGCTTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATc                    <  1:767535/150‑1 (MQ=255)                                                                                                                                 atGATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATc                    <  1:1087294/150‑1 (MQ=255)                                                                                                                                   gATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACTCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATcgcg                 <  1:599995/151‑1 (MQ=255)                                                                                                                                   gATCCACACGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATcgcg                 >  2:463941/1‑151 (MQ=255)                                                                                                                                        acacGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACAtt            <  1:668110/151‑1 (MQ=255)                                                                                                                                        acacGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACAtt            >  1:471959/1‑151 (MQ=255)                                                                                                                                         cacGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTc           <  2:176426/151‑1 (MQ=255)                                                                                                                                         cacGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTc           >  1:956780/1‑151 (MQ=255)                                                                                                                                           cGCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTg         <  1:632469/151‑1 (MQ=255)                                                                                                                                            gCCTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTg         >  1:1111846/1‑150 (MQ=255)                                                                                                                                             ccTTCCTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTgcc       >  1:157391/1‑149 (MQ=255)                                                                                                                                                tcctccGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTGACCCg    >  2:1133904/1‑151 (MQ=255)                                                                                                                                                  cTCCGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTGACCCGTc  >  2:339594/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTGACCCGTC  >  NC_000913/3614809‑3615104

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