Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,622,393	C→T	G354G (GGC→GGT)	ydeE →	dipeptide exporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,622,393	0	C	T	96.8%	201.5 / ‑2.3	62	G354G (GGC→GGT)	ydeE	dipeptide exporter
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (29/31);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TCGGTGAAATCATTTATGCGCCGGGCGAGTATATGTTGATTGACCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGCGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGTTAAAAGGGATTCGAGCAAGACCGTGGGGGCAGCCCGCGCTTTGTTGATTTAAGTCG  >  NC_000913/1622260‑1622528                                                                                                                                      |                                                                                                                                        tCGGTGAAATCATTTATGCGCCGGGCGAGTATATGTTGATTGACCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTg                                                                                                                                   >  2:627915/1‑140 (MQ=255)  cGGTGAAATCATTTATGCGCCGGGCGAGTATATGTTGATTGACCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGc                                                                                                                                  >  1:1361614/1‑140 (MQ=255)    gTGAAATCATTTATGCGCCGGGCGAGTATATGTTGATTGACCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTa                                                                                                                                >  2:1148177/1‑140 (MQ=255)    gTGAAATCATTTATGCGCCGGGCGAGTATATGTTGATTGACCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTa                                                                                                                                >  1:1297831/1‑140 (MQ=255)       aaaTCATTTATGCGCCGGGCGAGTATATGTTGATTGACCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAAcc                                                                                                                             <  2:1410234/140‑1 (MQ=255)        aaTCATTTATGCGCCGGGCGAGTATATGTTGATTGACCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCa                                                                                                                            <  2:880240/140‑1 (MQ=255)          tCATTTATGCGCCGGGCGAGTATATGTTGATTGACCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGc                                                                                                                          <  2:762037/140‑1 (MQ=255)               tATGCGCCGGGCGAGTATATGTTGATTGACCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTgcc                                                                                                                     <  1:1366553/140‑1 (MQ=255)                        ggCGAGTATATGTTGATTGACCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTc                                                                                                            >  2:775932/1‑140 (MQ=255)                         gCGAGTATATGTTGATTGACCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCg                                                                                                           <  1:452579/140‑1 (MQ=255)                          cGAGTATATGTTGATTGACCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGc                                                                                                          <  1:1058643/140‑1 (MQ=255)                          cGAGTATATGTTGATTGACCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGc                                                                                                          <  1:468623/140‑1 (MQ=255)                          cGAGTATATGTTGATTGACCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGc                                                                                                          <  1:1148177/140‑1 (MQ=255)                            aGTATATGTTGATTGACCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTg                                                                                                        <  1:292672/140‑1 (MQ=255)                               atatGTTGATTGACCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGttt                                                                                                     <  2:640847/140‑1 (MQ=255)                                tatGTTGATTGACCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTg                                                                                                    <  2:304490/140‑1 (MQ=255)                                       attgaCCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGt                                                                                                       <  2:537999/130‑1 (MQ=255)                                       attgaCCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGt                                                                                                       >  1:537999/1‑130 (MQ=255)                                       attgaCCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCtt                                                                                             >  1:1341431/1‑140 (MQ=255)                                       attgaCCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCtt                                                                                             >  2:1210737/1‑140 (MQ=255)                                        ttgaCCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTa                                                                                            >  2:539716/1‑140 (MQ=255)                                        ttgaCCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTa                                                                                            >  2:933397/1‑140 (MQ=255)                                         tgaCCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAg                                                                                           <  2:1282539/140‑1 (MQ=255)                                           aCCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCg                                                                                         >  2:135115/1‑140 (MQ=255)                                             cATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCtt                                                                                             <  1:232010/134‑1 (MQ=255)                                             cATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCtt                                                                                             >  2:232010/1‑134 (MQ=255)                                             cATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCg                                                                                         <  2:1064173/138‑1 (MQ=255)                                             cATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCg                                                                                         >  1:1064173/1‑138 (MQ=255)                                              atatTGCGCCGCCAGAAATGAAGGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTg                                                                                      <  1:1187450/140‑1 (MQ=255)                                               tatTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGCGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTgg                                                                                     >  2:847213/1‑140 (MQ=255)                                                atTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGt                                                                                    <  2:163710/140‑1 (MQ=255)                                                atTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGt                                                                                    <  1:539716/140‑1 (MQ=255)                                                  tGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGa                                                                                  <  1:775932/140‑1 (MQ=255)                                                      ccgccAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCAt                                                                                                >  2:291745/1‑122 (MQ=255)                                                      ccgccAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCAt                                                                                                <  1:291745/122‑1 (MQ=255)                                                       cgccAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCAtt                                                                             >  2:807288/1‑140 (MQ=255)                                                       cgccAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCAtt                                                                             <  2:1361614/140‑1 (MQ=255)                                                         ccAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATtgc                                                                           >  1:1055112/1‑140 (MQ=255)                                                                           tATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGtt                                                         >  1:513470/1‑140 (MQ=255)                                                                              ttttCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGTTaaa                                                      <  2:1297831/140‑1 (MQ=255)                                                                                    gCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGTTAAAAGGGAt                                                <  2:1197421/140‑1 (MQ=255)                                                                                     cccAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGTTAAAAGGGAtt                                               >  2:715397/1‑140 (MQ=255)                                                                                      ccAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGTTAAAAGGGATTc                                              <  1:399447/140‑1 (MQ=255)                                                                                             ttAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGCGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGTTAAAAGGGATTCGAGCAAg                                       <  1:847213/140‑1 (MQ=255)                                                                                                    ggctTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGTTAAAAGGGAt                                                >  2:537028/1‑124 (MQ=255)                                                                                                    ggctTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGTTAAAAGGGAt                                                <  1:537028/124‑1 (MQ=255)                                                                                                    ggctTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGTTAAAAGGGATTCGAGCAAGACCGTgg                                >  1:480356/1‑140 (MQ=255)                                                                                                    ggctTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGTTAAAAGGGATTCGAGCAAGACCGTgg                                >  2:1194718/1‑140 (MQ=255)                                                                                                      ctTGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTTCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGTTAAAAGGGATTCGAGCAAGACCGTgggg                              >  2:934272/1‑140 (MQ=255)                                                                                                       ttGGTGCCGCGATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGTTAAAAGGGATTCGAGCAAGACCGTggggg                             <  1:1210737/140‑1 (MQ=255)                                                                                                               gcgATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGTTAAAAGGGATTCGAGCAAGACCGTGGGGGCAGCCcgc                     >  1:360873/1‑140 (MQ=255)                                                                                                                cgATTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGTTAAAAGGGATTCGAGCAAGACCGTGGGGGCAGCAcgcg                    <  2:1084236/140‑1 (MQ=255)                                                                                                                  aTTAACCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGTTAAAAGGGATTCGAGCAAGACCGTGGGGGCAGCCCGCGCt                  <  2:1055112/140‑1 (MQ=255)                                                                                                                     aaCCCATTAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGTTAAAAGGGATTCGAGCAAGACCGTGGGGGCAGCCCGCGCTttg               >  1:505439/1‑140 (MQ=255)                                                                                                                            tAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGTTAAAAGGGATTCGAGCAAGACCGTGGGGGCAGCCCGCGCTTTGTTGAttt        <  1:135115/140‑1 (MQ=255)                                                                                                                            tAGTGAGTGGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGTTAAAAGGGATTCGAGCAAGACCGTGGGGGCAGCCCGCGCTTTGTTGAttt        <  1:627915/140‑1 (MQ=255)                                                                                                                                 agtgGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGTTAAAAGGGATTCGAGCAAGACCGTGGGGGCAGCCCGCGCTTTGTTGATTTAAGTc   >  2:121959/1‑140 (MQ=255)                                                                                                                                 agtgGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGTTAAAAGGGATTCGAGCAAGACCGTGGGGGCAGCCCGCGCTTTGTTGATTTAAGTc   >  2:600010/1‑140 (MQ=255)                                                                                                                                 agtgGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGTTAAAAGGGATTCGAGCAAGACCGTGGGGGCAGCCCGCGCTTTGTTGATTTAAGTc   >  1:228148/1‑140 (MQ=255)                                                                                                                                 agtgGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGTTAAAAGGGATTCGAGCAAGACCGTGGGGGCAGCCCGCGCTTTGTTGATTTAAGTc   >  2:1140521/1‑140 (MQ=255)                                                                                                                                 agtgGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGTTAAAAGGGATTCGAGCAAGACCGTGGGGGCAGCCCGCGCTTTGTTGATTTAAGTc   <  1:95367/140‑1 (MQ=255)                                                                                                                                 agtgGTGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGCGATCATTGCTGCGTGGGTGCTGATGTTAAAAGGGATTCGAGCAAGACCGTGGGGGCAGCCCGCGCTTTGTTGATTTAAGTc   >  1:1224045/1‑140 (MQ=255)                                                                                                                                  gtggtgtggtGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGTTAAAAGGGATTCGAGCAAGACCGTGGGGGCAGCCCGCGCTTTGTTGATTTAAGTCg  >  2:729150/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                        TCGGTGAAATCATTTATGCGCCGGGCGAGTATATGTTGATTGACCATATTGCGCCGCCAGAAATGAAAGCCAGCTATTTTTCCGCCCAGTCTTTAGGCTGGCTTGGTGCCGCGATTAACCCATTAGTGAGTGGCGTAGTGCTAACCAGCCTGCCGCCTTCCTCGCTGTTTGTCATCTTAGCGTTGGTGATCATTGCTGCGTGGGTGCTGATGTTAAAAGGGATTCGAGCAAGACCGTGGGGGCAGCCCGCGCTTTGTTGATTTAAGTCG  >  NC_000913/1622260‑1622528

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