Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,963,829	C→T	R134H (CGC→CAC)	gpp ←	guanosine pentaphosphatase/exopolyphosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,963,829	0	C	T	100.0%	260.8 / NA	81	R134H (CGC→CAC)	gpp	guanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (40/41);  total (40/41)

TTACGATCGGCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGCGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCA  >  NC_000913/3963703‑3963963                                                                                                                               |                                                                                                                                       ttACGATCGGCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGcac                                                                                                                            <  1:1568542/139‑1 (MQ=255)      aTCGGCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGcac                                                                                                                            >  1:1044759/1‑134 (MQ=255)      aTCGGCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGcac                                                                                                                            <  2:1044759/134‑1 (MQ=255)         ggCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTACGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAgtg                                                                                                                    >  1:1366145/1‑139 (MQ=255)          gCAAAATAGCGTTCCAGCCAGGTGACGTAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGACACCACCAgtgg                                                                                                                   >  1:1093672/1‑139 (MQ=255)          gCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTCCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGcac                                                                                                                            >  1:1166362/1‑130 (MQ=255)          gCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTCCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGcac                                                                                                                            <  2:1166362/130‑1 (MQ=255)                      tCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGc                                                                                                       <  2:1093672/139‑1 (MQ=255)                      tCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTACGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGc                                                                                                       <  2:1366145/139‑1 (MQ=255)                       ccagccagGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGcc                                                                                                      <  1:1580747/139‑1 (MQ=255)                        cagccagGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGc                                                                                                            <  1:688459/132‑1 (MQ=255)                        cagccagGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGc                                                                                                            >  2:688459/1‑132 (MQ=255)                        cagccagGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGccc                                                                                                     >  2:1143805/1‑139 (MQ=255)                          gccagGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTg                                                                                                   <  1:1445961/139‑1 (MQ=255)                                    gCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGa                                                                                         <  1:1143805/139‑1 (MQ=255)                                     cAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGAc                                                                                        <  2:1603425/139‑1 (MQ=255)                                        cccATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTg                                                                                     >  1:772205/1‑139 (MQ=255)                                         ccATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGc                                                                                    <  2:1723603/139‑1 (MQ=255)                                          cATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGt                                                                                      <  1:1136839/136‑1 (MQ=255)                                          cATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGt                                                                                      >  2:1136839/1‑136 (MQ=255)                                                            cGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCg                                                                                                                                >  1:1300600/1‑76 (MQ=255)                                                            cGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCg                                                                                                                                <  2:1300600/76‑1 (MQ=255)                                                             gAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAgt                                                                                                                     >  1:1316613/1‑86 (MQ=255)                                                             gAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAgt                                                                                                                     <  2:1316613/86‑1 (MQ=255)                                                                 tggtTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGcaccac                                                                                                                         >  1:1686259/1‑78 (MQ=255)                                                                 tggtTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGcaccac                                                                                                                         <  2:1686259/78‑1 (MQ=255)                                                                 tggtTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGt                                                                                      >  1:1699762/1‑113 (MQ=255)                                                                 tggtTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGt                                                                                      <  2:1699762/113‑1 (MQ=255)                                                                 tggtTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTg                                                            >  2:316549/1‑139 (MQ=255)                                                                       gtgCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCAc                                                               <  1:1637039/130‑1 (MQ=255)                                                                       gtgCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCAc                                                               >  2:1637039/1‑130 (MQ=255)                                                                        tgCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGc                                                                                                                              >  2:1392648/1‑66 (MQ=255)                                                                        tgCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGc                                                                                                                              <  1:1392648/66‑1 (MQ=255)                                                                        tgCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGa                                                     <  2:772205/139‑1 (MQ=255)                                                                                tGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCaccacc                                                                                                                        >  2:98973/1‑64 (MQ=255)                                                                                tGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCaccacc                                                                                                                        <  1:98973/64‑1 (MQ=255)                                                                                tGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCAc                                                               <  2:1314791/121‑1 (MQ=255)                                                                                tGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCAc                                                               >  1:1314791/1‑121 (MQ=255)                                                                                  ccGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGaa                                           >  2:695892/1‑139 (MQ=255)                                                                                       tACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCt                                                                                   >  2:974698/1‑94 (MQ=255)                                                                                       tACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCt                                                                                   <  1:974698/94‑1 (MQ=255)                                                                                       tACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTcc                                      <  1:695892/139‑1 (MQ=255)                                                                                                  tACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAg                                                                                                             <  2:371090/57‑1 (MQ=255)                                                                                                  tACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAg                                                                                                             >  1:371090/1‑57 (MQ=255)                                                                                                  tACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGc                                                                                                            >  1:422951/1‑58 (MQ=255)                                                                                                  tACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGc                                                                                                            <  2:422951/58‑1 (MQ=255)                                                                                                  tACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGc                           >  2:1539347/1‑139 (MQ=255)                                                                                                     tGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCAccacca                                                                                                                       <  2:1690873/44‑1 (MQ=255)                                                                                                     tGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCAccacca                                                                                                                       >  1:1690873/1‑44 (MQ=255)                                                                                                       gCTCCGCCTATATGCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATaa                      >  1:1203463/1‑139 (MQ=255)                                                                                                       gCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCa                                                                                                           >  1:1540232/1‑54 (MQ=255)                                                                                                       gCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCa                                                                                                           <  2:1540232/54‑1 (MQ=255)                                                                                                       gCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGa                                                                                         <  1:443749/72‑1 (MQ=255)                                                                                                       gCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGa                                                                                         >  2:443749/1‑72 (MQ=255)                                                                                                           cGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGc                                                                                                       <  1:1017601/54‑1 (MQ=255)                                                                                                           cGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGc                                                                                                       >  2:1017601/1‑54 (MQ=255)                                                                                                           cGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCAc                                                                         >  2:160743/1‑84 (MQ=255)                                                                                                           cGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCAc                                                                         <  1:160743/84‑1 (MQ=255)                                                                                                            gCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGcc                                                                                                      >  1:968750/1‑54 (MQ=255)                                                                                                            gCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGcc                                                                                                      <  2:968750/54‑1 (MQ=255)                                                                                                            gCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATaa                                                                                               >  1:1252509/1‑61 (MQ=255)                                                                                                            gCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATaa                                                                                               <  2:1252509/61‑1 (MQ=255)                                                                                                            gCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAAt                                                                                             <  1:1682160/63‑1 (MQ=255)                                                                                                            gCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAAt                                                                                             >  2:1682160/1‑63 (MQ=255)                                                                                                            gCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCa                 >  2:1357063/1‑139 (MQ=255)                                                                                                                 tatCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCAc                                                               <  1:1813355/88‑1 (MQ=255)                                                                                                                 tatCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCAc                                                               >  2:1813355/1‑88 (MQ=255)                                                                                                                    ccaccaccaGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCAtt         <  1:1539347/139‑1 (MQ=255)                                                                                                                     caccaccaGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCtt                               >  2:1048511/1‑116 (MQ=255)                                                                                                                     caccaccaGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCtt                               <  1:1048511/116‑1 (MQ=255)                                                                                                                      accaccaGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGa       >  2:1361131/1‑139 (MQ=255)                                                                                                                        caccaGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCAccc              >  1:594135/1‑130 (MQ=255)                                                                                                                        caccaGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCAccc              <  2:594135/130‑1 (MQ=255)                                                                                                                          ccaGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGc                                                                                    >  2:971028/1‑58 (MQ=255)                                                                                                                          ccaGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGc                                                                                    <  1:971028/58‑1 (MQ=255)                                                                                                                          ccaGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGc   >  1:1468791/1‑139 (MQ=255)                                                                                                                           caGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTg                                                                    <  2:231318/73‑1 (MQ=255)                                                                                                                           caGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTg                                                                    >  1:231318/1‑73 (MQ=255)                                                                                                                           caGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCAc                                                               <  2:690691/78‑1 (MQ=255)                                                                                                                           caGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCAc                                                               >  1:690691/1‑78 (MQ=255)                                                                                                                           caGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCa  <  1:1357063/139‑1 (MQ=255)                                                                                                                               |                                                                                                                                       TTACGATCGGCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGCGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCA  >  NC_000913/3963703‑3963963

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