Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	335,970	C→T	Q16* (CAA→TAA)	yahE →	DUF2877 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	335,970	0	C	T	100.0%	242.7 / NA	75	Q16* (CAA→TAA)	yahE	DUF2877 family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (43/32);  total (43/32)

GGAACGGGGCTTTGAAGAGATTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGCAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGCACTCACTCACTTTG  >  NC_000913/335834‑336102                                                                                                                                         |                                                                                                                                     ggAACGGGGCTTTGAAGAGATTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTaa                                                                                                                                    >  2:1126078/1‑139 (MQ=255) ggAACGGGGCTTTGAAGAGATTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTaa                                                                                                                                    >  2:319261/1‑139 (MQ=255) ggAACGGGGCTTTGAAGAGATTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTaa                                                                                                                                    >  2:832402/1‑139 (MQ=255)              gAAGAGATTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAAc                                                                                                                       <  2:1017482/139‑1 (MQ=255)              gAAGAGATTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAAc                                                                                                                       <  2:1616116/139‑1 (MQ=255)              gAAGAGATTGCCCAGTTACTGATTGTCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGTCGCAGCGGGGGTAAGGACAGGTTGAAc                                                                                                                       <  2:1161950/139‑1 (MQ=255)                  agaTTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGt                                                                                                                   <  1:623357/139‑1 (MQ=255)                      tGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCttt                                                                                                               >  1:1766444/1‑139 (MQ=255)                      tGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGAGCAGCGGGGGTAAGGACAGGTTGAACAGGTCttt                                                                                                               >  1:604423/1‑139 (MQ=255)                       gcgcAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTg                                                                                                              >  1:326267/1‑139 (MQ=255)                        cgcAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAgg                                                                                                                                  <  1:444236/118‑1 (MQ=255)                        cgcAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAgg                                                                                                                                  >  2:444236/1‑118 (MQ=255)                          cAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCa                                                                                                           <  2:74871/139‑1 (MQ=255)                            gTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACGGGTTGAACAGGTCTTTGCCaga                                                                                                         <  1:1126078/139‑1 (MQ=255)                            gTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCaga                                                                                                         <  2:1559793/139‑1 (MQ=255)                                     tGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAat                                                                                                >  1:681782/1‑139 (MQ=255)                                       ccGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAatat                                                                                              >  1:202409/1‑139 (MQ=255)                                       ccGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAatat                                                                                              >  1:737686/1‑139 (MQ=255)                                          cAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCagag                                                                                                        >  2:611961/1‑126 (MQ=255)                                          cAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCagag                                                                                                        <  1:611961/126‑1 (MQ=255)                                             gTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACt                                                                                        >  2:396764/1‑139 (MQ=255)                                                        gggAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGcc                                                                             >  1:950656/1‑139 (MQ=255)                                                         ggAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCca                                                                            >  1:529969/1‑139 (MQ=255)                                                            ggCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCcagca                                                                         <  2:202409/139‑1 (MQ=255)                                                               ttGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAgg                                                                                                                                  >  1:1345098/1‑79 (MQ=255)                                                               ttGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAgg                                                                                                                                  <  2:1345098/79‑1 (MQ=255)                                                               ttGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGtt                                                                      >  1:520914/1‑139 (MQ=255)                                                                     tCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCa                                                                                                           <  1:710241/96‑1 (MQ=255)                                                                     tCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCa                                                                                                           >  2:710241/1‑96 (MQ=255)                                                                      cAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACg                                                               <  2:326267/139‑1 (MQ=255)                                                                      cAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACg                                                               <  2:737686/139‑1 (MQ=255)                                                                      cAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACg                                                               <  2:950656/139‑1 (MQ=255)                                                                       aaCACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACgc                                                              >  2:1041376/1‑139 (MQ=255)                                                                        acacacCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACgct                                                             <  2:520914/139‑1 (MQ=255)                                                                            acCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAAc                                                                                                                       >  1:924964/1‑77 (MQ=255)                                                                            acCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAAc                                                                                                                       <  2:924964/77‑1 (MQ=255)                                                                            acCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCtt                                                         >  2:775096/1‑139 (MQ=255)                                                                                 aaaGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTaaa                                                                                                 >  1:50846/1‑94 (MQ=255)                                                                                 aaaGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTaaa                                                                                                 <  2:50846/94‑1 (MQ=255)                                                                                 aaaGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCGGCAGTTGCTGACGCTGCTTTGTGa                                                    >  2:1630783/1‑139 (MQ=255)                                                                                 aaaGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGa                                                    >  1:1048003/1‑139 (MQ=255)                                                                                 aaaGACGGTGTGTGTGGGCGCTAACTGCAGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGcc                                                                             <  1:391078/114‑1 (MQ=255)                                                                                 aaaGACGGTGTGTGTGGGCGCTAACTGCAGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGcc                                                                             >  2:391078/1‑114 (MQ=255)                                                                                      cGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGt                                               >  1:126351/1‑139 (MQ=255)                                                                                        gtgtgtgtGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACg                                                               <  2:799627/121‑1 (MQ=255)                                                                                        gtgtgtgtGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACg                                                               >  1:799627/1‑121 (MQ=255)                                                                                        gtgtgtgtGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTAc                                             >  2:948058/1‑139 (MQ=255)                                                                                           tgtgtgGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGAc                                          >  1:487327/1‑139 (MQ=255)                                                                                                 ggCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAgg                                                                                                                                  <  1:268178/45‑1 (MQ=255)                                                                                                 ggCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAgg                                                                                                                                  >  2:268178/1‑45 (MQ=255)                                                                                                  gcgcTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATgcgc                                   >  1:1626734/1‑139 (MQ=255)                                                                                                   cgcTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTtgtg                                                     <  1:380091/120‑1 (MQ=255)                                                                                                   cgcTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTtgtg                                                     >  2:380091/1‑120 (MQ=255)                                                                                                   cgcTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTtgcg                                                     >  2:725220/1‑118 (MQ=255)                                                                                                   cgcTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTtgcg                                                     <  1:725220/120‑3 (MQ=255)                                                                                                     cTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCaa                                >  1:1033070/1‑139 (MQ=255)                                                                                                       aaCTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAAc                              >  2:1063342/1‑139 (MQ=255)                                                                                                        aCTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCg                                                                                            <  1:734840/76‑1 (MQ=255)                                                                                                        aCTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCg                                                                                            >  2:734840/1‑76 (MQ=255)                                                                                                        aCTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACa                             >  2:1178872/1‑139 (MQ=255)                                                                                                          tgcggatgcggaTTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGt                           >  1:1125993/1‑139 (MQ=255)                                                                                                           gcggatgcggatTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACgctgct                                                          <  2:1131649/107‑1 (MQ=255)                                                                                                           gcggatgcggatTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACgctgct                                                          >  1:1131649/1‑107 (MQ=255)                                                                                                            cggatgcggatTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTg                         >  2:1779088/1‑139 (MQ=255)                                                                                                            cggatgcggatTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTg                         >  1:739968/1‑139 (MQ=255)                                                                                                                  cggatTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCa                                                                                                           <  2:1048551/51‑1 (MQ=255)                                                                                                                  cggatTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCa                                                                                                           >  1:1048551/1‑51 (MQ=255)                                                                                                                      ttttCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGca               <  2:739968/139‑1 (MQ=255)                                                                                                                      ttttCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGca               <  2:5986/139‑1 (MQ=255)                                                                                                                         tCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGcactc            >  2:1088308/1‑139 (MQ=255)                                                                                                                           tGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGcactcac          <  1:948058/139‑1 (MQ=255)                                                                                                                              cgcAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGcactcactca       <  2:673088/139‑1 (MQ=255)                                                                                                                              cgcAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGcactcactca       <  2:1626734/139‑1 (MQ=255)                                                                                                                               gcAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGcactcactcac      <  1:1178872/139‑1 (MQ=255)                                                                                                                                   cGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGCACTCACTCACTTtg  >  1:1486468/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                     GGAACGGGGCTTTGAAGAGATTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGCAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGCACTCACTCACTTTG  >  NC_000913/335834‑336102

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