Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	227,023	A→G	noncoding (1265/2904 nt)	rrlH →	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	227,023	0	A	G	96.8%	269.6 / ‑0.4	95	noncoding (1265/2904 nt)	rrlH	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (42/50);  total (44/51)
Fisher's exact test for biased strand distribution p-value = 5.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

TGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCA  >  NC_000913/226934‑227079                                                                                          |                                                         tGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTa                                                       >  2:3130895/1‑93 (MQ=255)  gCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTaa                                                      <  1:2374949/93‑1 (MQ=255)  gCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTaa                                                      <  1:2038147/93‑1 (MQ=255)   cGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAAc                                                     <  1:4494316/93‑1 (MQ=255)     ttgttgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGa                                                   >  1:2581859/1‑93 (MQ=37)     ttgttgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGa                                                   <  2:2144073/93‑1 (MQ=255)      tgttgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGTTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGAt                                                  <  2:3828508/93‑1 (MQ=255)      tgttgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGAt                                                  >  2:2696786/1‑93 (MQ=37)       gttgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATa                                                 <  1:3650062/93‑1 (MQ=255)        ttgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATaa                                                <  1:1102037/93‑1 (MQ=255)        ttgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATaa                                                >  2:3472409/1‑93 (MQ=37)         tgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATaaa                                               <  2:890007/93‑1 (MQ=255)         tgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATaaa                                               <  2:4361535/93‑1 (MQ=255)         tgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATaaa                                               <  2:4368590/93‑1 (MQ=255)         tgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATaaa                                               <  2:3213887/93‑1 (MQ=255)          gggTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAg                                              <  2:3998145/93‑1 (MQ=255)          gggTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAg                                              <  2:866455/93‑1 (MQ=255)          gggTAGGGGAGCGTTCTGTAAGCCTGAGAAGGGGTACTGTGAGGTATGCTGGAGGTATAAGAAGTGCGAATGCTGACATAGGTAACAATaaaa                                              >  2:4217412/1‑92 (MQ=11)           ggTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGc                                             <  2:645767/93‑1 (MQ=37)           ggTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGc                                             <  2:756491/93‑1 (MQ=255)           ggTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGc                                             >  2:4116466/1‑93 (MQ=37)            gTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCg                                            <  2:2928067/93‑1 (MQ=37)             tAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCgg                                           <  2:3169783/93‑1 (MQ=255)             tAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCgg                                           >  2:1744111/1‑93 (MQ=37)             tAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCgg                                           <  1:2961803/93‑1 (MQ=255)             tAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCgg                                           <  1:2975778/93‑1 (MQ=255)             tAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCgg                                           <  2:1009788/93‑1 (MQ=255)              aGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCggg                                          >  1:234764/1‑93 (MQ=37)              aGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCggg                                          <  2:477078/93‑1 (MQ=37)               ggggAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGt                                         <  1:3386961/93‑1 (MQ=37)               ggggAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGt                                         <  1:483444/93‑1 (MQ=255)               ggggAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGt                                         >  2:4311991/1‑93 (MQ=37)               ggggAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGt                                         >  2:4318676/1‑93 (MQ=37)                gggAGCGGTCTGTAAGCCTGTGAAGGTGTCCTGTGAGGTAGGCTGGAGGTAGCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTg                                        <  2:390107/93‑1 (MQ=11)                  gAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGaa                                      >  1:4091644/1‑93 (MQ=37)                  gAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGGATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGaa                                      >  2:1640900/1‑93 (MQ=255)                    gCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGaaaa                                    <  2:3383449/93‑1 (MQ=255)                    gCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGaaaa                                    <  1:1049202/93‑1 (MQ=37)                     cGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAg                                   <  1:2359413/93‑1 (MQ=255)                     cGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAg                                   <  2:1739310/93‑1 (MQ=37)                     cGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAg                                   <  2:3762330/93‑1 (MQ=255)                     cGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAg                                   <  1:1713166/93‑1 (MQ=255)                     cGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAg                                   <  1:3884713/93‑1 (MQ=37)                      gTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGc                                  >  1:74325/1‑93 (MQ=37)                      gTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGc                                  >  1:3871729/1‑93 (MQ=255)                      gTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCAAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGc                                  <  2:3603287/93‑1 (MQ=255)                       ttCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGcc                                 >  1:2873952/1‑93 (MQ=37)                       ttCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGc                                  >  1:3554884/1‑92 (MQ=37)                       ttCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGc                                  <  2:3554884/92‑1 (MQ=37)                       ttCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGcc                                 >  1:2670271/1‑93 (MQ=37)                        tCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGccc                                <  2:3387553/93‑1 (MQ=37)                        tCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGccc                                <  2:3482366/93‑1 (MQ=255)                          tGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGc                              <  2:3183323/93‑1 (MQ=255)                          tGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGc                              <  1:2557017/93‑1 (MQ=37)                          tGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGc                              >  1:2380000/1‑93 (MQ=37)                           gTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCt                             <  1:4487837/93‑1 (MQ=255)                            tAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTc                            >  2:1991290/1‑93 (MQ=37)                            tAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTc                            >  1:4516902/1‑93 (MQ=37)                            tAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGAGAAAGCGGGTGAAAAGCCCGCTc                            >  1:4547062/1‑93 (MQ=255)                             aaGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCg                           >  1:1297949/1‑93 (MQ=37)                               gCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGcc                         <  2:4370787/93‑1 (MQ=37)                               gCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGcc                         >  2:3729068/1‑93 (MQ=37)                               gCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGcc                         >  2:4464196/1‑93 (MQ=37)                               gCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGcc                         >  2:189333/1‑93 (MQ=37)                                 cTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCgg                       >  1:3005918/1‑93 (MQ=37)                                  tgtgAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGa                      >  2:3595958/1‑93 (MQ=255)                                  tgtgAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGa                      >  1:3500202/1‑93 (MQ=255)                                   gtgAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGaa                     <  2:2995181/93‑1 (MQ=37)                                    tgAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAg                    >  2:152873/1‑93 (MQ=255)                                    tgAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAg                    >  1:4365266/1‑93 (MQ=37)                                     gAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGa                   <  2:2386937/93‑1 (MQ=37)                                      aaGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGAc                  >  1:1496763/1‑93 (MQ=255)                                       aGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGAcc                 >  2:275528/1‑93 (MQ=255)                                       aGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGAcc                 >  2:2004799/1‑93 (MQ=255)                                       aGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGAcc                 >  1:2998608/1‑93 (MQ=255)                                       aGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGCGAAAAGCCCGCTCGCCGGAAGAcc                 >  1:4546470/1‑93 (MQ=255)                                        ggTGTACTGTGAGGTATTGTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCa                <  1:1382646/93‑1 (MQ=255)                                        ggTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCa                >  1:1058411/1‑93 (MQ=255)                                         gtgtACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCaa               <  2:2921155/93‑1 (MQ=37)                                         gtgtACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCaa               <  2:2117453/93‑1 (MQ=37)                                         gtgtACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCaa               <  2:3732441/93‑1 (MQ=37)                                           gtACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAgg             >  1:2416088/1‑93 (MQ=255)                                           gtACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAgg             >  1:1732977/1‑93 (MQ=255)                                            tACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAggg            >  1:1076764/1‑93 (MQ=25)                                             aCTGTGAGGTATGCTGTAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGt           <  2:4230140/93‑1 (MQ=255)                                              cTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGtt          >  2:803446/1‑93 (MQ=18)                                                gtgAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGGGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTcc        <  1:157642/93‑1 (MQ=16)                                                 tgAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCt       >  2:603836/1‑93 (MQ=18)                                                 tgAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCt       <  2:4216164/93‑1 (MQ=18)                                                  gAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTg      <  1:2334428/93‑1 (MQ=18)                                                   aGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGt     <  2:1323034/93‑1 (MQ=18)                                                    ggTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTc    >  2:4647918/1‑93 (MQ=18)                                                     gTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTcc   >  2:465489/1‑93 (MQ=17)                                                     gTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTcc   <  1:4283368/93‑1 (MQ=17)                                                      tATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCa  >  2:3721445/1‑93 (MQ=14)                                                      tATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCa  >  1:124459/1‑93 (MQ=14)                                                                                          |                                                         TGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCA  >  NC_000913/226934‑227079

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