Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,870,267	T→C	R297R (CGT→CGC)	cbrA →	colicin M resistance protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,870,267	0	T	C	100.0%	221.7 / NA	62	R297R (CGT→CGC)	cbrA	colicin M resistance protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/31);  total (31/31)

GCTGTTTCCCTCGCGCTGGCAGGATTTTGTCTGCGGTAAGGACAACGCCTTTCTGATTGGTGAAGCGGCGGGATTTATCAGCGCCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGTTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATG  >  NC_000913/3870132‑3870402                                                                                                                                         |                                                                                                                                        gCTGTTTCCCTCGCGCTGGCAGGATTTTGTCTGCGGTAAGGACAACGCCTTTCTGATTGGTGAAGCGGCGGGATTTATCAGCGCCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCgg                                                                                                                                     <  1:88809/140‑1 (MQ=255)   tGTTTCCCTCGCGCTGGCAGGATTTTGTCTGCGGTAAGGACAACGCCTTTCTGATTGGTGAAGCGGCGGGATTTATCAGCGCCAGCTCGCTGGAGGGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGgt                                                                                                                                    <  2:435982/140‑1 (MQ=255)        cccTCGCGCTGGCAGGATTTTGTCTGCGGTAAGGACAACGCCTTTCTGATTGGTGAAGCGGCGGGATTTATCAGCGCCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTc                                                                                                                                       >  1:603370/1‑131 (MQ=255)          cTCGCGCTGGCAGGATTTTGTCTGCGGTAAGGACAACGCCTTTCTGATTGGTGAAGCGGCGGGATTTATCAGCGCCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGa                                                                                                                            >  1:157214/1‑140 (MQ=255)          cTCGCGCTGGCAGGATTTTGTCTGCGGTAAGGACAACGCCTTTCTGATTGGTGAAGCGGCGGGATTTATCAGCGCCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGa                                                                                                                            >  2:83768/1‑140 (MQ=255)           tCGCGCTGGCAGGATTTTGTCTGCGGTAAGGACAACGCCTTTCTGATTGGTGAAGCGGCGGGATTTATCAGCGCCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGaa                                                                                                                           <  2:310414/140‑1 (MQ=255)           tCGCGCTGGCAGGATTTTGTCTGCGGTAAGGACAACGCCTTTCTGATTGGTGAAGCGGCGGGATTTATCAGCGCCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGaa                                                                                                                           <  2:157214/140‑1 (MQ=255)           tCGCGCTGGCAGGATTTTGTCTGCGGTAAGGACAACGCCTTTCTGATTGGTGAAGCGGCGGGATTTATCAGCGCCAGCTCGCTGGA‑GCGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGaa                                                                                                                           <  1:729016/140‑1 (MQ=255)                     aGGATTTTGTCTGCGGTAAGGACAACGCCTTTCTGATTGGTGAAGCGGCGGGATTTATCAGCGCCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCagag                                                                                                                 <  1:891518/140‑1 (MQ=255)                     aGGATTTTGTCTGCGGTAAGGACAACGCCTTTCTGATTGGTGAAGCGGCGGGATTTATCAGCGCCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCagag                                                                                                                 <  1:620989/140‑1 (MQ=255)                     aGGATTTTGTCTGCGGTAAGGACAACGCCTTTCTGATTGGTGAAGCGGCGGGATTTATCAGCGCCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCagag                                                                                                                 <  2:522879/140‑1 (MQ=255)                     aGGATTTTGTCTGCGGTAAGGACAACGCCTTTCTGATTGGTGAAGCGGCGGGATTTATCAGCGCCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCagag                                                                                                                 <  1:93140/140‑1 (MQ=255)                         ttttGTCTGCGGTAAGGACAACGCCTTTCTGATTGGTGAAGCGGCGGGATTTATCAGCGCCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGaa                                                                                                                           <  2:87089/126‑1 (MQ=255)                         ttttGTCTGCGGTAAGGACAACGCCTTTCTGATTGGTGAAGCGGCGGGATTTATCAGCGCCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGaa                                                                                                                           >  1:87089/1‑126 (MQ=255)                          tttGTCTGCGGTAAGGACAACGCCTTTCTGATTGGTGAAGCGGCGGGATTTATCAGCGCCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCt                                                                                                            <  1:549712/140‑1 (MQ=255)                                            aaCGCCTTTCTGATTGGTGAAGCGGCGGGATTTATCAGCGCCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGgcg                                                                                          >  2:479667/1‑140 (MQ=255)                                               gCCTTTCTGATTGGTGAAGCGGCGGGATTTATCAGCGCCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGgcgcgc                                                                                       <  2:484390/140‑1 (MQ=255)                                               gCCTTTCTGATTGGTGAAGCGGCGGGATTTATCAGCGCCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGgcgcgc                                                                                       >  1:484390/1‑140 (MQ=255)                                               gCCTTTCTGATTGGTGAAGCGGCGGGATTTATCAGCGCCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGgcgcgc                                                                                       <  1:83768/140‑1 (MQ=255)                                               gCCTTTCTGATTGGTGAAGCGGCGGGATTTATCAGCGCCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGgcgcgc                                                                                       <  1:811041/140‑1 (MQ=255)                                                 cTTTCTGATTGGTGAAGCGGCGGGATTTATCAGCGCCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCa                                                                                     >  1:642845/1‑140 (MQ=255)                                                    tCTGATTGGTGAAGCGGCGGGATTTATCAGCGCCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCAccc                                                                                  >  2:1065183/1‑140 (MQ=255)                                                         ttGGTGAAGCGGCGGGATTTATCAGCCCCAGCCCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCaaa                                                                             >  1:268657/1‑140 (MQ=255)                                                                       ggATTTATCAGCGCCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAAct                                                               >  1:1044572/1‑140 (MQ=255)                                                                              tCAGCGCCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGc                                                        >  2:303984/1‑140 (MQ=255)                                                                                aGCGCCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCaa                                                      <  1:479667/140‑1 (MQ=255)                                                                                 gcgcCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAg                                                     >  2:337588/1‑140 (MQ=255)                                                                                  cgcCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCgg                                                         <  1:1051520/135‑1 (MQ=255)                                                                                  cgcCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCgg                                                         >  2:1051520/1‑135 (MQ=255)                                                                                  cgcCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCGGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCgg                                                         >  1:927473/1‑135 (MQ=255)                                                                                  cgcCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCGGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCgg                                                         <  2:927473/135‑1 (MQ=255)                                                                                       gCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTG‑‑GCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGa                                                    >  2:238962/1‑133 (MQ=255)                                                                                       gCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTG‑‑GCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGa                                                    <  1:238962/133‑1 (MQ=255)                                                                                        cTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTaa                                              >  2:781574/1‑140 (MQ=255)                                                                                         tCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTaaa                                             <  2:889806/140‑1 (MQ=255)                                                                                             tGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGc                                         <  1:781574/140‑1 (MQ=255)                                                                                                   gggATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATg                                    >  2:92974/1‑140 (MQ=255)                                                                                                   gggATTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATg                                    <  1:92974/140‑1 (MQ=255)                                                                                                    ggATTAGCTACGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGAAAGTAAAAAGCCGATGc                                   <  2:896628/140‑1 (MQ=255)                                                                                                      aTTAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCt                                 >  1:597380/1‑140 (MQ=255)                                                                                                        tAGCTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGa                               >  1:777546/1‑140 (MQ=255)                                                                                                           cTATGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCg                            >  2:417406/1‑140 (MQ=255)                                                                                                             aTGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCa                          <  2:112145/140‑1 (MQ=255)                                                                                                              tGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGc                                         >  2:991471/1‑124 (MQ=255)                                                                                                              tGCGCTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGc                                         <  1:991471/124‑1 (MQ=255)                                                                                                               gcgcTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGc                                         <  2:263406/123‑1 (MQ=255)                                                                                                               gcgcTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGc                                         >  1:263406/1‑123 (MQ=255)                                                                                                               gcgcTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCAcc                        >  1:140832/1‑140 (MQ=255)                                                                                                               gcgcTGGATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCAcc                        >  1:931200/1‑140 (MQ=255)                                                                                                                    ggATAGCACAGACATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCtt                   <  2:931200/140‑1 (MQ=255)                                                                                                                               aCATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTa                 >  2:394198/1‑131 (MQ=255)                                                                                                                               aCATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTa                 <  1:394198/131‑1 (MQ=255)                                                                                                                               aCATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTgg        >  2:702084/1‑140 (MQ=255)                                                                                                                               aCATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTgg        >  1:373042/1‑140 (MQ=255)                                                                                                                               aCATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTgg        <  1:389305/140‑1 (MQ=255)                                                                                                                               aCATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTgg        >  1:432479/1‑140 (MQ=255)                                                                                                                               aCATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTgg        >  1:491308/1‑140 (MQ=255)                                                                                                                               aCATTCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTgg        >  2:237668/1‑140 (MQ=255)                                                                                                                                  ttCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGAtt     >  2:1003271/1‑140 (MQ=255)                                                                                                                                  ttCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGAtt     <  1:1065183/140‑1 (MQ=255)                                                                                                                                  ttCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGAtt     <  2:396682/140‑1 (MQ=255)                                                                                                                                   tCTGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTa    >  1:1033058/1‑140 (MQ=255)                                                                                                                                     tGCGCTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATg  <  2:307213/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                        GCTGTTTCCCTCGCGCTGGCAGGATTTTGTCTGCGGTAAGGACAACGCCTTTCTGATTGGTGAAGCGGCGGGATTTATCAGCGCCAGCTCGCTGGA‑GGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGTTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATG  >  NC_000913/3870132‑3870402

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