Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,294,366	A→T	intergenic (+227/‑115)	nrfG → / → gltP	heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfG/glutamate/aspartate:proton symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,294,366	0	A	T	100.0%	310.6 / NA	78	intergenic (+227/‑115)	nrfG/gltP	heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfG/glutamate/aspartate:proton symporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (43/35);  total (43/35)

CAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAATTCAGCCTGGCCTG  >  NC_000913/4294218‑4294509                                                                                                                                                     |                                                                                                                                                cAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTc                                                                                                                                               <  2:746701/151‑1 (MQ=255)  aGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTc                                                                                                                                               >  1:654785/1‑150 (MQ=255)   gTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCAt                                                                                                                                             >  1:392641/1‑151 (MQ=255)   gTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCAt                                                                                                                                             >  1:509523/1‑151 (MQ=255)     attattGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATa                                                                                                                                            >  1:771203/1‑150 (MQ=255)        attGCATATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATAATGTTCATAAg                                                                                                                                          <  2:279633/149‑1 (MQ=255)        attGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCg                                                                                                                                        <  2:978484/151‑1 (MQ=255)            cAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTa                                                                                                                                     >  2:296732/1‑150 (MQ=255)              gATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaa                                                                                                                                  <  2:903130/151‑1 (MQ=255)              gATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaa                                                                                                                                   >  1:791987/1‑150 (MQ=255)                  aaaCAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAAt                                                                                                                              <  2:572876/151‑1 (MQ=255)                       aaTAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTAc                                                                                                                          >  2:187145/1‑150 (MQ=255)                                  ttCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACg                                                                                                                 <  2:733897/148‑1 (MQ=255)                                  ttCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCa                                                                                                               >  1:885803/1‑150 (MQ=255)                                    ccATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAac                                                                                                             <  2:879494/150‑1 (MQ=255)                                     cATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAacac                                                                                                           >  1:863001/1‑151 (MQ=255)                                     cATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAacac                                                                                                           <  2:427781/151‑1 (MQ=255)                                      aTAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAacac                                                                                                           >  2:107043/1‑150 (MQ=255)                                      aTAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAacac                                                                                                           >  1:623839/1‑150 (MQ=255)                                        aaCAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAAt                                                                                                        >  2:815014/1‑151 (MQ=255)                                         aCAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAAtt                                                                                                       <  1:263380/151‑1 (MQ=255)                                            aaTGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCAt                                                                                                    <  1:481995/151‑1 (MQ=255)                                             aTGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATg                                                                                                   >  1:583245/1‑151 (MQ=255)                                                gTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGccc                                                                                                >  2:443101/1‑151 (MQ=255)                                                 ttATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCt                                                                                               >  1:830699/1‑151 (MQ=255)                                                 ttATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCt                                                                                               >  1:15370/1‑151 (MQ=255)                                                  tATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCt                                                                                               <  1:939367/150‑1 (MQ=255)                                                   aTTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTgg                                                                                             >  1:423261/1‑151 (MQ=255)                                                      cATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGc                                                                                            >  1:96479/1‑149 (MQ=255)                                                      cATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGc                                                                                            >  1:53470/1‑149 (MQ=255)                                                       aTTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGt                                                                                         <  2:19433/151‑1 (MQ=255)                                                         tAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTAt                                                                                       >  1:626827/1‑151 (MQ=255)                                                            tCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTc                                                                                    >  2:773811/1‑151 (MQ=255)                                                              cTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCAc                                                                                  >  1:904019/1‑151 (MQ=255)                                                               tGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACg                                                                                 >  2:933113/1‑151 (MQ=255)                                                                      aaCTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTcgc                                                                          >  2:382769/1‑151 (MQ=255)                                                                          cgcgTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGAGGAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTcg                                                                           >  1:111804/1‑146 (MQ=255)                                                                            cgTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCt                                                                    >  2:858228/1‑151 (MQ=255)                                                                             gTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTg                                                                   >  2:214318/1‑151 (MQ=255)                                                                             gTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTg                                                                   >  2:464366/1‑151 (MQ=255)                                                                              ttattaTGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTg                                                                   >  1:328319/1‑150 (MQ=255)                                                                              ttattaTGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTg                                                                   >  1:367351/1‑150 (MQ=255)                                                                                  tatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACg                                                               <  2:133805/150‑1 (MQ=255)                                                                                     gCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACgggg                                                            >  2:981527/1‑150 (MQ=255)                                                                                     gCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGa                                                           <  2:38184/151‑1 (MQ=255)                                                                                          aaTGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGgc                                                      >  1:841212/1‑151 (MQ=255)                                                                                             gCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGgcgc                                                    >  1:431034/1‑150 (MQ=255)                                                                                                  gAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCAtt                                              <  1:337731/151‑1 (MQ=255)                                                                                                  gAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCAtt                                              <  2:424448/151‑1 (MQ=255)                                                                                                   aaaaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTg                                             <  2:250711/151‑1 (MQ=255)                                                                                                   aaaaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTg                                             >  1:888712/1‑151 (MQ=255)                                                                                                    aaaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTg                                             >  1:308243/1‑150 (MQ=255)                                                                                                     aaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAg                                           <  1:198244/151‑1 (MQ=255)                                                                                                      aGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAgg                                          <  2:870544/151‑1 (MQ=255)                                                                                                            gttgttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGtc                                     <  1:970834/150‑1 (MQ=255)                                                                                                                 tAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGtcattca                                <  1:818337/150‑1 (MQ=255)                                                                                                                  aaaGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCata                              <  1:95608/151‑1 (MQ=255)                                                                                                                    aGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATg                            <  2:4526/151‑1 (MQ=255)                                                                                                                       gTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGaa                          >  2:728387/1‑150 (MQ=255)                                                                                                                          gcgcAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGaaaaa                       <  1:666805/150‑1 (MQ=255)                                                                                                                              aaCATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAGCGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATaa                  <  1:819185/151‑1 (MQ=255)                                                                                                                               aCATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATaaa                 <  2:509531/151‑1 (MQ=255)                                                                                                                                 aTACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAAt               >  1:16742/1‑151 (MQ=255)                                                                                                                                  tACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAAtt              >  2:89967/1‑151 (MQ=255)                                                                                                                                   aCCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAAtt              <  1:187145/150‑1 (MQ=255)                                                                                                                                    ccGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAATTCa            <  2:691960/151‑1 (MQ=255)                                                                                                                                     cgcgCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAATTc             >  2:741187/1‑149 (MQ=255)                                                                                                                                     cgcgCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAATTc             >  2:145348/1‑149 (MQ=255)                                                                                                                                     cgcgCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATGCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAATTCAg           <  1:562315/151‑1 (MQ=255)                                                                                                                                        gcAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAATTCAgcct        >  1:589199/1‑151 (MQ=255)                                                                                                                                        gcAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAATTCAgcct        <  2:244316/151‑1 (MQ=255)                                                                                                                                         cAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAATTCAgcctg       <  2:931096/151‑1 (MQ=255)                                                                                                                                          aaaTGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAATTCAgcctgg      <  1:971069/151‑1 (MQ=255)                                                                                                                                          aaaTGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAATTCAgcctgg      >  1:449253/1‑151 (MQ=255)                                                                                                                                          aaaTGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAATTCAgcct        <  1:69728/149‑1 (MQ=255)                                                                                                                                              gATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAATTCAgcctggcctg  <  2:623854/151‑1 (MQ=255)                                                                                                                                              gATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAATTCAgcctggcctg  >  2:61419/1‑151 (MQ=255)                                                                                                                                              gATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAATTCAgcctggcctg  <  1:9209/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                CAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAATTCAGCCTGGCCTG  >  NC_000913/4294218‑4294509

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