Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,186,605	A→C	H419P (CAC→CCC)	rpoC →	RNA polymerase, beta prime subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,186,605	0	A	C	100.0%	229.9 / NA	60	H419P (CAC→CCC)	rpoC	RNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (28/32);  total (28/32)

AGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACC  >  NC_000913/4186458‑4186742                                                                                                                                                    |                                                                                                                                          agaaAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACcccc                                                                                                                                        <  1:701355/151‑1 (MQ=255) agaaAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACcccc                                                                                                                                        >  2:48419/1‑151 (MQ=255)                aGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCg                                                                                                                         >  2:608536/1‑151 (MQ=255)                    gTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTg                                                                                                                       <  2:670995/149‑1 (MQ=255)                    gTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGc                                                                                                                      >  1:908761/1‑150 (MQ=255)                      tCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCAc                                                                                                                    <  2:931550/150‑1 (MQ=255)                           ccGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGActc                                                                                                              <  1:869553/151‑1 (MQ=255)                             gTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTg                                                                                                            <  1:887058/151‑1 (MQ=255)                                cATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCAc                                                                                                         >  2:647879/1‑151 (MQ=255)                                     aCGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTc                                                                                                     <  1:597585/150‑1 (MQ=255)                                     aCGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCt                                                                                                    <  1:474105/151‑1 (MQ=255)                                      cGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTg                                                                                                   <  1:175988/151‑1 (MQ=255)                                        gCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTggg                                                                                                 >  2:271124/1‑151 (MQ=255)                                          aaGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGt                                                                                                >  2:594561/1‑150 (MQ=255)                                              tGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATcc                                                                                            >  1:832014/1‑150 (MQ=255)                                               ggAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATcc                                                                                            <  2:624877/149‑1 (MQ=255)                                               ggAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCa                                                                                           <  2:276490/150‑1 (MQ=255)                                                gAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAGGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAgg                                                                                         >  2:102875/1‑151 (MQ=255)                                                gAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAgg                                                                                         <  1:802278/151‑1 (MQ=255)                                                  aCTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGc                                                                                        >  2:886180/1‑150 (MQ=255)                                                     gCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAttt                                                                                    >  1:282571/1‑151 (MQ=255)                                                      cGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGt                                                                                                      >  1:514628/1‑132 (MQ=255)                                                      cGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGt                                                                                                      <  2:514639/132‑1 (MQ=255)                                                         ggTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGaa                                                                                 <  1:611174/150‑1 (MQ=255)                                                               gCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTa                                                                           <  1:71554/150‑1 (MQ=255)                                                                 tACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTg                                                                        >  1:162827/1‑151 (MQ=255)                                                                  accaccATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACt                                                                         <  2:111187/149‑1 (MQ=255)                                                                  accaccATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTg                                                                        <  1:61761/150‑1 (MQ=255)                                                                    caccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATc                                                                     <  1:467998/151‑1 (MQ=255)                                                                       caTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGaa                                                                  <  2:29010/151‑1 (MQ=255)                                                                             aGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTaa                                                             >  2:147281/1‑150 (MQ=255)                                                                               cTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGc                                                          >  1:631443/1‑151 (MQ=255)                                                                                  cGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTAt                                                       <  1:608521/151‑1 (MQ=255)                                                                                   gaagaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATc                                                      <  2:682651/151‑1 (MQ=255)                                                                                        aaaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCt                                                 >  1:260593/1‑151 (MQ=255)                                                                                        aaaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCt                                                 <  1:267043/151‑1 (MQ=255)                                                                                        aaaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCt                                                 >  2:325391/1‑151 (MQ=255)                                                                                        aaaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCt                                                 >  2:315951/1‑151 (MQ=255)                                                                                        aaaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCt                                                 >  1:639702/1‑151 (MQ=255)                                                                                         aaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCt                                                 >  1:937515/1‑150 (MQ=255)                                                                                         aaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTg                                                >  1:236808/1‑151 (MQ=255)                                                                                            ggTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCAc                                             >  2:246447/1‑151 (MQ=255)                                                                                              ttGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTAGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCAccc                                           >  2:380770/1‑151 (MQ=255)                                                                                               tGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCAccc                                           <  1:688860/150‑1 (MQ=255)                                                                                                  gcgcGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTg                                       <  2:908776/151‑1 (MQ=255)                                                                                                   cgcgAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTgg                                      <  2:91028/151‑1 (MQ=255)                                                                                                     cgAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGtt                                    <  2:380453/151‑1 (MQ=255)                                                                                                            gCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCgg                             <  1:152939/151‑1 (MQ=255)                                                                                                                 cGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCatat                        >  1:770068/1‑151 (MQ=255)                                                                                                                  gTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCatat                        <  1:936883/150‑1 (MQ=255)                                                                                                                     tGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAAc                     <  1:78490/150‑1 (MQ=255)                                                                                                                        gATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCg                 <  2:16668/151‑1 (MQ=255)                                                                                                                         atatCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGa                <  2:918056/151‑1 (MQ=255)                                                                                                                             ccTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTc            <  1:48421/151‑1 (MQ=255)                                                                                                                             ccTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTc            >  2:332825/1‑151 (MQ=255)                                                                                                                               tGGACGAAGTTATCCGCGAACCCCCGGTGCTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTc            <  1:385606/149‑1 (MQ=255)                                                                                                                               tGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCg           >  2:361075/1‑150 (MQ=255)                                                                                                                                   cGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGt      <  2:10776/151‑1 (MQ=255)                                                                                                                                      aGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGAc   >  1:963791/1‑151 (MQ=255)                                                                                                                                      aGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGAc   >  2:503702/1‑151 (MQ=255)                                                                                                                                        ttATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGAcc  >  2:510040/1‑150 (MQ=255)                                                                                                                                         tATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGAcc  >  2:726417/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                          AGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACC  >  NC_000913/4186458‑4186742

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