Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,220,971	G→A	pseudogene (1371/2648 nt)	ycgH →	putative transporter component, N‑terminal fragment

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,220,971	0	G	A	100.0%	245.5 / NA	69	pseudogene (1371/2648 nt)	ycgH	putative transporter component, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (44/25);  total (44/25)

ATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAATCTCCAGCAG  >  NC_000913/1220835‑1221107                                                                                                                                         |                                                                                                                                         aTTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCt                                                                                                                                       <  2:154934/140‑1 (MQ=255)    cAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACg                                                                                                                                    >  2:480009/1‑140 (MQ=255)    cAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACg                                                                                                                                    >  2:589588/1‑140 (MQ=255)      gACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGa                                                                                                                                  <  2:721505/140‑1 (MQ=255)       aCAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGAc                                                                                                                                 >  1:1111082/1‑140 (MQ=255)         aGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAg                                                                                                                               <  1:589588/140‑1 (MQ=255)           tattatCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGAc                                                                                                                             >  2:523762/1‑140 (MQ=255)                   aTCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTAttgtt                                                                                                                     <  1:1056039/140‑1 (MQ=255)                     cTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTAttgt                                                                                                                      >  2:698635/1‑137 (MQ=255)                     cTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTAttgt                                                                                                                      <  1:698635/137‑1 (MQ=255)                     cTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGa                                                                                                                   >  2:108485/1‑140 (MQ=255)                     cTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGa                                                                                                                   >  1:723068/1‑140 (MQ=255)                                  atgGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCg                                                                                                      >  2:52234/1‑140 (MQ=255)                                  atgGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCg                                                                                                      >  1:167243/1‑140 (MQ=255)                                  atgGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCg                                                                                                      <  2:295164/140‑1 (MQ=255)                                  atgGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCg                                                                                                      >  1:768739/1‑140 (MQ=255)                                   tgGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCgg                                                                                                     >  2:977274/1‑140 (MQ=255)                                      tCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGtt                                                                                                  <  1:334768/140‑1 (MQ=255)                                        aaaCATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCa                                                                                                >  2:809353/1‑140 (MQ=255)                                          aCATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCggg                                                                                                    >  2:935178/1‑134 (MQ=255)                                          aCATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCggg                                                                                                    >  2:1086965/1‑134 (MQ=255)                                          aCATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCggg                                                                                                    <  1:935178/134‑1 (MQ=255)                                          aCATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCggg                                                                                                    <  1:1086965/134‑1 (MQ=255)                                               tACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCg                                                                                         >  1:36438/1‑140 (MQ=255)                                                      ggAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGtt                                                                                                  <  2:455958/124‑1 (MQ=255)                                                      ggAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGtt                                                                                                  >  1:455958/1‑124 (MQ=255)                                                      ggAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTAtgt                                                                                  >  1:377699/1‑140 (MQ=255)                                                      ggAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTAtgt                                                                                  >  1:1126423/1‑140 (MQ=255)                                                      ggAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTAtgt                                                                                  >  2:994860/1‑140 (MQ=255)                                                      ggAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACCTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTAtgt                                                                                  >  1:147763/1‑140 (MQ=255)                                                                  aCCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTgg                                                                      >  2:430139/1‑140 (MQ=255)                                                                        aaaGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAg                                                                <  1:279571/140‑1 (MQ=255)                                                                            gTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAgcggc                                                            >  2:586734/1‑140 (MQ=255)                                                                                ttATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAgcggcgcgg                                                        <  1:586734/140‑1 (MQ=255)                                                                                     tCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACg                                                   <  2:138847/140‑1 (MQ=255)                                                                                     tCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTCATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACg                                                   <  1:147368/140‑1 (MQ=255)                                                                                       ggTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCt                                                 >  2:302532/1‑140 (MQ=255)                                                                                         taataaTGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTaa                                               >  2:1049453/1‑140 (MQ=255)                                                                                          aataatGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTaaa                                              >  1:1136843/1‑140 (MQ=255)                                                                                            taatGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACg                                            <  1:523762/140‑1 (MQ=255)                                                                                               tGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGa                                         >  1:637447/1‑140 (MQ=255)                                                                                                 gTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATc                                       <  2:1189794/140‑1 (MQ=255)                                                                                                 gTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATc                                       >  1:1189794/1‑140 (MQ=255)                                                                                                  tCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCg                                      >  2:1004301/1‑140 (MQ=255)                                                                                                  tCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCg                                      >  1:757001/1‑140 (MQ=255)                                                                                                   cAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGa                                     >  2:870318/1‑140 (MQ=255)                                                                                                    aGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGaa                                    >  1:410902/1‑140 (MQ=255)                                                                                                      cTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACt                                  >  1:980110/1‑140 (MQ=255)                                                                                                      cTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACt                                  >  1:688155/1‑140 (MQ=255)                                                                                                        tAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCa                                <  2:1111082/140‑1 (MQ=255)                                                                                                        tAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCa                                <  2:688155/140‑1 (MQ=255)                                                                                                        tAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCa                                <  1:302532/140‑1 (MQ=255)                                                                                                                gCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGa                                         <  1:978403/123‑1 (MQ=255)                                                                                                                gCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGa                                         >  2:978403/1‑123 (MQ=255)                                                                                                                gCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCa                                <  1:279458/132‑1 (MQ=255)                                                                                                                gCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCa                                >  2:279458/1‑132 (MQ=255)                                                                                                                gCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACtgtg                        >  1:1040403/1‑140 (MQ=255)                                                                                                                   ccGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAAt                     <  1:108485/140‑1 (MQ=255)                                                                                                                      tATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACt                                  >  1:835362/1‑124 (MQ=255)                                                                                                                      tATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGc                  >  1:499100/1‑140 (MQ=255)                                                                                                                      tATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGc                  >  1:263900/1‑140 (MQ=255)                                                                                                                      tATTAGGTGATGATAAATCAGCTACGAACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACt                                  <  2:835362/124‑1 (MQ=255)                                                                                                                              gatgatAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGc                  >  1:163202/1‑132 (MQ=255)                                                                                                                              gatgatAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGc                  <  2:163202/132‑1 (MQ=255)                                                                                                                               atgatAAATCAGCTGCGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAAtct         >  2:572722/1‑140 (MQ=255)                                                                                                                                tgatAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAAtctc        >  1:445893/1‑140 (MQ=255)                                                                                                                                tgatAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAAtctc        >  2:455610/1‑140 (MQ=255)                                                                                                                                tgatAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGACTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAAtctc        >  2:225513/1‑140 (MQ=255)                                                                                                                                      aTCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAATCTCcagcag  <  1:809353/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                         ATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAATCTCCAGCAG  >  NC_000913/1220835‑1221107

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