Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,287,910	T→C	H49H (CAT→CAC)	nrfA →	cytochrome c552 nitrite reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,287,910	0	T	C	100.0%	284.6 / NA	79	H49H (CAT→CAC)	nrfA	cytochrome c552 nitrite reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/50);  total (29/50)

ATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCG  >  NC_000913/4287773‑4288036                                                                                                                                          |                                                                                                                               aTAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAATACCTGTAAGTGTGGAAGCGAAGAATCAAACTTTTGCCCCGCAGCAccc                                                                                                                              >  2:792637/1‑140 (MQ=255)  tAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCAcccc                                                                                                                             <  1:148225/140‑1 (MQ=255)   aaaaaTAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCg                                                                                                                            <  2:1040564/140‑1 (MQ=255)     aaaTAAACGCACGCCGTATCTTGAGTTTAGTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGAt                                                                                                                          <  1:810034/140‑1 (MQ=255)     aaaTAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGCAGCGAAGCATGGAACCTTTGCACCGCAGCACCCCGAt                                                                                                                          >  1:498150/1‑140 (MQ=255)      aaTAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATc                                                                                                                         <  2:60844/140‑1 (MQ=255)           gacagcgcCGTATCTTCAGCTTATTGATTCCATTTTTCTCTTTCACTTATGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAatat                                                                                                                    <  2:585020/135‑1 (MQ=255)           acgcacgcCGTATCTTCCGCTTTTTGATTCCCTTTTTCTGTTTCACTTTTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAatat                                                                                                                    <  2:798416/140‑1 (MQ=255)           acgcacgcCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAatat                                                                                                                    <  1:1763374/140‑1 (MQ=255)           acgcacgcCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAatat                                                                                                                    <  1:1160538/140‑1 (MQ=255)           acgcacgcCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAatat                                                                                                                    <  1:1668920/140‑1 (MQ=255)           acgcacgcCGTATCTTCAGCTTATTGATTCCATTTATCTATTTCACTTATGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAatat                                                                                                                    <  2:647506/140‑1 (MQ=255)           acgcacgcCGTATCTCCAGCTTCTTTATTTCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAatat                                                                                                                    <  1:540159/140‑1 (MQ=255)           acgcacgcCGTATCTACAGCTTATTAATTTCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAatat                                                                                                                    <  1:711682/140‑1 (MQ=255)              cacgcCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATAtctc                                                                                                                 <  1:163023/140‑1 (MQ=255)                cgcCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCGTTGCCCCGCAGCACCCCGATCAATAtctctc                                                                                                               >  2:625235/1‑140 (MQ=255)                 gcCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTcc                                                                                                              <  1:1575951/140‑1 (MQ=255)                  ccGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAACGAAGAATAAAACCTTAGCCCCGCAGCACCCCGATCAATATCTCTCCt                                                                                                             >  2:793324/1‑140 (MQ=255)                         ttCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGc                                                                                                      <  1:420428/140‑1 (MQ=255)                           cAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCGCCCGGAGCAATATGTCTCCTGGaactgtc                                                                                                    >  1:585020/1‑135 (MQ=255)                           cAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCCCCCGGACCAATATGTCTCCTGGAAAGCCa                                                                                                    >  1:627447/1‑140 (MQ=255)                             gCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCAcc                                                                                                  <  1:1548386/140‑1 (MQ=255)                                tATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCg                                                                                               >  1:57113/1‑140 (MQ=255)                                 aTTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCgg                                                                                              >  2:1828389/1‑140 (MQ=255)                                 aTTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCGTTGCCCCGAAGCTCCCCGATCAATACCTCTTCTGGAAAGCCACCTCgg                                                                                              >  2:571204/1‑140 (MQ=255)                                   tGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGATCAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGaa                                                                                            <  1:696434/140‑1 (MQ=255)                                   tGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGGATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGaa                                                                                            <  2:1156607/140‑1 (MQ=255)                                        ccTTTTTTCTTTTTCGCTTCTGGTCGCGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTc                                                                                       <  1:793324/140‑1 (MQ=255)                                        ccTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTc                                                                                       >  1:95412/1‑140 (MQ=255)                                        ccTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTc                                                                                       >  2:30851/1‑140 (MQ=255)                                        ccTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTc                                                                                       >  1:1679774/1‑140 (MQ=255)                                        ccTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTc                                                                                       >  2:1906350/1‑140 (MQ=255)                                           tttttctttttcACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCaga                                                                                    >  2:1877570/1‑140 (MQ=255)                                               tctttttcACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCgt                                                                                >  1:1101899/1‑140 (MQ=255)                                                   tttcACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTg                                                                            <  2:1682875/140‑1 (MQ=255)                                                   tttcACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTg                                                                            >  2:1737297/1‑140 (MQ=255)                                                   tttcACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTg                                                                            <  2:1752040/140‑1 (MQ=255)                                                   tttcACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTg                                                                            >  1:1263349/1‑140 (MQ=255)                                                        cTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGcc                                                                       >  2:1223770/1‑140 (MQ=255)                                                         ttCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGCAACTGTGAAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACATTCAGAGCGTGTTGACGccc                                                                      >  2:1661670/1‑140 (MQ=255)                                                          tCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCt                                                                     <  1:1906350/140‑1 (MQ=255)                                                            tGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTgg                                                                   <  2:132113/140‑1 (MQ=255)                                                                cACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGa                                                               <  1:278269/140‑1 (MQ=255)                                                                cACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGa                                                               <  2:1263349/140‑1 (MQ=255)                                                                cACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGa                                                               <  1:1877570/140‑1 (MQ=255)                                                                   gCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCAACTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGa                                                            >  2:706294/1‑140 (MQ=255)                                                                        gcAAACGGCTGCTCCCGTAAAACCTGTAACTTTGGCAGCCAAGAATGAGACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCAc                                                       <  2:542154/139‑1 (MQ=255)                                                                        aCAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCAc                                                       <  1:1632387/140‑1 (MQ=255)                                                                        aCAAACGGCTGCTCCCGCAAAACCTGTAACTGAGGAAGAGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATGAAAATCTCTTCTGGAAAGCCACCTCGGAACAGTTATAGCGTGTTGACGCCCTGGAGGAAGATCGAc                                                       >  2:1169750/1‑140 (MQ=255)                                                                         cAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACg                                                      <  1:1072498/140‑1 (MQ=255)                                                                         cAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACg                                                      <  1:1223770/140‑1 (MQ=255)                                                                           aaCGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGc                                                    <  1:1737297/140‑1 (MQ=255)                                                                           aaCGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGc                                                    <  2:320010/140‑1 (MQ=255)                                                                           aaCGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGc                                                    <  1:1828389/140‑1 (MQ=255)                                                                           aaCGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGAGTGTTGACGCCCTGGCGGAAGATCCACGGc                                                    <  1:1319355/140‑1 (MQ=255)                                                                              ggCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTgg                                                 <  2:1569281/140‑1 (MQ=255)                                                                                ctgctCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTg                                               <  1:30851/140‑1 (MQ=255)                                                                                 tgctCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGa                                              >  1:1836455/1‑140 (MQ=255)                                                                                 tgctCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGa                                              >  1:1217103/1‑140 (MQ=255)                                                                                    tCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGAGTTGACGGCCTCGCGGAAGATACACGGCTGGTGATcc                                           >  1:542154/1‑140 (MQ=255)                                                                                          aaaaCCTGTAACTGTGGAAGCGGAGGATGGAACCTTTTCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTggg                                     <  1:686515/140‑1 (MQ=255)                                                                                          aaaaCCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTggg                                     <  2:95412/140‑1 (MQ=255)                                                                                          aaaaCCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTggg                                     <  2:254690/140‑1 (MQ=255)                                                                                            aaCCTGTAACTGTGGAAGCGAACAATGAATCCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCg                                   <  1:613868/140‑1 (MQ=255)                                                                                             aCCTGTAACTGTGGAAGAGAAGTATGAAACCTTTGCCCCCCAGCACCCAGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCgg                                  <  2:522096/140‑1 (MQ=255)                                                                                              ccTGTAACTGTGGAAGCGAAGACTGTAACTTTTGCCCTGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCggg                                 <  1:571204/140‑1 (MQ=255)                                                                                                tGTAACTGTGGAAGCGAAGAATGAAGCCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGt                               <  1:625235/140‑1 (MQ=255)                                                                                                        tgGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCgggg                                >  1:1096270/1‑131 (MQ=255)                                                                                                        tgGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCgggg                                <  2:1096270/131‑1 (MQ=255)                                                                                                        tgGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTtc                       >  2:249921/1‑140 (MQ=255)                                                                                                        tgGAAGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATC‑ATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTtct                      >  1:1598039/1‑140 (MQ=255)                                                                                                           aaGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTcg                    <  2:260859/140‑1 (MQ=255)                                                                                                           aaGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTcg                    <  1:249921/140‑1 (MQ=255)                                                                                                            aGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTcgc                   <  2:1101899/140‑1 (MQ=255)                                                                                                            aGCGAAGAATGAAACCTTTGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTcgc                   <  1:706294/140‑1 (MQ=255)                                                                                                             gCGAAGAATCAAACCTTTGCCCCCCAGCACACCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTcgcg                  <  1:1169750/140‑1 (MQ=255)                                                                                                                aagaaTGAAACCTTTGCCCCGCGGCACCCCGCTCATTATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGa               <  2:627447/140‑1 (MQ=255)                                                                                                                             ttGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTcgc                   <  2:1769346/123‑1 (MQ=255)                                                                                                                             ttGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTcgc                   >  1:1769346/1‑123 (MQ=255)                                                                                                                             ttGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCcg  >  2:1206655/1‑140 (MQ=255)                                                                                                                             ttGCCCCGCAGCACCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCcg  >  1:1749049/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                               ATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCG  >  NC_000913/4287773‑4288036

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑