Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,587,702	A→G	T252A (ACT→GCT)	mrr →	methylated adenine and cytosine restriction protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,587,702	0	A	G	100.0%	273.5 / NA	75	T252A (ACT→GCT)	mrr	methylated adenine and cytosine restriction protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/41);  total (34/41)

ATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGG  >  NC_000913/4587566‑4587831                                                                                                                                         |                                                                                                                                  atCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCtt                                                                                                                                >  1:1117692/1‑140 (MQ=255)    gCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTg                                                                                                                             >  2:192292/1‑140 (MQ=255)    gCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTg                                                                                                                             >  1:1177030/1‑140 (MQ=255)     cTTGATAAACTTGGCCTGGGGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAATGTGGGGTGTTTATTACCGCTTCTgg                                                                                                                            >  1:1116638/1‑140 (MQ=255)     cTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTgg                                                                                                                            >  1:1154015/1‑140 (MQ=255)     cTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTgg                                                                                                                            >  2:251386/1‑140 (MQ=255)     cTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTgg                                                                                                                            >  2:259544/1‑140 (MQ=255)     cTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTgg                                                                                                                            >  1:273564/1‑140 (MQ=255)     cTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTgg                                                                                                                            >  1:214222/1‑140 (MQ=255)         aTAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTc                                                                                                                               >  1:271938/1‑133 (MQ=255)         aTAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTc                                                                                                                               <  2:271938/133‑1 (MQ=255)          tAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCt                                                                                                                                 >  1:995992/1‑130 (MQ=255)              cTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTtc                                                                                                                   >  2:924175/1‑140 (MQ=255)                  gCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAg                                                                                                               >  1:915525/1‑140 (MQ=255)                  gCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAg                                                                                                               <  2:1177030/140‑1 (MQ=255)                            aGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACt                                                                                                     >  1:186215/1‑140 (MQ=255)                                     tCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAATGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAAt                                                                                            <  2:1116638/140‑1 (MQ=255)                                     tCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAAt                                                                                            <  2:214222/140‑1 (MQ=255)                                     tCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAAt                                                                                            <  1:1225296/140‑1 (MQ=255)                                     tCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAAt                                                                                            <  1:259544/140‑1 (MQ=255)                                     tCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAAt                                                                                            <  1:251386/140‑1 (MQ=255)                                       aGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATcc                                                                                          >  1:844022/1‑140 (MQ=255)                                        ggCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCg                                                                                         >  1:732113/1‑140 (MQ=255)                                         gCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAAt                                                                                            <  2:77131/136‑1 (MQ=255)                                         gCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAAt                                                                                            >  1:77131/1‑136 (MQ=255)                                         gCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGt                                                                                        <  1:825869/140‑1 (MQ=255)                                         gCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGt                                                                                        <  1:63949/140‑1 (MQ=255)                                           aaaaCGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCg                                                                                      <  2:273564/140‑1 (MQ=255)                                           aaaaCGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCg                                                                                      <  1:192292/140‑1 (MQ=255)                                           aaaaCGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCg                                                                                      <  1:789384/140‑1 (MQ=255)                                            aaaCGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGa                                                                                     <  2:1117692/140‑1 (MQ=255)                                                gTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGt                                                                                 <  2:732113/140‑1 (MQ=255)                                                  tGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAAt                                                                                            >  2:1205365/1‑127 (MQ=255)                                                  tGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAAt                                                                                            <  1:1205365/127‑1 (MQ=255)                                                  tGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTAt                                                                               <  2:90368/140‑1 (MQ=255)                                                   ggCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGa                                                                                     <  1:321628/133‑1 (MQ=255)                                                   ggCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGa                                                                                     >  2:321628/1‑133 (MQ=255)                                                       gAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGgtgt                                                                          >  2:630217/1‑140 (MQ=255)                                                             tGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTgttggttg                                                                    <  2:197025/140‑1 (MQ=255)                                                              gTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGa                                                                   >  1:529275/1‑140 (MQ=255)                                                              gTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGa                                                                   <  2:1089362/140‑1 (MQ=255)                                                                  gcaggcCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATggg                                                               >  2:119297/1‑140 (MQ=255)                                                                  gcaggcCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATggg                                                               >  1:561359/1‑140 (MQ=255)                                                                   caggcCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATgggg                                                              <  2:561359/140‑1 (MQ=255)                                                                   caggcCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATgggg                                                              <  1:636165/140‑1 (MQ=255)                                                                    aggcCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGa                                                             >  1:1152767/1‑140 (MQ=255)                                                                    aggcCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGa                                                             <  2:186215/140‑1 (MQ=255)                                                                    aggcCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGa                                                             <  2:1152767/140‑1 (MQ=255)                                                                     ggcCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGaa                                                            >  1:443927/1‑140 (MQ=255)                                                                      gcCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAAc                                                           >  1:938849/1‑140 (MQ=255)                                                                      gcCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAAc                                                           >  2:1006790/1‑140 (MQ=255)                                                                      gcCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAAc                                                           >  1:958768/1‑140 (MQ=255)                                                                      gcCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAAc                                                           >  1:39874/1‑140 (MQ=255)                                                                       ccAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAAc                                                           >  1:12429/1‑139 (MQ=255)                                                                       ccAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAAc                                                           <  2:12429/139‑1 (MQ=255)                                                                       ccAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACg                                                          >  1:32851/1‑140 (MQ=255)                                                                        cAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGc                                                         <  2:141756/140‑1 (MQ=255)                                                                            aTTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTgg                                                     <  2:616763/140‑1 (MQ=255)                                                                            aTTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTgg                                                     <  1:1019921/140‑1 (MQ=255)                                                                            aTTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTgg                                                     <  1:631496/140‑1 (MQ=255)                                                                            aTTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCACGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTgg                                                     <  1:937293/140‑1 (MQ=255)                                                                             ttACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGt                                                    <  2:915525/140‑1 (MQ=255)                                                                              tACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTg                                                   <  2:1154015/140‑1 (MQ=255)                                                                              tACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTg                                                   <  1:280346/140‑1 (MQ=255)                                                                                cAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATgggg                                                              >  2:594625/1‑127 (MQ=255)                                                                                cAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATgggg                                                              <  1:594625/127‑1 (MQ=255)                                                                                cAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCa                                                 <  2:181984/140‑1 (MQ=255)                                                                                          aCGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATc                                       <  1:462326/140‑1 (MQ=255)                                                                                               gcACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGaaaa                                  <  1:119297/140‑1 (MQ=255)                                                                                                         ggCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAggg                        <  2:32851/140‑1 (MQ=255)                                                                                                          gCAAAAAGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAgggg                       >  2:1013973/1‑140 (MQ=255)                                                                                                               aaGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTttct                  >  2:123928/1‑140 (MQ=255)                                                                                                                aGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTttctt                 <  1:123928/140‑1 (MQ=255)                                                                                                                aGCGAAACGTGGGGTGTTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTttctt                 <  1:1006790/140‑1 (MQ=255)                                                                                                                               gtTTATTACCGCTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAgg  <  2:529275/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                  ATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGG  >  NC_000913/4587566‑4587831

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