Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	664,102	C→T	*363* (TAG→TAA)	rlpA ←	septal ring protein, suppressor of prc, minor lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	664,102	0	C	T	100.0%	254.6 / NA	78	*363* (TAG→TAA)	rlpA	septal ring protein, suppressor of prc, minor lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (37/41);  total (37/41)

CTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGA  >  NC_000913/663967‑664234                                                                                                                                        |                                                                                                                                     cTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTAc                                                                                                                                   >  2:446566/1‑139 (MQ=255)   gAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTg                                                                                                                                 <  1:446566/139‑1 (MQ=255)    aaCTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTgc                                                                                                                                >  2:1425778/1‑139 (MQ=255)    aaCTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTgc                                                                                                                                >  2:619035/1‑139 (MQ=255)     aCTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTgcg                                                                                                                               <  1:847101/139‑1 (MQ=255)     aCTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTgcg                                                                                                                               <  2:1283164/139‑1 (MQ=255)     aCTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTgcg                                                                                                                               <  2:442157/139‑1 (MQ=255)      cTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTgcgc                                                                                                                              <  2:1274668/139‑1 (MQ=255)        aCGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACTGGCAGGCATCTGACTTTCCGCGTGAGTTTGTTAACGTCATTTACAGAAATCGACACATCGGATGCCTGCTTTACGTTACTGCGCgg                                                                                                                            <  1:682160/139‑1 (MQ=255)         cGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTgc                                                                                                                                <  2:1221163/134‑1 (MQ=255)         cGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTgc                                                                                                                                >  1:1221163/1‑134 (MQ=255)         cGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGgt                                                                                                                           <  1:619035/139‑1 (MQ=255)          gACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGgta                                                                                                                          >  2:53005/1‑139 (MQ=255)           aCATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGgtag                                                                                                                         <  2:152240/139‑1 (MQ=255)           aCATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGgtag                                                                                                                         <  1:1425778/139‑1 (MQ=255)                         aTTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTgcg                                                                                                                               >  2:638581/1‑119 (MQ=255)                         aTTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTgcg                                                                                                                               <  1:638581/119‑1 (MQ=255)                          ttAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTa                                                                                                          >  2:395705/1‑139 (MQ=255)                                        tactaTAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTgcgc                                                                                                                              <  1:716825/105‑1 (MQ=255)                                        tactaTAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTgcgc                                                                                                                              >  2:716825/1‑105 (MQ=255)                                               gCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGgtag                                                                                                                         <  1:219469/103‑1 (MQ=255)                                               gCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGgtag                                                                                                                         >  2:219469/1‑103 (MQ=255)                                                   aTGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTtg                                                                                 <  1:554421/139‑1 (MQ=255)                                                    tGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTgcg                                                                                                                               >  2:42577/1‑92 (MQ=255)                                                    tGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTgcg                                                                                                                               <  1:42577/92‑1 (MQ=255)                                                    tGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTgc                                                                                                                                >  1:608186/1‑91 (MQ=255)                                                    tGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTgc                                                                                                                                <  2:608186/91‑1 (MQ=255)                                                    tGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCgg                                                                                                                            <  1:763440/95‑1 (MQ=255)                                                    tGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCgg                                                                                                                            >  2:763440/1‑95 (MQ=255)                                                    tGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAAt                                                                                                                 >  2:237279/1‑106 (MQ=255)                                                    tGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAAt                                                                                                                 <  1:237279/106‑1 (MQ=255)                                                    tGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGAc                                                                                                              >  1:674267/1‑109 (MQ=255)                                                    tGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGAc                                                                                                              <  2:674267/109‑1 (MQ=255)                                                    tGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTtgc                                                                                >  2:716257/1‑139 (MQ=255)                                                    tGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTtgc                                                                                >  2:1302488/1‑139 (MQ=255)                                                    tGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTtgc                                                                                >  1:1120958/1‑139 (MQ=255)                                                    tGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTtgc                                                                                >  1:111565/1‑139 (MQ=255)                                                    tGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTtg                                                                                 >  1:978971/1‑138 (MQ=255)                                                    tGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTtg                                                                                 <  2:978971/138‑1 (MQ=255)                                                    tGCACTACCGACAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTtg                                                                                 <  2:919607/138‑1 (MQ=255)                                                    tGCACTACCGACAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTtg                                                                                 >  1:919607/1‑138 (MQ=255)                                                     gCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGgtagta                                                                                                                       >  1:1106029/1‑99 (MQ=255)                                                     gCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGgtagta                                                                                                                       <  2:1106029/99‑1 (MQ=255)                                                       aCTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGa                                                                                                               >  2:811850/1‑105 (MQ=255)                                                       aCTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGa                                                                                                               <  1:811850/105‑1 (MQ=255)                                                        cTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTTGCTGCa                                                                            <  1:650601/139‑1 (MQ=255)                                                         tACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACg                                                                                    <  1:739086/130‑1 (MQ=255)                                                         tACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACg                                                                                    >  2:739086/1‑130 (MQ=255)                                                         tACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGt                                                                                   <  1:1604802/131‑1 (MQ=255)                                                         tACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGt                                                                                   >  2:1604802/1‑131 (MQ=255)                                                          aCCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCgg                                                                                              >  1:782297/1‑119 (MQ=255)                                                          aCCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCgg                                                                                              <  2:782297/119‑1 (MQ=255)                                                          aCCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCg                                                                                               >  1:748203/1‑118 (MQ=255)                                                          aCCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCg                                                                                               <  2:748203/118‑1 (MQ=255)                                                             ggcaggcaTCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTTGCTGCAAGGTa                                                                       <  1:1302488/139‑1 (MQ=255)                                                                aggcaTCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTTGCTGCAAGGTACTg                                                                    >  1:444725/1‑139 (MQ=255)                                                                     tCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGgtag                                                                                                                         <  2:1090751/81‑1 (MQ=255)                                                                     tCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGgtag                                                                                                                         >  1:1090751/1‑81 (MQ=255)                                                                             tCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTg                                                                                                                                 <  2:1467615/65‑1 (MQ=255)                                                                             tCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTg                                                                                                                                 >  1:1467615/1‑65 (MQ=255)                                                                             tCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGc                                                       >  1:697671/1‑139 (MQ=255)                                                                               cgcgTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGtaataa                                                                                                                   <  2:1335725/77‑1 (MQ=255)                                                                               cgcgTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGtaataa                                                                                                                   >  1:1335725/1‑77 (MQ=255)                                                                                    gACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTTGCTGCa                                                                            <  2:142862/111‑1 (MQ=255)                                                                                    gACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTTGCTGCa                                                                            >  1:142862/1‑111 (MQ=255)                                                                                        ttGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATggg                                            <  1:716257/139‑1 (MQ=255)                                                                                        ttGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATggg                                            <  1:395705/139‑1 (MQ=255)                                                                                             aaCGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAg                                       <  2:1120958/139‑1 (MQ=255)                                                                                                 tCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTgg                                   <  2:111565/139‑1 (MQ=255)                                                                                                    tttACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGgtag                                                                                                                         <  1:1452985/50‑1 (MQ=255)                                                                                                    tttACAGAAATTGACACATCAGATGCCTGCTTTACGTTACTGCGCGgtag                                                                                                                         >  2:1452985/1‑50 (MQ=255)                                                                                                             aTTGACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGAcc                       >  1:440078/1‑139 (MQ=255)                                                                                                                gACACATCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACcgcg                    <  2:440078/139‑1 (MQ=255)                                                                                                                     aTCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAg                                       >  1:1485990/1‑115 (MQ=255)                                                                                                                     aTCAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAg                                       <  2:1485990/115‑1 (MQ=255)                                                                                                                       cAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTTGCTGCAAGGTACTGGCTTCGGCt                                                           <  1:148801/93‑1 (MQ=255)                                                                                                                       cAGATGCCTGCTTTACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTTGCTGCAAGGTACTGGCTTCGGCt                                                           >  2:148801/1‑93 (MQ=255)                                                                                                                                  tttACGTTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGa  >  1:245685/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                     CTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGA  >  NC_000913/663967‑664234

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