Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,241,611	G→C	10.4%	noncoding (35/76 nt)	PP_t25 →	tRNA‑Ala

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,241,611	0	G	C	10.4%	190.5 / 11.5	65	noncoding (35/76 nt)	PP_t25	tRNA‑Ala
Reads supporting (aligned to +/- strand):  ref base G (28/30);  new base C (4/3);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

AAATTAGTTTGACAGCCCTGCAATCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGTGGCCTAGGACACCGC  >  NC_002947/2241465‑2241755                                                                                                                                                   |                                                                                                                                                 aaaTTAGTTTGACAGCCCTGCAATCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTgcatggcat                                                                                                                                               >  2:284957/1‑150 (MQ=255)  aaTTAGTTTGACAGCCCTGCAATCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTgcatggcatg                                                                                                                                              >  2:244695/1‑150 (MQ=255)     tAGTTTGACAGCCCTGCAATCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCag                                                                                                                                           <  1:2040173/150‑1 (MQ=255)     tAGTTTGACAGCCCTGCAATCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCag                                                                                                                                           <  2:1766255/150‑1 (MQ=255)            aCAGCCCTGCAATCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCa                                                                                                                                    <  1:2073257/150‑1 (MQ=255)              aGCCCTGCAGTACGAAATAAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGc                                                                                                                                  >  2:1332928/1‑150 (MQ=255)                   tGCAATCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTc                                                                                                                             <  1:431100/150‑1 (MQ=255)                   tGCAATCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGATATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTc                                                                                                                             <  1:1921170/150‑1 (MQ=255)                      aaTCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGag                                                                                                                          >  1:855720/1‑149 (MQ=255)                          cGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGAAGAGCAGCGGTTAGATCACg                                                                                                                      >  2:1890934/1‑150 (MQ=38)                               ttAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAg                                                                                                                 >  1:1422353/1‑150 (MQ=255)                                   tatGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTcc                                                                                                             >  2:467790/1‑150 (MQ=255)                                     tGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTccac                                                                                                           >  1:49770/1‑150 (MQ=255)                                     tGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGAATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTccac                                                                                                           <  2:1663729/150‑1 (MQ=255)                                       cgcCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTccacca                                                                                                         >  2:848028/1‑150 (MQ=255)                                       cgcCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTccacca                                                                                                         >  2:892999/1‑150 (MQ=255)                                        gcCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCaa                                                                                                        >  1:915767/1‑150 (MQ=255)                                         ccccGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCaaa                                                                                                       >  2:1143269/1‑150 (MQ=255)                                          cccGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAAt                                                                                                      >  2:391210/1‑150 (MQ=255)                                           ccGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCaggag                                                                                                                                        >  2:2342645/1‑115 (MQ=255)                                           ccGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCaggag                                                                                                                                        <  1:2342645/115‑1 (MQ=255)                                            cGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCaaa                                                                                                       <  1:1440026/147‑1 (MQ=255)                                            cGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCaaa                                                                                                       >  2:1440026/1‑147 (MQ=255)                                              tCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCg                                                                                                  <  1:391210/150‑1 (MQ=255)                                               ccACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGc                                                                                                 <  2:915767/150‑1 (MQ=255)                                                  caGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCg                                                                                                                                 >  2:1634956/1‑115 (MQ=255)                                                  caGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCg                                                                                                                                 <  1:1634956/115‑1 (MQ=255)                                                    gCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTagc                                                                                            >  1:1801278/1‑150 (MQ=255)                                                     cGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTagca                                                                                           <  1:892999/150‑1 (MQ=255)                                                        tGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCt                                                                                        >  2:881539/1‑150 (MQ=255)                                                          aacgaaaacgaaACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTaa                                                                                      <  1:467790/150‑1 (MQ=255)                                                            cgaaaacgaaaCGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAgg                                                                                    <  2:1153391/150‑1 (MQ=255)                                                             gaaaacgaaaCGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCaaa                                                                                                       >  2:875393/1‑130 (MQ=255)                                                             gaaaacgaaaCGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCaaa                                                                                                       <  1:875393/130‑1 (MQ=255)                                                             gaaaacgaaaCGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGc                                                                                   >  1:157636/1‑150 (MQ=255)                                                               aaacgaaacgCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTagca                                                                                           >  1:2086070/1‑140 (MQ=255)                                                               aaacgaaacgCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTagca                                                                                           <  2:2086070/140‑1 (MQ=255)                                                                aacgaaacgCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGa                                                                                <  1:848028/150‑1 (MQ=255)                                                                     aacgCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCAGTTCTATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGc                                                                           <  1:1240626/150‑1 (MQ=255)                                                                          cAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGATCAGCTGGGAGAGCACAGGTTTAGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGAc                                                                      <  1:1332928/150‑1 (MQ=37)                                                                            ggCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCg                                                                    <  1:1822333/150‑1 (MQ=255)                                                                               aCCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATa                                                                              <  2:163103/137‑1 (MQ=255)                                                                               aCCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATa                                                                              >  1:163103/1‑137 (MQ=255)                                                                                cccAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTc                                                                <  1:867718/150‑1 (MQ=255)                                                                                     ccagcGATTCAGTTCAAAGCTACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCa                                                           >  2:572105/1‑150 (MQ=255)                                                                                      cagcGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCaa                                                          >  2:1859533/1‑150 (MQ=255)                                                                                       agcGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAAc                                                         >  1:1556831/1‑150 (MQ=255)                                                                                                  tCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGaa                                              <  1:1541887/150‑1 (MQ=255)                                                                                                  tCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGaa                                              <  2:1556831/150‑1 (MQ=255)                                                                                                    aaaGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAgg                                            >  1:2056994/1‑150 (MQ=255)                                                                                                      aGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTgcatggcatg                                                                                                                                              >  1:1037995/1‑50 (MQ=255)                                                                                                      aGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTgcatggcatg                                                                                                                                              <  2:1037995/50‑1 (MQ=255)                                                                                                            aTTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAgg                                    >  1:159715/1‑150 (MQ=255)                                                                                                             ttGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCTGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGt                                   <  1:1143269/150‑1 (MQ=255)                                                                                                               gtgGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTc                                 <  2:1801278/150‑1 (MQ=255)                                                                                                               gtgGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTc                                 <  2:159715/150‑1 (MQ=255)                                                                                                               gtgGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTc                                 <  1:244695/150‑1 (MQ=255)                                                                                                                      tataGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCt                          <  2:49770/150‑1 (MQ=255)                                                                                                                            tCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCt                    >  1:1457198/1‑150 (MQ=255)                                                                                                                            tCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCt                    >  1:1445392/1‑150 (MQ=255)                                                                                                                                      gagCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGTGGCCTAg          <  1:572105/150‑1 (MQ=255)                                                                                                                                         ctccTGCATGCCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCgg                                                                   >  2:785435/3‑90 (MQ=255)                                                                                                                                         ctccTGCATGCCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCgg                                                                   <  1:785435/88‑1 (MQ=255)                                                                                                                                         ctccTGCATGCCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGt                 <  2:1014759/138‑1 (MQ=255)                                                                                                                                         ctccTGCATGCCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGt                 >  1:1014759/3‑140 (MQ=255)                                                                                                                                          tccTGCATGCCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGaa                                              >  2:1955216/2‑110 (MQ=255)                                                                                                                                          tccTGCATGCCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGaa                                              <  1:1955216/109‑1 (MQ=255)                                                                                                                                          tccTGCATGCCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGTGGCCTAGGaca      >  2:1384004/2‑150 (MQ=255)                                                                                                                                            cTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGTGGCCTAGGACAcc    >  2:1275241/1‑150 (MQ=255)                                                                                                                                              gcatgccatgCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGTGGCCTAGGACACCGc  >  2:1986780/1‑150 (MQ=255)                                                                                                                                                   gcatgCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGTTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGt                 <  2:2391470/130‑1 (MQ=255)                                                                                                                                                   gcatgCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGTTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGt                 >  1:2391470/1‑130 (MQ=255)                                                                                                                                                   |                                                                                                                                                 AAATTAGTTTGACAGCCCTGCAATCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGTGGCCTAGGACACCGC  >  NC_002947/2241465‑2241755

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