Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,348,795	T→A	5.4%	intergenic (+27/+27)	yfaE → / ← inaA	ferredoxin involved with ribonucleotide reductase diferric‑tyrosyl radical (Y*) cofactor maintenance/acid‑inducible Kdo/WaaP family putative kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,348,795	0	T	A	5.4%	267.0 / 5.0	74	intergenic (+27/+27)	yfaE/inaA	ferredoxin involved with ribonucleotide reductase diferric‑tyrosyl radical (Y*) cofactor maintenance/acid‑inducible Kdo/WaaP family putative kinase
Reads supporting (aligned to +/- strand):  ref base T (34/36);  new base A (2/2);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

CGCACACGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGCAAGCGGCG  >  NC_000913/2348649‑2348938                                                                                                                                                   |                                                                                                                                                cGCACACGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCtt                                                                                                                                                <  4:271562/148‑1 (MQ=255)  gCACACGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTa                                                                                                                                               <  4:246907/148‑1 (MQ=255)     cacGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGc                                                                                                                                            >  2:238608/1‑148 (MQ=255)     cacGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGc                                                                                                                                            >  1:142998/1‑148 (MQ=255)     cacGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGc                                                                                                                                            >  3:20325/1‑148 (MQ=255)       cGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTa                                                                                                                                               <  3:630582/143‑1 (MQ=255)       cGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTa                                                                                                                                               >  4:630582/1‑143 (MQ=255)         ccTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCt                                                                                                                                        >  2:62286/1‑148 (MQ=255)           tGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTa                                                                                                                                      <  3:559883/148‑1 (MQ=255)                 cAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGc                                                                                                                                <  4:20325/148‑1 (MQ=255)                 cAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGc                                                                                                                                <  2:255392/148‑1 (MQ=255)                 cAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGc                                                                                                                                <  2:142998/148‑1 (MQ=255)                  aGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGcc                                                                                                                               <  1:62286/148‑1 (MQ=255)                         gtcgACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTa                                                                                                                                               >  4:898770/4‑125 (MQ=255)                         gTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTa                                                                                                                                               <  3:898770/125‑1 (MQ=255)                              cTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTa                                                                                                                                               >  3:1066137/1‑120 (MQ=255)                              cTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTa                                                                                                                                               <  4:1066137/120‑1 (MQ=255)                                ggATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATAAGACGCATAAGCGTTGTATCTGGca                                                                                                                               >  3:973720/1‑133 (MQ=38)                                ggATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATAAGACGCATAAGCGTTGTATCTGGca                                                                                                                               <  4:973720/134‑2 (MQ=38)                                              ttAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGa                                                                                                                       <  4:956077/128‑1 (MQ=255)                                              ttAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGa                                                                                                                       >  3:956077/1‑128 (MQ=255)                                                    tttATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCa                                                                                             <  1:238608/148‑1 (MQ=255)                                                                   gAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCAcc                                                                                           >  1:1097026/1‑135 (MQ=255)                                                                   gAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCAcc                                                                                           <  2:1097026/135‑1 (MQ=255)                                                                    aaaTTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGc                                                                             >  3:487739/1‑148 (MQ=255)                                                                        tttGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTg                                                                         <  4:487739/148‑1 (MQ=255)                                                                              cTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGa                                                                                                                       >  2:964757/1‑96 (MQ=255)                                                                              cTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGAtt                                                                   >  1:1042816/1‑148 (MQ=255)                                                                              cTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGACGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGa                                                                                                                       <  1:964757/96‑1 (MQ=255)                                                                               tgttgtTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTAc                                                                                                                   >  2:156949/1‑99 (MQ=255)                                                                               tgttgtTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTAc                                                                                                                   <  1:156949/99‑1 (MQ=255)                                                                               tgttgtTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTg                                                                  >  1:533248/1‑148 (MQ=255)                                                                                 ttgttgCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCt                                                                                       >  3:929579/1‑125 (MQ=255)                                                                                 ttgttgCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCt                                                                                       <  4:929579/125‑1 (MQ=255)                                                                                   gttgCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTc                                                              <  2:1042816/148‑1 (MQ=255)                                                                                                      tatTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTgg                                                                        >  1:332729/1‑119 (MQ=255)                                                                                                      tatTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTgg                                                                        <  2:332729/119‑1 (MQ=255)                                                                                                             aTCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTgg                                                                                                                                 <  4:580114/55‑1 (MQ=255)                                                                                                             aTCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTgg                                                                                                                                 >  3:580114/1‑55 (MQ=255)                                                                                                             aTCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTc                                                                                              >  3:227200/1‑90 (MQ=255)                                                                                                             aTCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTc                                                                                              <  4:227200/90‑1 (MQ=255)                                                                                                               cGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTccc                                                                                    >  4:170808/1‑98 (MQ=255)                                                                                                               cGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTccc                                                                                    <  3:170808/98‑1 (MQ=255)                                                                                                                         aTTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCg                                                                              >  4:138312/1‑94 (MQ=255)                                                                                                                         aTTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCg                                                                              <  3:138312/94‑1 (MQ=255)                                                                                                                         aTTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTcc                                                             >  3:231236/1‑111 (MQ=255)                                                                                                                         aTTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTcc                                                             <  4:231236/111‑1 (MQ=255)                                                                                                                           tGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCg                                         >  3:582553/1‑129 (MQ=255)                                                                                                                           tGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCg                                         <  2:977481/129‑1 (MQ=255)                                                                                                                           tGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCg                                         >  1:977481/1‑129 (MQ=255)                                                                                                                           tGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCg                                         <  4:582553/129‑1 (MQ=255)                                                                                                                           tGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCt                      <  2:533248/148‑1 (MQ=255)                                                                                                                                 tcgTGCCAGATACAATGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTc                                                                                      <  2:807755/76‑1 (MQ=255)                                                                                                                                 tAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTc                                                                                      >  1:807755/1‑78 (MQ=255)                                                                                                                                      ccAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGa                                                                     <  1:95518/90‑1 (MQ=255)                                                                                                                                      ccAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGa                                                                     >  2:95518/1‑90 (MQ=255)                                                                                                                                       cAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTg                                                                         >  2:1108076/1‑85 (MQ=255)                                                                                                                                       cAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTg                                                                         <  1:1108076/85‑1 (MQ=255)                                                                                                                                         gataagaCGCATAAGCGTTGTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTa                            <  1:779013/122‑1 (MQ=255)                                                                                                                                         gataagaCGCATAAGCGTTGTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTa                            >  2:779013/7‑128 (MQ=255)                                                                                                                                         gATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATg                                >  3:861179/1‑124 (MQ=255)                                                                                                                                         gATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATg                                <  4:861179/124‑1 (MQ=255)                                                                                                                                               aCGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCg                                                                                                             <  1:1105334/41‑1 (MQ=255)                                                                                                                                               aCGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCg                                                                                                             >  2:1105334/1‑41 (MQ=255)                                                                                                                                               aCGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGAt                                                                    >  2:351296/1‑82 (MQ=255)                                                                                                                                               aCGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGAt                                                                    <  1:351296/82‑1 (MQ=255)                                                                                                                                               aCGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATg                                <  1:323292/118‑1 (MQ=255)                                                                                                                                               aCGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATg                                >  2:323292/1‑118 (MQ=255)                                                                                                                                               aCGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATc                 <  3:844160/133‑1 (MQ=255)                                                                                                                                               aCGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGCAAgcggcg  >  1:137836/1‑148 (MQ=255)                                                                                                                                               aCGCTTATGCGTCTGATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATc                 >  4:844160/1‑133 (MQ=255)                                                                                                                                                cGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATcgcg              <  2:98438/135‑1 (MQ=255)                                                                                                                                                cGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATcgcg              >  1:98438/1‑135 (MQ=255)                                                                                                                                                  cTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCt                                                                                       <  3:982703/60‑1 (MQ=255)                                                                                                                                                  cTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCt                                                                                       >  4:982703/1‑60 (MQ=255)                                                                                                                                                   |                                                                                                                                                CGCACACGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGCAAGCGGCG  >  NC_000913/2348649‑2348938

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