Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,496,777	T→C	intergenic (+146/‑79)	ydcJ → / → opgD	putative metalloenzyme/OPG biosynthetic periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,496,777	0	T	C	100.0%	263.5 / NA	69	intergenic (+146/‑79)	ydcJ/opgD	putative metalloenzyme/OPG biosynthetic periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/38);  total (31/38)

AATCTGACCATCCGCCTTTGCAAAAATTTGCCTGATTTTACAAACGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGTCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGGCCGCCGTGTGCGGTACCAGCGGCATTG  >  NC_000913/1496632‑1496922                                                                                                                                                  |                                                                                                                                                  aaTCTGACCATCCGCCTTTGCAAAAATTTGCCTGATTTTACAAACGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCataa                                                                                                                                              <  1:162292/151‑1 (MQ=255)     tGACCATCCGCCTTTGCAAAAATTTGCCTGATTTTACAAACGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCataaata                                                                                                                                           >  2:12694/1‑150 (MQ=255)     tGACCATCCGCCTTTGCAAAAATTTGCCTGATTTTACAAACGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCataaata                                                                                                                                           >  1:380236/1‑150 (MQ=255)      gACCATCCGCCTTTGCAAAAATTTGCCTGATTTTACAAACGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCataaa                                                                                                                                             <  1:196906/147‑1 (MQ=255)      gACCATCCGCCTTTGCAAAAATTTGCCTGATTTTACAAACGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAg                                                                                                                                         <  1:574017/151‑1 (MQ=255)        ccATCCGCCTTTGCAAAAATTTGCCTGATTTTACAAACGAATCAGGCTCATCCCATCGACAT‑AAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCAc                                                                                                                                      >  1:683486/1‑151 (MQ=255)          aTCCGCCTTTGCAAAAATTTGCCTGATTTTACAAACGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACt                                                                                                                                     >  1:1190720/1‑151 (MQ=255)            ccGCCTTTGCAAAAATTTGCCTGATTTTACAAACGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACttt                                                                                                                                   >  1:389008/1‑151 (MQ=255)               ccTTTGCAAAAATTTGCCTGATTTTACAAACGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCtt                                                                                                                                <  2:840186/151‑1 (MQ=255)                cTTTGCAAAAATTTGCCTGATTTTACAAACGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAGCAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTg                                                                                                                               >  2:952898/1‑151 (MQ=255)                     cAAAAATTTGCCTGATTTTACAAACGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCt                                                                                                                          <  1:779337/151‑1 (MQ=255)                     cAAAAATTTGCCTGATTTTACAAACGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCt                                                                                                                          <  2:930351/151‑1 (MQ=255)                     cAAAAATTTGCCTGATTTTACAAACGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCt                                                                                                                          <  2:963818/151‑1 (MQ=255)                      aaaaaTTTGCCTGATTTTACAAACGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCt                                                                                                                          >  2:824944/1‑150 (MQ=255)                      aaaaaTTTGCCTGATTTTACAAACGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTg                                                                                                                         >  2:1145034/1‑151 (MQ=255)                        aaaTTTGCCTGATTTTACAAACGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGaa                                                                                                                       <  1:290804/151‑1 (MQ=255)                        aaaTTTGCCTGATTTTACAAACGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGa                                                                                                                        >  1:337156/1‑150 (MQ=255)                          aTTTGCCTGATTTTACAAACGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGaaaa                                                                                                                     <  2:494354/151‑1 (MQ=255)                           tttGCCTGATTTTACAAACGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAAc                                                                                                                    >  1:795881/1‑151 (MQ=255)                             tGCCTGATTTTACAAACGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACg                                                                                                                   <  1:938995/150‑1 (MQ=255)                                  gATTTTACAAACGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCACAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGGTGAAAACGACCAgc                                                                                                             <  2:1105961/151‑1 (MQ=255)                                    ttttACAAACGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAgcgc                                                                                                           >  1:6693/1‑151 (MQ=255)                                       tACAAACGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTTAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCgg                                                                                                         <  2:909706/150‑1 (MQ=255)                                         cAAACGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGAc                                                                                                      >  1:518701/1‑151 (MQ=255)                                            aCGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCAt                                                                                                   <  1:178399/151‑1 (MQ=255)                                              gAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTc                                                                                                 <  2:863634/151‑1 (MQ=255)                                                  cAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAATGACCAGCGCGGGACCATTCacaa                                                                                             <  1:613110/151‑1 (MQ=255)                                                   aGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCacaac                                                                                            >  2:887852/1‑151 (MQ=255)                                                        cATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAg                                                                                       <  1:396645/151‑1 (MQ=255)                                                         aTCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACAcc                                                                                         >  1:607627/1‑148 (MQ=255)                                                                      aaaaaaaaTGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACttt                                                                          >  2:1242905/1‑150 (MQ=255)                                                                      aaaaaaaaTGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTc                                                                         >  2:784438/1‑151 (MQ=255)                                                                        aaaaaaTGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTc                                                                         >  1:1012012/1‑149 (MQ=255)                                                                         aaaaaTGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAgg                                                                      >  1:524491/1‑151 (MQ=255)                                                                            aaTGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTAt                                                                   >  1:1014770/1‑151 (MQ=255)                                                                             aTGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATg                                                                  <  1:540637/151‑1 (MQ=255)                                                                               gCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGa                                                                <  2:387004/151‑1 (MQ=255)                                                                                 cGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATc                                                              <  2:215989/151‑1 (MQ=255)                                                                                     ttATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGt                                                            >  1:1208583/1‑149 (MQ=255)                                                                                       aTGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGaa                                                        <  2:6800/151‑2 (MQ=255)                                                                                       aTGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGAc                                                        >  1:811920/1‑151 (MQ=255)                                                                                       aTGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGAc                                                        <  1:470869/151‑1 (MQ=255)                                                                                          cATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGAt                                                     <  1:70052/151‑1 (MQ=255)                                                                                          cATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGAt                                                     >  2:253553/1‑151 (MQ=255)                                                                                                  ctcAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATttatt                                                <  1:390984/148‑1 (MQ=255)                                                                                                    cAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTaaa                                             <  1:525308/149‑1 (MQ=255)                                                                                                    cAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTaa                                              <  1:784999/148‑1 (MQ=255)                                                                                                     aGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTAAAGGt                                          >  1:850150/1‑151 (MQ=255)                                                                                                          aacaaTTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGGATGGATCGTAGACGATTTATTAAAGGTTCaa                                      <  1:552644/150‑1 (MQ=255)                                                                                                            caaTTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATgg                                   >  1:1075040/1‑151 (MQ=255)                                                                                                             aaTTGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGc                                  >  2:1117654/1‑151 (MQ=255)                                                                                                                tGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTa                                <  1:767409/150‑1 (MQ=255)                                                                                                                tGGATTATTAATAAATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTAt                               >  1:476266/1‑151 (MQ=255)                                                                                                                        taataaATATTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGGccgccg                       >  1:1049859/1‑151 (MQ=255)                                                                                                                         aataaataTTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGGccgccg                       <  2:474623/150‑1 (MQ=255)                                                                                                                          ataaataTTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGGCCGCCgtg                     <  1:422051/151‑1 (MQ=255)                                                                                                                            aaataTTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGGCCGCCGtgtg                   <  2:761885/151‑1 (MQ=255)                                                                                                                            aaataTTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGGCCGCCGtgtg                   <  1:952901/151‑1 (MQ=255)                                                                                                                             aataTTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGGCCGCCGTGTGc                  <  2:52137/151‑1 (MQ=255)                                                                                                                               tatTGTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGGCCGCCGTGTGCgg                >  2:916859/1‑151 (MQ=255)                                                                                                                                  tGTTTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGCTTTATTAAAGGATCAATAGCTATGGCCGCCGAGTGCGGTAc             >  2:858299/1‑151 (MQ=255)                                                                                                                                   gTCTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGGCCGCCGTGTGCGGTAcc            <  2:388993/151‑1 (MQ=255)                                                                                                                                     cTAGAGTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGGCCGCCGTGTGCGGTAcc            >  2:601875/1‑149 (MQ=255)                                                                                                                                       agagTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGGCCGCCGTGTGCGGTACCa           <  2:951611/148‑1 (MQ=255)                                                                                                                                       agagTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGGCCGCCGTGTGCGGTACCAGc         <  1:1242891/150‑1 (MQ=255)                                                                                                                                       agagTGAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGGCCGCCGTGTGCGGTACCAGCg        <  2:869459/151‑1 (MQ=255)                                                                                                                                            gAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGGCCGCCGTGTGCGGTACCAGCGGCAtt   >  1:846494/1‑151 (MQ=255)                                                                                                                                            gAGCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGGCCGCCGTGTGCGGTACCAGCGGCAt    <  1:669161/150‑1 (MQ=255)                                                                                                                                              gCGGCCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGGCCGCCGTGTGCGGTACCAGCGGCATTg  <  1:784444/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  AATCTGACCATCCGCCTTTGCAAAAATTTGCCTGATTTTACAAACGAATCAGGCTCATCCCATCGACATAAAAAAAATGCCGATTTATGCATATTCTCTCAGTTCAACAATTGGATTATTAATAAATATTGTCTAGAGTGAGCGGTCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGGCCGCCGTGTGCGGTACCAGCGGCATTG  >  NC_000913/1496632‑1496922

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