Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	486,641	T→C	S36P (TCA→CCA)	mscK →	mechanosensitive channel protein, intermediate conductance, K+ regulated

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	486,641	0	T	C	100.0%	298.6 / NA	82	S36P (TCA→CCA)	mscK	mechanosensitive channel protein, intermediate conductance, K+ regulated
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (46/36);  total (46/36)

TTCTTCAGGTTCAGAAACCTTCATTCATCATGACTATGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGATACATTAGCCACC  >  NC_000913/486507‑486766                                                                                                                                       |                                                                                                                              ttcttcAGGTTCAGAAACCTTCATTCATCATGACTATGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATc                                                                                                                           >  1:56792/1‑139 (MQ=255)    ttcAGGTTCAGAAACCTTCATTCATCATGACTATGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGaa                                                                                                                        <  2:834702/139‑1 (MQ=255)         gTTCAGAAACCTTCATTCATCATGACTATGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTg                                                                                                                   <  1:48174/139‑1 (MQ=255)            cAGAAACCTTCATTCATCATGACTATGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATc                                                                                                                           <  2:1780422/128‑1 (MQ=255)            cAGAAACCTTCATTCATCATGACTATGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATc                                                                                                                           >  1:1780422/1‑128 (MQ=255)            cAGAAACCTTCATTCATCATGACTATGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATc                                                                                                                >  1:665388/1‑139 (MQ=255)                    ttcattcaTCATGACTATGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACa                                                                                                        >  2:465466/1‑139 (MQ=255)                      cattcatCATAACTATGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTC‑TGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACaaaa                                                                                                     >  1:659265/1‑139 (MQ=255)                             catGACTATGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACTCGGCGTTTGCGCGGTCGCCATCGAATGGTGATCTGCCGACAAAAGCGGAc                                                                                               >  2:897966/1‑139 (MQ=255)                             catGACTATGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGAc                                                                                                         <  2:142832/129‑1 (MQ=255)                             catGACTATGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGAc                                                                                                         >  1:142832/1‑129 (MQ=255)                             catGACTATGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCgg                                                                                                 >  2:297183/1‑137 (MQ=255)                             catGACTATGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCgg                                                                                                 <  1:297183/137‑1 (MQ=255)                                    aTGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAgg                                                                                        <  2:56792/139‑1 (MQ=255)                                     tGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAg                                                                                                    >  2:469677/1‑126 (MQ=255)                                     tGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAg                                                                                                    <  1:469677/126‑1 (MQ=255)                                     tGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGgc                                                                                       <  2:665388/139‑1 (MQ=255)                                     tGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGgc                                                                                       >  2:492353/1‑139 (MQ=255)                                     tGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGgc                                                                                       >  1:1068705/1‑139 (MQ=255)                                            gTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCaa                                                                                   >  1:1625660/1‑136 (MQ=255)                                            gTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCaa                                                                                   <  2:1625660/136‑1 (MQ=255)                                            gTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACtt                                                                                >  2:202970/1‑139 (MQ=255)                                               ttACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACaaaa                                                                                                     >  1:1495344/1‑115 (MQ=255)                                               ttACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACaaaa                                                                                                     <  2:1495344/115‑1 (MQ=255)                                                tACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACt                                                                            >  1:453398/1‑139 (MQ=255)                                                 aCAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTc                                                                           >  1:439069/1‑139 (MQ=255)                                                       gATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGAt                                                                                                                 <  1:713718/95‑1 (MQ=255)                                                       gATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGAt                                                                                                                 >  2:713718/1‑95 (MQ=255)                                                       gATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACtaaa                                                                     >  2:1318729/1‑139 (MQ=255)                                                       gATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACtaaa                                                                     >  2:1271740/1‑139 (MQ=255)                                                        aTCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACtaaat                                                                    >  2:591942/1‑139 (MQ=255)                                                         tCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGaa                                                                                                                        >  2:192360/1‑86 (MQ=255)                                                         tCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGaa                                                                                                                        <  1:192360/86‑1 (MQ=255)                                                         tCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTg                                                                               >  2:1172908/1‑127 (MQ=255)                                                         tCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTg                                                                               <  1:1172908/127‑1 (MQ=255)                                                         tCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACtaaata                                                                   >  2:26712/1‑139 (MQ=255)                                                          cACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACtaaataa                                                                  <  1:465466/139‑1 (MQ=255)                                                             ggCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACa                                                               >  2:898871/1‑139 (MQ=255)                                                                  tttGTTTTTTCGGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAgg                                                                                        <  2:1035922/109‑1 (MQ=255)                                                                  tttGTTTTTTCGGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAgg                                                                                        >  1:1035922/1‑109 (MQ=255)                                                                    tGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCa                                                                                                                             <  1:653510/70‑1 (MQ=255)                                                                    tGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCa                                                                                                                             >  2:653510/1‑70 (MQ=255)                                                                        ttttCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCttt                                                    <  2:526299/139‑1 (MQ=255)                                                                             aGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCt                                               <  1:1063158/139‑1 (MQ=255)                                                                                  ttATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGa                                          <  1:911535/139‑1 (MQ=255)                                                                                               ttttgtCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGc                                                                                                   >  2:954479/1‑69 (MQ=255)                                                                                               ttttgtCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGc                                                                                                   <  1:954479/69‑1 (MQ=255)                                                                                                tttgtCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACaa                                                              >  1:1008379/1‑105 (MQ=255)                                                                                                tttgtCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACaa                                                              <  2:1008379/105‑1 (MQ=255)                                                                                                   tgtctGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAgg                         <  2:659265/135‑1 (MQ=255)                                                                                                   gtCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAgg                         >  2:422892/1‑139 (MQ=255)                                                                                                     cTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGAt                       <  2:453398/139‑1 (MQ=255)                                                                                                        gTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTg                    >  1:549413/1‑139 (MQ=255)                                                                                                        gTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTg                    >  1:233063/1‑139 (MQ=255)                                                                                                         ttATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGa                   <  1:591942/139‑1 (MQ=255)                                                                                                            tGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAg                >  2:1756860/1‑139 (MQ=255)                                                                                                             gCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAGGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGa               <  1:1271740/139‑1 (MQ=255)                                                                                                              ccAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGAt              >  2:891175/1‑139 (MQ=255)                                                                                                              ccAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGAt              >  1:134377/1‑139 (MQ=255)                                                                                                                 gAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGATACa           <  1:26712/139‑1 (MQ=255)                                                                                                                  aaCACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCt                                               <  2:640248/102‑1 (MQ=255)                                                                                                                  aaCACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCt                                               >  1:640248/1‑102 (MQ=255)                                                                                                                  aaCACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGATACAt          >  1:1219156/1‑139 (MQ=255)                                                                                                                  aaCACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGATACAt          >  1:1664315/1‑139 (MQ=255)                                                                                                                   acacGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTc                                                                           <  1:554112/73‑1 (MQ=255)                                                                                                                   acacGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTc                                                                           >  2:554112/1‑73 (MQ=255)                                                                                                                   acacGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCt                     <  2:1025615/127‑1 (MQ=255)                                                                                                                   acacGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCt                     >  1:1025615/1‑127 (MQ=255)                                                                                                                      cGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGa                   <  2:1697594/126‑1 (MQ=255)                                                                                                                      cGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGa                   >  1:1697594/1‑126 (MQ=255)                                                                                                                       ggCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACtaaata                                                                   >  1:932975/1‑77 (MQ=255)                                                                                                                       ggCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACtaaata                                                                   <  2:932975/77‑1 (MQ=255)                                                                                                                       ggCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCttt                                                    >  1:71423/1‑92 (MQ=255)                                                                                                                       ggCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCttt                                                    <  2:71423/92‑1 (MQ=255)                                                                                                                       ggCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGATACATTAGcc     >  2:728408/1‑139 (MQ=255)                                                                                                                       ggCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGATACATTAGcc     >  1:131549/1‑139 (MQ=255)                                                                                                                         cGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTAACAGATACATTAGCCAc   >  2:749538/1‑139 (MQ=255)                                                                                                                          gTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGATACATTAGCCAcc  <  2:233063/139‑1 (MQ=255)                                                                                                                          gTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGATACATTAGCCAcc  <  2:1664315/139‑1 (MQ=255)                                                                                                                          gTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGATACATTAGCCAcc  <  1:1756860/139‑1 (MQ=255)                                                                                                                                   ggCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCt                                               <  2:1180641/85‑1 (MQ=255)                                                                                                                                   ggCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCt                                               >  1:1180641/1‑85 (MQ=255)                                                                                                                                       |                                                                                                                              TTCTTCAGGTTCAGAAACCTTCATTCATCATGACTATGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGATACATTAGCCACC  >  NC_000913/486507‑486766

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