Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,639,617	T→C	G147G (GGA→GGG)	rrrQ ←	Qin prophage; putative lysozyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,639,617	0	T	C	100.0%	218.9 / NA	61	G147G (GGA→GGG)	rrrQ	Qin prophage; putative lysozyme
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (37/24);  total (37/24)

ATTTCCGACGATAAATCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGTCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTTGTGGTTCG  >  NC_000913/1639471‑1639764                                                                                                                                                   |                                                                                                                                                    aTTTCCGACGATAAATCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGc                                                                                                                                                 >  1:637232/1‑151 (MQ=255)      cgacgaAAAACCCTAACAACACGAGAAAAATCCCCATGAATATCTGATTCCCTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCa                                                                                                                                               <  2:600764/148‑1 (MQ=255)        acgaTAAATCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCtt                                                                                                                                          >  1:786213/1‑151 (MQ=255)            tAAATCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATc                                                                                                                                      >  2:429024/1‑151 (MQ=255)                 ccTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATccacc                                                                                                                                  >  1:120521/1‑150 (MQ=255)                 ccTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATccacc                                                                                                                                  >  1:283042/1‑150 (MQ=255)                    gACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGc                                                                                                                               <  1:379964/150‑1 (MQ=255)                    gACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCg                                                                                                                              <  1:6558/151‑1 (MQ=255)                     acaacaCGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGa                                                                                                                             <  2:62156/151‑1 (MQ=255)                         cacGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCg                                                                                                                         >  2:233156/1‑151 (MQ=255)                           cGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCt                                                                                                                       >  2:279974/1‑151 (MQ=255)                                aaaaTCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCa                                                                                                                  <  1:191393/151‑1 (MQ=255)                                 aaaTCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCAt                                                                                                                 <  1:938216/151‑1 (MQ=255)                                  aaTCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATg                                                                                                                >  2:19862/1‑151 (MQ=255)                                  aaTCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATg                                                                                                                >  2:573809/1‑151 (MQ=255)                                     cACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCa                                                                                                              <  1:22728/150‑1 (MQ=255)                                       ccATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCt                                                                                                           <  2:995836/151‑1 (MQ=255)                                       ccATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCt                                                                                                           >  2:765907/1‑151 (MQ=255)                                           gAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTAc                                                                                                       >  1:282108/1‑151 (MQ=255)                                              tatCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGAt                                                                                                    >  1:972902/1‑151 (MQ=255)                                               atCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCACTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGAt                                                                                                    >  2:486508/1‑150 (MQ=255)                                                  tGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCAcc                                                                                                >  2:1025873/1‑151 (MQ=255)                                                  tGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCAc                                                                                                 >  1:412967/1‑150 (MQ=255)                                                       cACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCa                                                                                            <  2:596038/150‑1 (MQ=255)                                                        aCTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGc                                                                                             <  2:547160/148‑1 (MQ=255)                                                             tCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGcc                                                                                     >  1:267306/1‑151 (MQ=255)                                                                aTCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGcc                                                                                     <  2:48607/148‑1 (MQ=255)                                                                      cagcagGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATaa                                                                             >  2:1224704/1‑150 (MQ=255)                                                                         cagGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACg                                                                         >  2:489859/1‑151 (MQ=255)                                                                         cagGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACg                                                                         >  1:48060/1‑151 (MQ=255)                                                                              tAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGAc                                                                    >  2:1012516/1‑151 (MQ=255)                                                                               aaTGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGAc                                                                    <  2:1174316/150‑1 (MQ=255)                                                                                  gcgcTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACggg                                                                 >  2:433337/1‑150 (MQ=255)                                                                                    gcTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGaa                                                               >  1:657423/1‑150 (MQ=255)                                                                                      tctcCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAac                                                             >  2:786729/1‑150 (MQ=255)                                                                                         ccTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCCCCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACtt                                                         <  1:284303/151‑1 (MQ=255)                                                                                         ccTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAaca                                                            <  1:1025867/148‑1 (MQ=255)                                                                                         ccTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACtt                                                         <  2:514410/151‑1 (MQ=255)                                                                                         ccTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACt                                                          <  2:1248637/150‑1 (MQ=255)                                                                                          cTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACtt                                                         >  2:621545/1‑150 (MQ=255)                                                                                            ggTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTAcc                                                      >  1:583490/1‑151 (MQ=255)                                                                                               cACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCggg                                                   >  2:788910/1‑151 (MQ=255)                                                                                                acgacgAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCgggg                                                  >  1:769457/1‑151 (MQ=255)                                                                                                 cgacgaATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCgggg                                                  <  2:786207/150‑1 (MQ=255)                                                                                                     gaATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCaa                                              <  1:1218538/150‑1 (MQ=255)                                                                                                       aTAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGt                                           >  2:359923/1‑151 (MQ=255)                                                                                                       aTAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGt                                           >  2:748614/1‑151 (MQ=255)                                                                                                              gACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAggg                                    <  1:1105343/151‑1 (MQ=255)                                                                                                                  gTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACaa                                <  1:233142/151‑1 (MQ=255)                                                                                                                       aGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGa                           <  2:430035/151‑1 (MQ=255)                                                                                                                          tATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGAcgc                        >  1:733982/1‑151 (MQ=255)                                                                                                                            tttGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGAcgcg                       >  2:458461/1‑150 (MQ=255)                                                                                                                            tttGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGAcgcg                       >  1:50621/1‑150 (MQ=255)                                                                                                                              tGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGa                      >  2:1074208/1‑149 (MQ=255)                                                                                                                                aaCGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATAcc                  >  2:1114840/1‑151 (MQ=255)                                                                                                                                 aCGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATAccc                 >  2:579587/1‑151 (MQ=255)                                                                                                                                    aaTGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCt              >  2:504657/1‑151 (MQ=255)                                                                                                                                         ggCAATCGGGCCCGCCATCCTTAATCCACCAGCGAATCGCTGCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACTCTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCtt             <  2:966540/147‑1 (MQ=255)                                                                                                                                          gCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTtgt        >  1:29515/1‑151 (MQ=255)                                                                                                                                           cAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTtgtg       <  2:552439/151‑1 (MQ=255)                                                                                                                                           cAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTtgtg       <  2:120522/151‑1 (MQ=255)                                                                                                                                                gcgcCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTTGTGGTTCg  >  1:167186/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                    ATTTCCGACGATAAATCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGTCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTTGTGGTTCG  >  NC_000913/1639471‑1639764

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