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%	252.0 / NA	77	Q16* (CAA→TAA)	yahE	DUF2877 family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (40/37);  total (40/37)

TTTGAAGAGATTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGCAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGCACTCACTCACTTTG  >  NC_000913/335844‑336102                                                                                                                                |                                                                                                                                     tttGAAGAGATTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTg                                                                                                                          <  1:840464/139‑1 (MQ=255)    gAAGAGATTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAAc                                                                                                                       <  1:498027/139‑1 (MQ=255)    gAAGAGATTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAAc                                                                                                                       <  2:538664/139‑1 (MQ=255)          aTTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCt                                                                                                                 >  2:463985/1‑139 (MQ=255)          aTTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCt                                                                                                                 >  1:296918/1‑139 (MQ=255)                    tACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGc                                                                                                       >  2:981912/1‑139 (MQ=255)                           tGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCaga                                                                                                         >  1:157106/1‑130 (MQ=255)                           tGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCaga                                                                                                         <  2:157106/130‑1 (MQ=255)                             ccGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCTGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAatat                                                                                              >  2:107149/1‑139 (MQ=255)                             ccGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAgg                                                                                                                             >  2:1007708/1‑108 (MQ=255)                             ccGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAgg                                                                                                                             <  1:1007708/108‑1 (MQ=255)                             ccGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAatat                                                                                              >  2:1177636/1‑139 (MQ=255)                             ccGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAatat                                                                                              >  2:278803/1‑139 (MQ=255)                             ccGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAatat                                                                                              >  2:1030028/1‑139 (MQ=255)                               gCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCg                                                                                            <  1:278803/139‑1 (MQ=255)                                   gTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGc                                                                                                             <  2:763937/118‑1 (MQ=255)                                   gTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGc                                                                                                             >  1:763937/1‑118 (MQ=255)                                          aaCCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCt                                                                                 <  2:296918/139‑1 (MQ=255)                                              gggAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGcc                                                                             <  2:1004989/139‑1 (MQ=255)                                               ggAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCACTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAAc                                                                                                                       <  1:985625/96‑1 (MQ=255)                                               ggAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCACTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAAc                                                                                                                       >  2:985625/1‑96 (MQ=255)                                                gAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGc                                                                              <  2:10464/136‑1 (MQ=255)                                                gAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGc                                                                              >  1:10464/1‑136 (MQ=255)                                                gAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATATTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCcag                                                                           >  1:286446/1‑139 (MQ=255)                                                       gCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCg                                                                                            <  2:357796/115‑1 (MQ=255)                                                       gCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCg                                                                                            >  1:357796/1‑115 (MQ=255)                                                           tCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGAc                                                                >  1:1009431/1‑139 (MQ=255)                                                                 cacCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCACAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACgctgct                                                          >  2:43325/1‑139 (MQ=255)                                                                  acCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCtt                                                         >  2:700120/1‑139 (MQ=255)                                                                  acCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCtt                                                         >  1:257029/1‑139 (MQ=255)                                                                  acCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCtt                                                         >  2:879170/1‑139 (MQ=255)                                                                  acCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATCTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCtt                                                         >  1:237427/1‑139 (MQ=255)                                                                    cAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTtg                                                       >  1:988323/1‑139 (MQ=255)                                                                      gAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGc                                                                    <  2:42095/124‑1 (MQ=255)                                                                      gAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGc                                                                    >  1:42095/1‑124 (MQ=255)                                                                      gAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTtgtg                                                     <  2:1009431/139‑1 (MQ=255)                                                                      gAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTtgtg                                                     <  1:1012665/139‑1 (MQ=255)                                                                      gAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTtgtg                                                     <  1:700120/139‑1 (MQ=255)                                                                        aaGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGaa                                                   <  2:988323/139‑1 (MQ=255)                                                                         aGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACa                                                                                                                               <  2:687296/62‑1 (MQ=255)                                                                         aGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACa                                                                                                                               >  1:687296/1‑62 (MQ=255)                                                                            cGGTGTGTGTGGGCGCTAACTGCGGATGTGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCaga                                                                                                         <  2:822043/81‑1 (MQ=255)                                                                            cGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCaga                                                                                                         >  1:822043/1‑81 (MQ=255)                                                                            cGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCa                                                                                                           >  2:841500/1‑79 (MQ=255)                                                                            cGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCa                                                                                                           <  1:841500/79‑1 (MQ=255)                                                                            cGGTGTGGGTGGGCGCGAACTGCGGATGCGGATTTTCTGCCGCATC‑GGTGGTAAGGACAGGTTGAACATGTCGTTGCCAGACCGGTATAGATAGAACTCCCGGCACGCTAGCAGTTACTGACGCTGCttctttcagact                                               >  1:327965/1‑129 (MQ=255)                                                                              gtgtgtgtGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTa                                              >  1:773702/1‑139 (MQ=255)                                                                              gtgtgtgtGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCa                                                                                                           >  1:344591/1‑77 (MQ=255)                                                                              gtgtgtgtGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCa                                                                                                           <  2:344591/77‑1 (MQ=255)                                                                              gtgtgtgtGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTAc                                             >  2:130162/1‑139 (MQ=255)                                                                                       ggCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATgcg                                    >  2:1149942/1‑139 (MQ=255)                                                                                         cgcTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAgg                                                                                                                                  >  1:680426/1‑43 (MQ=255)                                                                                         cgcTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAgg                                                                                                                                  <  2:680426/43‑1 (MQ=255)                                                                                         cgcTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCa                                                                                                           >  1:459478/1‑66 (MQ=255)                                                                                         cgcTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCa                                                                                                           <  2:459478/66‑1 (MQ=255)                                                                                         cgcTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGcc                                  >  1:972855/1‑139 (MQ=255)                                                                                                  cggatgcggatTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAAc                                                                                                                       >  2:414468/1‑45 (MQ=255)                                                                                                  cggatgcggatTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAAc                                                                                                                       <  1:414468/45‑1 (MQ=255)                                                                                                  cggatgcggatTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCaaa                               >  2:704127/1‑133 (MQ=255)                                                                                                  cggatgcggatTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTACTTTGTGAAGAGTACGACAATGCGCCaaa                               <  1:704127/133‑1 (MQ=255)                                                                                                      tgcggatTTTCTGGCGCAGC‑GGGGGTATGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACt                                                                                        >  1:530590/1‑72 (MQ=255)                                                                                                      tgcggatTTTCTGGCGCAGC‑GGGGGTATGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACt                                                                                        <  2:530590/72‑1 (MQ=255)                                                                                                      tgcggatTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAAc                                                                                                                       >  2:918206/1‑41 (MQ=255)                                                                                                      tgcggatTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAAc                                                                                                                       <  1:918206/41‑1 (MQ=255)                                                                                                      tgcggatTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCgg                     <  1:130162/139‑1 (MQ=255)                                                                                                          gatTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTgg                 <  1:1030028/139‑1 (MQ=255)                                                                                                          gatTTTCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTgg                 <  1:755604/139‑1 (MQ=255)                                                                                                            ttttCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCaaa                               <  1:920297/123‑1 (MQ=255)                                                                                                            ttttCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCaaa                               >  2:920297/1‑123 (MQ=255)                                                                                                            ttttCTGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGca               <  1:879170/139‑1 (MQ=255)                                                                                                                 tGGCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGcactcac          <  1:1149942/139‑1 (MQ=255)                                                                                                                  ggCGCAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGcactcact         <  1:463985/139‑1 (MQ=255)                                                                                                                  ggCACAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGcactcact         <  1:43325/139‑1 (MQ=255)                                                                                                                      cAGC‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGcactcactcact     <  1:1177636/139‑1 (MQ=255)                                                                                                                         c‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCaga                                                                                                         <  2:93347/36‑1 (MQ=255)                                                                                                                         c‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCaga                                                                                                         >  1:93347/1‑36 (MQ=255)                                                                                                                         c‑GGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGCACTCACTCACTTtg  >  1:887600/1‑139 (MQ=255)                                                                                                                                |                                                                                                                                     TTTGAAGAGATTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGC‑GGGGGCAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGCACTCACTCACTTTG  >  NC_000913/335844‑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: ‑