Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,398,875	C→T	M1M (GTG→ATG) † *368* (TAG→TAA)	mreD ← mreC ←	cell wall structural complex MreBCD transmembrane component MreD cell wall structural complex MreBCD transmembrane component MreC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,398,875	0	C	T	100.0%	240.0 / NA	73	M1M (GTG→ATG) † *368* (TAG→TAA)	mreD mreC	cell wall structural complex MreBCD transmembrane component MreD cell wall structural complex MreBCD transmembrane component MreC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (41/32);  total (41/32)

AATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCT  >  NC_000913/3398740‑3399008                                                                                                                                        |                                                                                                                                      aaTACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTAt                                                                                                                                    >  1:797910/1‑139 (MQ=255)          aGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCACCATTATTGCCCTccc                                                                                                                           <  1:1586112/139‑1 (MQ=255)                gTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGTAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCa                                                                                                                     <  1:964333/139‑1 (MQ=255)                gTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCa                                                                                                                     <  2:797910/139‑1 (MQ=255)                  aaCACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCAcg                                                                                                                   >  2:1029779/1‑139 (MQ=255)                  aaCACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCAcg                                                                                                                   >  2:1415363/1‑139 (MQ=255)                     acCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACgcgc                                                                                                                >  1:611658/1‑139 (MQ=255)                         aGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCGCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGc                                                                                                            >  1:1599991/1‑139 (MQ=255)                         aGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGc                                                                                                            >  2:1567727/1‑139 (MQ=255)                         aGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAA‑GAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGcc                                                                                                                                <  1:705636/118‑1 (MQ=255)                         aGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAA‑GAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGcc                                                                                                                                >  2:705636/1‑118 (MQ=255)                                    aaaaCAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAg                                                                                                              >  2:484265/1‑126 (MQ=255)                                    aaaaCAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAg                                                                                                              <  1:484265/126‑1 (MQ=255)                                    aaaaCAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGcacc                                                                                                 >  1:522516/1‑139 (MQ=255)                                             aGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTgcgg                                                                                        <  2:596788/139‑1 (MQ=255)                                              ggTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTgcggc                                                                                       >  2:76751/1‑139 (MQ=255)                                                  aTCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCAGCGCACGCGCAGGCGGTTGAGCACCACTTTgcggcggcg                                                                                   >  1:454071/1‑139 (MQ=255)                                                    ccGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTa                                                                                 <  1:1576526/139‑1 (MQ=255)                                                    ccGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTa                                                                                 <  1:1326488/139‑1 (MQ=255)                                                      ggCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGc                                                                               <  2:522516/139‑1 (MQ=255)                                                        ccAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGccc                                                                             >  2:783769/1‑139 (MQ=255)                                                         cAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCt                                                                            <  2:1200804/139‑1 (MQ=255)                                                           gggCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCttt                                                                          <  1:783769/139‑1 (MQ=255)                                                            ggCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTgcggcggcg                                                                                   <  2:294133/129‑1 (MQ=255)                                                            ggCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTgcggcggcg                                                                                   >  1:294133/1‑129 (MQ=255)                                                                            aCAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTagca                                                         >  1:688050/1‑139 (MQ=255)                                                                                  cAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGc                                                   <  1:1415363/139‑1 (MQ=255)                                                                                     tGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAgg                                                                                                             >  1:358188/1‑78 (MQ=255)                                                                                     tGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAgg                                                                                                             <  2:358188/78‑1 (MQ=255)                                                                                      gAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCtgtg                                               <  2:1599991/139‑1 (MQ=255)                                                                                          aaaGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTgcgg                                                                                        >  2:1415224/1‑94 (MQ=255)                                                                                          aaaGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTgcgg                                                                                        <  1:1415224/94‑1 (MQ=255)                                                                                            agagagCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGgcg                                         <  2:611658/139‑1 (MQ=255)                                                                                              agagCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCgca                                       <  1:1567727/139‑1 (MQ=255)                                                                                               gagCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCAcgcg                                                                                                                 <  2:38308/64‑1 (MQ=255)                                                                                               gagCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCAcgcg                                                                                                                 >  1:38308/1‑64 (MQ=255)                                                                                                 gCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTgcgg                                                                                        <  2:1516835/87‑1 (MQ=255)                                                                                                 gCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTgcgg                                                                                        >  1:1516835/1‑87 (MQ=255)                                                                                                 gCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCgcagca                                    >  2:978728/1‑139 (MQ=255)                                                                                                 gCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCgcagca                                    >  1:481791/1‑139 (MQ=255)                                                                                                  ccAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCgg                                                                                                          >  1:1050602/1‑68 (MQ=255)                                                                                                  ccAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCgg                                                                                                          <  2:1050602/68‑1 (MQ=255)                                                                                                    aGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCa                                                                                                               <  1:1557054/61‑1 (MQ=255)                                                                                                    aGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCa                                                                                                               >  2:1557054/1‑61 (MQ=255)                                                                                                         aCCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGt                            <  1:358576/139‑1 (MQ=255)                                                                                                          cccAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTa                           >  2:767737/1‑139 (MQ=255)                                                                                                               cGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGc                      >  2:598267/1‑139 (MQ=255)                                                                                                               cGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGc                      >  2:472775/1‑139 (MQ=255)                                                                                                                gTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGgcg                                >  2:819195/1‑128 (MQ=255)                                                                                                                gTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGgcg                                <  1:819195/128‑1 (MQ=255)                                                                                                                  cccTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGc                         >  2:1166781/1‑133 (MQ=255)                                                                                                                  cccTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGc                         <  1:1166781/133‑1 (MQ=255)                                                                                                                    cTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCgcagc                                     >  1:639719/1‑119 (MQ=255)                                                                                                                    cTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCgcagc                                     <  2:639719/119‑1 (MQ=255)                                                                                                                     tGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCATAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTcc                >  2:82488/1‑139 (MQ=255)                                                                                                                           cGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGc          >  1:1268074/1‑139 (MQ=255)                                                                                                                            gATAGCTCGCCATTATTGCCCTCCCGTCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGcc         >  1:1330446/1‑139 (MQ=255)                                                                                                                            gATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGga                                                                     >  2:232035/1‑79 (MQ=255)                                                                                                                            gATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGga                                                                     <  1:232035/79‑1 (MQ=255)                                                                                                                            gATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGcc         <  2:543964/139‑1 (MQ=255)                                                                                                                            gATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGcc         <  2:1472213/139‑1 (MQ=255)                                                                                                                                 gtCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtg    >  1:1123942/2‑139 (MQ=255)                                                                                                                                 cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGcc         >  1:447896/1‑134 (MQ=255)                                                                                                                                 cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGcc         <  2:447896/134‑1 (MQ=255)                                                                                                                                 cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtg    >  2:1629414/1‑139 (MQ=255)                                                                                                                                 cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtg    >  2:1125825/1‑139 (MQ=255)                                                                                                                                 cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtg    >  2:1261128/1‑139 (MQ=255)                                                                                                                                 cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtg    >  1:672277/1‑139 (MQ=255)                                                                                                                                 cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtg    >  1:665069/1‑139 (MQ=255)                                                                                                                                 cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCCGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtg    >  1:861636/1‑139 (MQ=255)                                                                                                                                  tCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGgcg                                         <  1:1482086/101‑1 (MQ=255)                                                                                                                                  tCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGgcg                                         >  2:1482086/1‑101 (MQ=255)                                                                                                                                   cGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtgct  <  1:1629414/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                      AATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCT  >  NC_000913/3398740‑3399008

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