Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,578,851	A→G	S163P (TCA→CCA)	mcrB ←	5‑methylcytosine‑specific restriction enzyme subunit McrB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,578,851	0	A	G	100.0%	234.6 / NA	66	S163P (TCA→CCA)	mcrB	5‑methylcytosine‑specific restriction enzyme subunit McrB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/31);  total (35/31)

GTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGT  >  NC_000913/4578715‑4578984                                                                                                                                         |                                                                                                                                      gTTTTTCCAACGCCGAGCGGCCGCTAGAGGATAATATGTTTTTTGATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCAt                                                                                                                                    <  1:378164/140‑1 (MQ=255)    tttCCAACGCCGGGCGGGCCCTTGAGGATAATATTCTTTTTGATGGTTTATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAgg                                                                                                                                 <  2:271446/140‑1 (MQ=255)    tttCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAgg                                                                                                                                 <  1:525206/140‑1 (MQ=255)      tCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTg                                                                                                                               <  2:173128/140‑1 (MQ=255)           cGCCGGGCGGCCCCTGGAGGATCATATTTTTTTTGGTGGTTAATCGTTCGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATa                                                                                                                          <  2:303088/140‑1 (MQ=255)           cGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATa                                                                                                                          >  1:182800/1‑140 (MQ=255)                ggcggcCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACAct                                                                                                                     <  1:358848/140‑1 (MQ=255)                                taataTTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAat                                                                                                     >  1:687646/1‑140 (MQ=255)                                taataTTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAat                                                                                                     >  2:702479/1‑140 (MQ=255)                                         ttttGATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTta                                                                                            <  2:480578/140‑1 (MQ=255)                                                       cGTTTGAGGATCGTCTCAATTGGGGTTTCATGGATTAAC‑TATCATTTAGCGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTAtg                                                                             <  2:403282/140‑1 (MQ=255)                                                        gTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTAtg                                                                             <  1:702479/140‑1 (MQ=255)                                                          ttGAGTATCGTCTCTTTTGTGGCTTCAGGAATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTAtgtt                                                                           <  1:368206/140‑1 (MQ=255)                                                           tGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTAtgttg                                                                          >  2:61054/1‑140 (MQ=255)                                                                  cGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCa                                                                   <  2:182800/140‑1 (MQ=255)                                                                    tctcTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTTCCAGAATTCAATCTTAATTTATAGTTGTTGATTATGTTGTCCAGCATa                                                                 >  2:347082/1‑140 (MQ=255)                                                                    tctcTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATa                                                                 >  1:790935/1‑140 (MQ=255)                                                                               gTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCg                                                      >  1:800301/1‑140 (MQ=255)                                                                               gTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCg                                                      >  1:593948/1‑140 (MQ=255)                                                                               gTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTACAGTCGTTAATTATGATGTCCAACATAGAGGCAAATCg                                                      >  1:255115/1‑140 (MQ=255)                                                                                    aGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGgtgt                                                 <  2:687646/140‑1 (MQ=255)                                                                                    aGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGgtgt                                                 <  2:800301/140‑1 (MQ=255)                                                                                    aGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGgtgt                                                 >  2:86719/1‑140 (MQ=255)                                                                                    aGGGATAAACAAATCATATAACCCATCTTCCTGACCGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGgtgt                                                 <  2:255115/140‑1 (MQ=255)                                                                                         tAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTGTAGTCGTTAATTATGTATTCCAGCATAGAGACAAATCGGGTGTAATCa                                            >  2:403035/1‑140 (MQ=255)                                                                                                aTCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCgg                                                     >  1:9919/1‑124 (MQ=255)                                                                                                aTCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCgg                                                     <  2:9919/124‑1 (MQ=255)                                                                                                aTCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCt                                     >  2:730681/1‑140 (MQ=255)                                                                                                  cATTTAACGCCTCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCtgt                                   >  2:613650/1‑140 (MQ=255)                                                                                                   aTTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCtgtg                                  <  1:86719/140‑1 (MQ=255)                                                                                                    tttAACGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCACATCGGGTGTGATCAATACCCTGTga                                 >  1:371218/1‑140 (MQ=255)                                                                                                       aaCGCATCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTCGATTACGTTGTCCAGCATCGAGCCAAATCGGTTGTAATCAATACCCTGTGAGAc                              >  2:290707/1‑140 (MQ=255)                                                                                                            aTCTTCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGGAAATCGGGTGTAATCAATACCCTGTGAGACttttt                         >  2:289920/1‑140 (MQ=255)                                                                                                               ttCCAGACAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTggg                      >  2:205679/1‑140 (MQ=255)                                                                                                                     caggTATGATTCAGTTTTTGGCCTAGGTGGACTAACACTCTTGCCAGACTTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAgg                <  2:209101/137‑1 (MQ=255)                                                                                                                     aCAGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAgg                <  2:790935/140‑1 (MQ=255)                                                                                                                       aGTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCg              <  1:205679/140‑1 (MQ=255)                                                                                                                        gTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAgg                <  2:583769/137‑1 (MQ=255)                                                                                                                        gTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAgg                >  1:583769/1‑137 (MQ=255)                                                                                                                        gTATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGACAAATCGGTTGTAATCAATACCCTGCGAGAGTTTTTGGGAACAGGCGt             >  2:342363/1‑140 (MQ=255)                                                                                                                         tATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAgg                >  1:731075/1‑136 (MQ=255)                                                                                                                         tATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAgg                <  2:731075/136‑1 (MQ=255)                                                                                                                         tATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTa            >  2:752987/1‑140 (MQ=255)                                                                                                                         tATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTa            >  2:593115/1‑140 (MQ=255)                                                                                                                         tATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTa            >  2:146925/1‑140 (MQ=255)                                                                                                                         tATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTa            >  2:415893/1‑140 (MQ=255)                                                                                                                         tATGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTa            >  2:444955/1‑140 (MQ=255)                                                                                                                          aTGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTaa           <  2:593951/140‑1 (MQ=255)                                                                                                                           tGATTCAGTTTTTGGCGTAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAAt          >  1:537230/1‑140 (MQ=255)                                                                                                                           tGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGa                     >  1:110676/1‑129 (MQ=255)                                                                                                                           tGATTCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGa                     <  2:110676/129‑1 (MQ=255)                                                                                                                              ttCAGTTTTTGGCATTGGTGGATTACCACTCTTGCCAAAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGAGAGACTTTTTGGGAACAGGCGTAATAgg       <  1:403035/140‑1 (MQ=255)                                                                                                                               tCAGTTTTTGGCATAGGTGGAATAACGCTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGa      <  1:613650/140‑1 (MQ=255)                                                                                                                               tCAGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGa      <  1:444955/140‑1 (MQ=255)                                                                                                                                 aGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGa                     >  2:200318/1‑123 (MQ=255)                                                                                                                                 aGTTTTTGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGa                     <  1:200318/123‑1 (MQ=255)                                                                                                                                   tttttGGCATAGGTGGAATCACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGaa                    <  1:288189/122‑1 (MQ=255)                                                                                                                                   tttttGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACGCTGTGAGGCTTGTTGGGAACATGCGTAATAGGAcctg  >  1:209101/1‑137 (MQ=255)                                                                                                                                   tttttGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGaa                    >  2:288189/1‑122 (MQ=255)                                                                                                                                   tttttGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAgg                >  2:203716/1‑126 (MQ=255)                                                                                                                                   tttttGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAgg                <  1:203716/126‑1 (MQ=255)                                                                                                                                   tttttGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAAt          >  2:585411/1‑132 (MQ=255)                                                                                                                                   tttttGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAAt          <  1:585411/132‑1 (MQ=255)                                                                                                                                   tttttGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGt  >  1:523724/1‑140 (MQ=255)                                                                                                                                     tttGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCg              >  2:200001/1‑126 (MQ=255)                                                                                                                                     tttGGCATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCg              <  1:200001/126‑1 (MQ=255)                                                                                                                                         |                                                                                                                                      GTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGT  >  NC_000913/4578715‑4578984

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: ‑