Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	123,475	C→T	F153F (TTC→TTT)	aceE →	pyruvate dehydrogenase, decarboxylase component E1, thiamine triphosphate‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	123,475	0	C	T	100.0%	229.5 / NA	71	F153F (TTC→TTT)	aceE	pyruvate dehydrogenase, decarboxylase component E1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (37/34);  total (37/34)

GCGTCCTTCCAGTCTTCCGCAACCATTTATGATGTGTGCTTTAACCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTCCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATC  >  NC_000913/123341‑123597                                                                                                                                       |                                                                                                                           gCGTCCTTCCAGTCTTCCGCAACCATTTATGATGTGTGCTTTAACCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTgg                                                                                                                        <  2:991565/139‑1 (MQ=255)    tCCTTCCAGTCTTCCGCAACCATTTATGATGTGTGCTTTAACCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAg                                                                                                                     >  1:1364574/1‑139 (MQ=255)      cTTCCAGTCTTCCGCAACCATTTATGATGTGTGCTTTAACCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGgt                                                                                                                   >  2:326836/1‑139 (MQ=255)          cAGTCTTCCGCAACCATTTATGATGTGTGCTTTAACCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGgtcgtc                                                                                                               <  2:1610020/139‑1 (MQ=255)           aGTCTTCCGCAACCATTTATGATGTGTGCTTTAACCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCt                                                                                                              <  1:985247/139‑1 (MQ=255)           aGTCTTCCGCAACCATTTATGATGTGTGCTTTAACCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCt                                                                                                              <  2:319308/139‑1 (MQ=255)            gTCTTCCGCAACCATTTATGATGTGTGCTTTAACCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTg                                                                                                             >  1:686147/1‑139 (MQ=255)             tCTTCCGCAACCATTTATGATGTGTGCTTTAACCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGa                                                                                                            >  2:1191859/1‑139 (MQ=255)                 ccGCAACCATTTATGATGTGTGCTTTAACCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCa                                                                                                        >  1:1062628/1‑139 (MQ=255)                                  gtgtgCTTTAACCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACt                                                                                       >  1:476431/1‑139 (MQ=255)                                   tgtgCTTTAACCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACt                                                                                                          <  2:318461/119‑1 (MQ=255)                                   tgtgCTTTAACCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACt                                                                                                          >  1:318461/1‑119 (MQ=255)                                   tgtgCTTTAACCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACtt                                                                                      >  1:641024/1‑139 (MQ=255)                                   tgtgCTTTAACCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACtt                                                                                      <  2:1364574/139‑1 (MQ=255)                                    gtgCTTTAACCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTc                                                                                     >  2:848110/1‑139 (MQ=255)                                     tgCTTTAACCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTcc                                                                                    >  1:1197385/1‑139 (MQ=255)                                     tgCTTTAACCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTcc                                                                                    >  2:88498/1‑139 (MQ=255)                                     tgCTTTAACCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTcc                                                                                    >  1:62193/1‑139 (MQ=255)                                       cTTTAACCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGt                                                                                  >  1:1615181/1‑139 (MQ=255)                                          tAACCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGagcagc                                                                                                >  1:315776/1‑122 (MQ=255)                                          tAACCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGagcagc                                                                                                <  2:315776/122‑1 (MQ=255)                                          tAACCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAg                                                                               >  2:587076/1‑139 (MQ=255)                                           aaCCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAgg                                                                              >  1:278993/1‑139 (MQ=255)                                            aCCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGa                                                                             >  2:803921/1‑139 (MQ=255)                                                 ttcttcCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGTGCGTGTACGCTCGTGCTTTTCTGGAAGgtcg                                                                                                                 >  2:1496797/1‑98 (MQ=255)                                                 ttcttcCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGgtcg                                                                                                                 <  1:1496797/98‑1 (MQ=255)                                                        gTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCaa                                                                 <  2:1162602/139‑1 (MQ=255)                                                         tGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAAt                                                                <  2:1062628/139‑1 (MQ=255)                                                           cACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGgtcgt                                                                                                                >  2:770275/1‑89 (MQ=255)                                                           cACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGgtcgt                                                                                                                <  1:770275/89‑1 (MQ=255)                                                            aCGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGc                                                             <  2:1197385/139‑1 (MQ=255)                                                            aCGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGc                                                             <  1:88498/139‑1 (MQ=255)                                                            aCGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGc                                                             <  2:191736/139‑1 (MQ=255)                                                              gCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGTGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCct                                                           <  1:326836/139‑1 (MQ=255)                                                              gCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCct                                                           <  1:1499424/139‑1 (MQ=255)                                                                    aGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTc                                                     <  1:803921/139‑1 (MQ=255)                                                                    aGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTc                                                     >  1:654982/1‑139 (MQ=255)                                                                          aTGGCGGCGACCAGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAAt                                                                >  2:925173/1‑122 (MQ=255)                                                                          aTGGCGGCGACCAGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACACTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAAt                                                                <  1:925173/122‑1 (MQ=255)                                                                            ggcggcGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGc                                             <  1:1191859/139‑1 (MQ=255)                                                                              cggcgACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGTTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCg                                                                                   >  2:1101735/1‑99 (MQ=255)                                                                              cggcgACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGTTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCg                                                                                   <  1:1101735/99‑1 (MQ=255)                                                                              cggcgACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATcc                                               <  1:886821/135‑1 (MQ=255)                                                                              cggcgACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATcc                                               >  2:886821/1‑135 (MQ=255)                                                                                   aCCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGaa                                      >  1:1081476/1‑139 (MQ=255)                                                                                   aCCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGaa                                      >  1:132679/1‑139 (MQ=255)                                                                                        gTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGa                                 <  2:132679/139‑1 (MQ=255)                                                                                              ttCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCgg                           <  1:848110/139‑1 (MQ=255)                                                                                                  aGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAAtt                       >  2:1607142/1‑139 (MQ=255)                                                                                                  aGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAAtt                       <  1:1540545/139‑1 (MQ=255)                                                                                                      ccACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTgg                   <  1:587076/139‑1 (MQ=255)                                                                                                      ccACATCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTgg                   <  2:676301/139‑1 (MQ=255)                                                                                                         caTCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATcc                                               <  1:626675/108‑1 (MQ=255)                                                                                                         caTCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATcc                                               >  2:626675/1‑108 (MQ=255)                                                                                                          aTCTCCCCGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGt               >  1:61425/1‑139 (MQ=255)                                                                                                                ccGGGCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGa         <  1:1607142/139‑1 (MQ=255)                                                                                                                    gCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACtt                                                                                      >  2:386892/1‑58 (MQ=255)                                                                                                                    gCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACtt                                                                                      <  1:386892/58‑1 (MQ=255)                                                                                                                    gCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACtt                                                                                      <  1:260442/58‑1 (MQ=255)                                                                                                                    gCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACtt                                                                                      >  2:260442/1‑58 (MQ=255)                                                                                                                    gCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGt     <  2:686147/139‑1 (MQ=255)                                                                                                                    gCGTGTACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGt     >  1:156352/1‑139 (MQ=255)                                                                                                                      gtgtACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTAt   >  1:1127100/1‑139 (MQ=255)                                                                                                                       tgtACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATc  >  1:1184352/1‑139 (MQ=255)                                                                                                                       tgtACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATc  >  2:1514899/1‑139 (MQ=255)                                                                                                                         tACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGAcc       <  1:192459/132‑1 (MQ=255)                                                                                                                         tACGCTCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGAcc       >  2:192459/1‑132 (MQ=255)                                                                                                                              tCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACtt                                                                                      >  2:1498799/1‑48 (MQ=255)                                                                                                                              tCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACtt                                                                                      <  1:1498799/48‑1 (MQ=255)                                                                                                                              tCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAgg                                                                              <  2:1414588/56‑1 (MQ=255)                                                                                                                              tCGTGCTTTTCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAgg                                                                              >  1:1414588/1‑56 (MQ=255)                                                                                                                                       |                                                                                                                           GCGTCCTTCCAGTCTTCCGCAACCATTTATGATGTGTGCTTTAACCACTTCTTCCGTGCACGCAACGAGCAGGATGGCGGCGACCTGGTTTACTTCCAGGGCCACATCTCCCCGGGCGTGTACGCTCGTGCTTTCCTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATC  >  NC_000913/123341‑123597

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