Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,673,783	G→A	7.2%	intergenic (+15/+33)	yphA → / ← yphB	DoxX family inner membrane protein/mutarotase superfamily protein, YphB family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,673,783	0	G	A	7.2%	152.7 / 2.9	56	intergenic (+15/+33)	yphA/yphB	DoxX family inner membrane protein/mutarotase superfamily protein, YphB family
Reads supporting (aligned to +/- strand):  ref base G (22/29);  new base A (2/2);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCATTACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCATCGGTTC  >  NC_000913/2673634‑2673932                                                                                                                                                       |                                                                                                                                                      cACCATTACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGTTg                                                                                                                                                       >  2:1262333/1‑150 (MQ=255) cACCATTACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATg                                                                                                                                                       <  1:307567/150‑1 (MQ=255)  aCCATTACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgc                                                                                                                                                      >  2:50880/1‑150 (MQ=255)  aCCATTACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgc                                                                                                                                                      >  2:612293/1‑150 (MQ=255)         cTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGa                                                                                                                                               <  2:445969/150‑1 (MQ=255)           gggATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAAc                                                                                                                                             >  2:264950/1‑150 (MQ=255)            gtatatGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTAATGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACg                                                                                                                                            <  2:76562/148‑1 (MQ=255)              atatGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGcc                                                                                                                                          <  2:659803/150‑1 (MQ=255)              atatGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGcc                                                                                                                                          >  2:442222/1‑150 (MQ=255)               tatGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCt                                                                                                                                         <  2:1037942/150‑1 (MQ=255)               tatGACCGGCGATGCGGGAGCGCCAAATATGCTTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCAGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGcc                                                                                                                                          <  2:1019577/150‑1 (MQ=255)                 tGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCATTTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                       >  2:792506/1‑150 (MQ=255)                 tGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTAGCTATTGGCAAATACCGCGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGTTGCGGCGTGAACGCCtga                                                                                                                                       >  1:634396/1‑148 (MQ=255)                    ccGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                    >  2:343185/1‑150 (MQ=255)                    ccGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                    <  1:1301391/150‑1 (MQ=255)                       gCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCCGTCTCCCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                 <  1:521847/150‑1 (MQ=255)                         gATGCGGTTGGGCCAAATATGATTAATATCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTCGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                                                                               <  1:783316/150‑1 (MQ=255)                           tGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGCGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATTCGGCGTGAACGCCTTATCCGGCAttt                                                                                                                             >  2:29518/1‑150 (MQ=255)                              ggTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaat                                                                                                                          <  1:251287/150‑1 (MQ=255)                               gTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaata                                                                                                                         >  2:1092404/1‑150 (MQ=255)                                   ggCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAAtt                                                                                                                     >  2:762281/1‑150 (MQ=255)                                      cAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACa                                                                                                                  <  2:724819/150‑1 (MQ=255)                                        aaTATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGc                                                                                                                <  2:359344/150‑1 (MQ=255)                                          tatGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTTATAAATTACAGCCa                                                                                                              >  1:124325/1‑150 (MQ=255)                                           atGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCAc                                                                                                             >  1:483067/1‑150 (MQ=255)                                            tGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCAct                                                                                                            >  2:464976/1‑150 (MQ=255)                                            tGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCAct                                                                                                            <  1:442222/150‑1 (MQ=255)                                              aTTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGACAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCActca                                                                                                          >  2:567633/1‑149 (MQ=255)                                                 aaTTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAcc                                                                                                       <  2:5218/150‑1 (MQ=255)                                                   tttCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCg                                                                                                     <  2:901224/150‑1 (MQ=255)                                                    ttCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGc                                                                                                    <  1:922699/150‑1 (MQ=255)                                                    ttCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGc                                                                                                    <  2:563315/150‑1 (MQ=255)                                                        ggAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACg                                                                                                >  2:1216730/1‑150 (MQ=255)                                                          aagaaTGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGAc                                                                                              >  2:1059690/1‑150 (MQ=255)                                                               tGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGGC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAAATACAGCCACTCTACCCGCAACGACAtctc                                                                                         <  2:885644/150‑1 (MQ=255)                                                                     tATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgt                                                                                   <  1:521926/150‑1 (MQ=255)                                                                             gcgcGTTCTTGCTATTGGCAATTACCGGGCC‑GCGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTGCCAGCAACGACATCTCTGAAGTTGTAGagt                                                                           >  1:97146/1‑148 (MQ=255)                                                                              cgcgTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTc                                                                          <  2:897470/150‑1 (MQ=255)                                                                              cgcgTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGACACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTc                                                                          >  1:676614/1‑150 (MQ=255)                                                                                  ttCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCt                                                                      >  1:345707/1‑150 (MQ=255)                                                                                  tgcTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGGGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCt                                                                      <  2:105180/148‑1 (MQ=255)                                                                                          aTTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGCCATCTCTGAAGTTGTTGCTTCCCCTGGTGCTAg                                                              >  2:756352/1‑150 (MQ=255)                                                                                                     ccgggcc‑gggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGAt                                                   <  1:1059301/150‑1 (MQ=255)                                                                                                        ggcc‑gggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCAc                                                >  1:426550/1‑150 (MQ=255)                                                                                                              ggggCAATTTCTCTCGATCGGCGTTAGGAAAGATTCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCt                                           >  1:1095832/1‑150 (MQ=255)                                                                                                                  cAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCg                                       <  1:464976/150‑1 (MQ=255)                                                                                                                     tttCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGc                                    <  2:1045556/150‑1 (MQ=255)                                                                                                                          ctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCgatga                               >  2:759723/1‑150 (MQ=255)                                                                                                                          ctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGCTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCgatga                               <  1:50880/150‑1 (MQ=255)                                                                                                                             gATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCgatgatga                            >  1:287612/1‑150 (MQ=255)                                                                                                                                cGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGatca                         <  1:882634/150‑1 (MQ=255)                                                                                                                                cGGAGTTAGTAAAGATGCCGGATGCGGGGTGAACGACTTATCCGGCATTTAATAAATTACAGCCTCTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGAGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGatca                         <  1:211943/150‑1 (MQ=255)                                                                                                                                  gCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCTGCATTTAATAAAATACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGatcatc                       <  2:1301957/150‑1 (MQ=255)                                                                                                                                       aGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGgcg                  <  2:352136/150‑1 (MQ=255)                                                                                                                                        ggAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGgcgc                 >  1:913915/1‑150 (MQ=255)                                                                                                                                                 ccGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCGTCCGGCGCATGGCTCAt        >  2:181923/1‑150 (MQ=255)                                                                                                                                                  cGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACATCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCATc       <  2:267842/150‑1 (MQ=255)                                                                                                                                                  cGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCATc       <  1:612293/150‑1 (MQ=255)                                                                                                                                                  cGGATAAGGCGTTCACGCCGCATCCGGCATTTAATTAATTACAGCCACTCTACCCGCAACGACATCTCTGATGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCATc       >  1:1299409/1‑150 (MQ=255)                                                                                                                                                    gATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCATCgg     <  1:249598/150‑1 (MQ=255)                                                                                                                                                       gcggcgTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCATCGGTTc  >  1:848622/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                      CACCATTACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCATCGGTTC  >  NC_000913/2673634‑2673932

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