Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,513,835	T→C	E158E (GAA→GAG)	fecB ←	ferric citrate ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,513,835	0	T	C	100.0%	334.3 / NA	83	E158E (GAA→GAG)	fecB	ferric citrate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (43/40);  total (43/40)

ACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGG  >  NC_000913/4513689‑4513982                                                                                                                                                   |                                                                                                                                                    aCGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCTCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCa                                                                                                                                                 <  1:372163/151‑1 (MQ=255)   gTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTcc                                                                                                                                                  >  2:1074696/1‑148 (MQ=255)   gTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGa                                                                                                                                               <  2:546950/151‑1 (MQ=255)      cAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGt                                                                                                                                            >  2:1137306/1‑151 (MQ=255)         ccccAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGcc                                                                                                                                         >  2:1060580/1‑151 (MQ=255)            cAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGc                                                                                                                                      <  1:32153/151‑1 (MQ=255)            cAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGc                                                                                                                                      >  1:261407/1‑151 (MQ=255)             agagagGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGc                                                                                                                                      >  2:10352/1‑150 (MQ=255)             agagagGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGc                                                                                                                                      >  2:225402/1‑150 (MQ=255)                  ggCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCAtctct                                                                                                                                <  1:1321542/151‑1 (MQ=255)                      aGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCt                                                                                                                            <  1:398643/151‑1 (MQ=255)                       gCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCtt                                                                                                                           <  1:1273127/151‑1 (MQ=255)                       gCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCtt                                                                                                                           <  2:1158710/151‑1 (MQ=255)                         aCGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCtttt                                                                                                                         <  2:1237436/151‑1 (MQ=255)                         aCGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCtttt                                                                                                                         <  2:701835/151‑1 (MQ=255)                             tGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTAcc                                                                                                                     >  1:54619/1‑151 (MQ=255)                             tGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGAGCCTGCATCTCTCGCTTTTTAcc                                                                                                                     >  2:1076639/1‑151 (MQ=255)                              gCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTAc                                                                                                                      <  1:1125921/149‑1 (MQ=255)                               ccGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTAccc                                                                                                                    >  2:550946/1‑150 (MQ=255)                               ccGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTAccc                                                                                                                    >  2:1183845/1‑150 (MQ=255)                                 ggTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACccacc                                                                                                                 <  1:870653/151‑1 (MQ=255)                                 ggTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACccacc                                                                                                                 >  2:447150/1‑151 (MQ=255)                                 ggTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACccacc                                                                                                                 >  1:1237578/1‑151 (MQ=255)                                  gTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACccacc                                                                                                                 >  2:358528/1‑150 (MQ=255)                                    ccAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCAtt                                                                                                              >  2:450827/1‑151 (MQ=255)                                     cAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGGCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCAttt                                                                                                             >  2:520637/1‑151 (MQ=255)                                     cAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCAttt                                                                                                             >  1:1092305/1‑151 (MQ=255)                                      aGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTc                                                                                                            >  2:530342/1‑151 (MQ=255)                                       ggTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTc                                                                                                            <  2:1383964/150‑1 (MQ=255)                                        gTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGc                                                                                                          <  2:486626/151‑1 (MQ=255)                                            ccTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCga                                                                                                       <  2:1282757/150‑1 (MQ=255)                                            ccTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCga                                                                                                       <  2:254381/150‑1 (MQ=255)                                             cTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCgatg                                                                                                     >  1:137123/1‑151 (MQ=255)                                              tGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCgatg                                                                                                     >  1:417247/1‑150 (MQ=255)                                              tGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCgatg                                                                                                     >  2:36944/1‑150 (MQ=255)                                                     gCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCg                                                                                              <  1:210141/150‑1 (MQ=255)                                                       aGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAg                                                                                           <  2:586193/151‑1 (MQ=255)                                                       aGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAg                                                                                           <  1:463895/151‑1 (MQ=255)                                                        ggTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGa                                                                                          <  1:12441/151‑1 (MQ=255)                                                         gTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGAt                                                                                         <  2:388294/151‑1 (MQ=255)                                                             aaTTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCa                                                                                     <  1:1296440/151‑1 (MQ=255)                                                                 gCTGTTCGCGTGATGCGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAAtt                                                                                 <  1:762339/151‑1 (MQ=255)                                                                    gTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTc                                                                              <  2:277145/151‑1 (MQ=255)                                                                    gTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTc                                                                              >  1:1175236/1‑151 (MQ=255)                                                                      tCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCa                                                                             <  1:276088/150‑1 (MQ=255)                                                                       cgcgTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGc                                                                           >  1:83295/1‑151 (MQ=255)                                                                        gcgTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCg                                                                          <  1:34316/151‑1 (MQ=255)                                                                             aTGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGt                                                                     >  1:80118/1‑151 (MQ=255)                                                                               gtgCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGttt                                                                   <  2:1205738/151‑1 (MQ=255)                                                                                  ccAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGt                                                                <  2:1462196/151‑1 (MQ=255)                                                                                    aaaGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTg                                                              <  1:1082198/151‑1 (MQ=255)                                                                                     aaGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGtt                                                               <  2:757022/149‑1 (MQ=255)                                                                                          cACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGAc                                                        <  1:582941/151‑1 (MQ=255)                                                                                            cgcgTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACtt                                                      <  2:54612/151‑1 (MQ=255)                                                                                                tgtCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAagc                                                  >  1:1205460/1‑151 (MQ=255)                                                                                                 gtCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAagc                                                  >  1:1364665/1‑150 (MQ=255)                                                                                                    ccTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAgcagca                                              <  2:781982/151‑1 (MQ=255)                                                                                                      tttGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGt                                            >  1:977491/1‑151 (MQ=255)                                                                                                      tttGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGt                                            >  2:494100/1‑151 (MQ=255)                                                                                                       ttGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGcagcag                                             <  2:807342/149‑1 (MQ=255)                                                                                                          ggAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCg                                        <  2:428621/151‑1 (MQ=255)                                                                                                           gAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCgg                                       >  1:87364/1‑151 (MQ=255)                                                                                                                   ctggcCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGAttt                               >  2:1272873/1‑151 (MQ=255)                                                                                                                      gcCCACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTtgct                            <  1:520695/151‑1 (MQ=255)                                                                                                                        ccACTGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTtgctgc                          >  2:468133/1‑151 (MQ=255)                                                                                                                            tGCGCCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAgg                      >  2:471427/1‑151 (MQ=255)                                                                                                                             gcgcCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAgg                      <  1:1074797/150‑1 (MQ=255)                                                                                                                             gcgcCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAg                       <  2:1047512/149‑1 (MQ=255)                                                                                                                              cgcCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCg                    >  2:1296408/1‑151 (MQ=255)                                                                                                                              cgcCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCg                    >  1:708643/1‑151 (MQ=255)                                                                                                                              cgcCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCg                    >  2:226057/1‑151 (MQ=255)                                                                                                                               gcCATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCg                    >  1:548810/1‑150 (MQ=255)                                                                                                                                ccATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGAt                  <  1:351233/151‑1 (MQ=255)                                                                                                                                 cATCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATg                 >  1:970341/1‑151 (MQ=255)                                                                                                                                   tCCTCTCTTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGt                >  2:1062229/1‑150 (MQ=255)                                                                                                                                     ctctctTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTaa              >  1:998043/1‑150 (MQ=255)                                                                                                                                     ctctctTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTaa              >  1:236688/1‑150 (MQ=255)                                                                                                                                         ctTTATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAAcccc         >  1:999255/1‑151 (MQ=255)                                                                                                                                           ttATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAAccccc        <  2:1465294/150‑1 (MQ=255)                                                                                                                                            tATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGc      >  2:518260/1‑151 (MQ=255)                                                                                                                                            tATGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGc      >  2:457658/1‑151 (MQ=255)                                                                                                                                             aTGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGc      >  1:575637/1‑150 (MQ=255)                                                                                                                                             aTGTTGCTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCa     <  1:225424/151‑1 (MQ=255)                                                                                                                                                ttgctCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATgg  >  1:884915/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                    ACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGG  >  NC_000913/4513689‑4513982

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