Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,673,784	C→A	7.1%	intergenic (+16/+32)	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,784	0	C	A	7.1%	137.5 / 2.5	56	intergenic (+16/+32)	yphA/yphB	DoxX family inner membrane protein/mutarotase superfamily protein, YphB family
Reads supporting (aligned to +/- strand):  ref base C (23/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 = 6.73e-01

ACCATTACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCATCGGTTC  >  NC_000913/2673635‑2673932                                                                                                                                                       |                                                                                                                                                     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:442222/1‑150 (MQ=255)             atatGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGcc                                                                                                                                          <  2:659803/150‑1 (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                                                                                                                                    <  1:1301391/150‑1 (MQ=255)                   ccGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                    >  2:343185/1‑150 (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)                                                                                                                                                       |                                                                                                                                                     ACCATTACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCC‑GGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCATCGGTTC  >  NC_000913/2673635‑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: ‑