Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,507,424	G→C	pseudogene (228/255 nt)	insO →	IS911B regulator fragment

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,507,424	0	G	C	100.0%	141.2 / NA	43	pseudogene (228/255 nt)	insO	IS911B regulator fragment
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (22/21);  total (22/21)

CAGGTGCTCCAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTGTTCTCATTTGTTCCCCGCTGCCAAGGACTCTGAGGA  >  NC_000913/4507186‑4507670                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       caGGTGCTCCAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATgg                                                                                                                                                                                                                                              >  1:196895/1‑249 (MQ=255)      gCTCCAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGa                                                                                                                                                                                                                                       <  1:420081/251‑1 (MQ=255)        tCCAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGAAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGaa                                                                                                                                                                                                                                      >  2:255959/1‑250 (MQ=255)           aaTGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAata                                                                                                                                                                                                                                   <  1:344012/250‑1 (MQ=255)           aaTGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAata                                                                                                                                                                                                                                   >  2:344012/1‑250 (MQ=255)                        aaTTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCtccc                                                                                                                                                                                                                     <  2:359941/251‑4 (MQ=255)                         aTTTCAGCGCAGAGTTTAAACGCGAAGCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccg                                                                                                                                                                                                                    >  2:170281‑M2/1‑247 (MQ=255)                          tttCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGCAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATAAAATATTAAAAAGGCTAcagc                                                                                                                                                                                                                   >  2:220094/1‑248 (MQ=255)                          tttCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAATGCTaccgc                                                                                                                                                                                                                   <  1:362316‑M2/251‑6 (MQ=255)                                 gcAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttg                                                                                                                                                                                                            <  1:382048‑M2/251‑13 (MQ=255)                                  cAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAACCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGTAGAATGAAATATTAAAAAGGCTaccgcgctcttg                                                                                                                                                                                                            <  2:124424‑M2/250‑13 (MQ=255)                                         tAAACGCGAATCCGCTCAACTGGTTGTGGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTGGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcag                                                                                                                                                                                                     >  1:391066‑M2/1‑231 (MQ=255)                                          aaaCGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAAATGCGTGATGAGCGTAAGGGCAAAACACCAAAAGCCTCTCCGAAAACACCAGAACAAATCGAAATACGCACGCTGAGGAAAAAGCTACAAAGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaacgcgctcttgatgtcagac                                                                                                                                                                                                   <  1:401614‑M2/251‑22 (MQ=255)                                                     cGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacag                                                                                                                                                                                        >  1:213978‑M2/1‑219 (MQ=255)                                                           aCTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAACTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcg                                                                                                                                                                                  >  1:141995‑M2/1‑213 (MQ=255)                                                                         aGAAATACACGGTGGCAGATGCCGCGAAAGATATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaact                                                                                                                                                                    >  2:254863‑M2/1‑199 (MQ=255)                                                                                           aTGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgt                                                                                                                                                  >  1:356294‑M2/1‑181 (MQ=255)                                                                                               cgcgAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcggg                                                                                                                                                                         >  1:260354‑M2/1‑177 (MQ=255)                                                                                               cgcgAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcggg                                                                                                                                                                         <  2:260354‑M2/224‑48 (MQ=255)                                                                                                       cTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctg                                                                                                                                      <  2:277005‑M2/251‑83 (MQ=255)                                                                                                            gATGTTGGCCTTTCCACAATGTCAAGATGGGTCAAACAACTGCGTGAGGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAAACGAAATACGTAAGCTGAGGAAAAAGCTACAACGCAGTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataagcgggaaactcagagcgcattatcctgtggtcacactctgccatg                                                                                                                                 >  2:258563‑M2/1‑164 (MQ=255)                                                                                                             aTGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgt                                                                                                                                <  1:28037‑M2/251‑89 (MQ=255)                                                                                                               gTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAACCCTCTCCGCTAACACCAGAACAAATCGAAATACGTAAGCTGTGGCAAAAGCTACAACGCATTGAACTGGAGAATGCAATATTAAAAACCCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactccgagcgcattatcctgtggtcacactctgccatgtgt                                                                                                                              <  1:57868‑M2/251‑91 (MQ=255)                                                                                                                 gggcCTTTCCCCAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAAGACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataagcgggaaactcagagcgcattatcctgtggtcacactgtgccatgtggtc                                                                                                                            >  2:479923‑M2/2‑159 (MQ=255)                                                                                                                  agcCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAACGCTACAACGCATTGAAATGGAGCATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactacctgaacagttctcgctaatcgggaacctcagagcgcattatcctggggtcacactctgccatgtgttgg                                                                                                                           <  2:120163‑M2/250‑94 (MQ=255)                                                                                                                                  cAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgtttgggtttcatcgcagcag                                                                                                           >  1:389409‑M2/1‑142 (MQ=255)                                                                                                                                  cAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcag                                                                                                           >  1:392168‑M2/1‑142 (MQ=255)                                                                                                                                  cAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcag                                                                                                           >  2:278177‑M2/1‑142 (MQ=255)                                                                                                                                            cAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagata                                                                                                  >  2:69958‑M2/1‑132 (MQ=255)                                                                                                                                                 aacTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCGGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacacgctgccatgtgttcggggttcatcgcagcagctacagatactggaa                                                                                            >  2:98098‑M2/1‑127 (MQ=255)                                                                                                                                                               cGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaa                                                                               <  1:355013‑M2/250‑138 (MQ=255)                                                                                                                                                                    gggCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaacca                                                                         >  1:400229‑M2/1‑108 (MQ=255)                                                                                                                                                                    gggCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaacc                                                                          <  1:160142‑M2/250‑143 (MQ=255)                                                                                                                                                                                       tctcCGATAACACCAGAACACATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagacgggctgtat                                                      <  2:292144‑M2/251‑163 (MQ=255)                                                                                                                                                                                                 cacCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcct                                                                                                                                                    <  1:34077‑M2/147‑69 (MQ=255)                                                                                                                                                                                                 cacCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcct                                                                                                                                                    >  2:34077‑M2/1‑79 (MQ=255)                                                                                                                                                                                                   ccAGAACAAATCGAAATACGTAAGCTGTGGAAAATGCTACAACGCATTGAAATGGAGAATGCAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttgatcgcagcagctacagataccggcaagaccgtcctgaaaaaccagacggcagacgggctgtatgacgcagtcagg                                          <  2:186989‑M2/251‑175 (MQ=255)                                                                                                                                                                                                         cAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagacgggctgtattacgcagtcaggtacttg                                    >  1:159350‑M2/1‑71 (MQ=255)                                                                                                                                                                                                              cGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagacgggctgtattacgcagtcaggtacttgagct                                <  1:30963‑M2/250‑185 (MQ=255)                                                                                                                                                                                                                      gTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagacgggctgtattacgcagtcaggaacttgagctacatggcat                       <  2:119704‑M2/251‑194 (MQ=255)                                                                                                                                                                                                                         aGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagacgggctgtattacgcagtcaggtacttgagctacatggcatc                      <  1:194710‑M2/249‑195 (MQ=255)                                                                                                                                                                                                                             gAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctctttatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgc                                                                                                                                     >  1:283296‑M2/1‑51 (MQ=255)                                                                                                                                                                                                                                           aaCGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagacgggctgtattacgcagtcaggtacttgagctacatggcatcagccacggttcggccggagc  <  1:259680‑M2/251‑215 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       CAGGTGCTCCAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTGTTCTCATTTGTTCCCCGCTGCCAAGGACTCTGAGGA  >  NC_000913/4507186‑4507670

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