Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,347,323	A→G	S2P (TCA→CCA)	fumB ←	anaerobic class I fumarate hydratase (fumarase B)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,347,323	0	A	G	100.0%	211.0 / NA	60	S2P (TCA→CCA)	fumB	anaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/33);  total (27/33)

CAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGCCGAAGACGCACGATACGCTCACACCAATCAAC  >  NC_000913/4347177‑4347470                                                                                                                                                   |                                                                                                                                                    cagcagGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATa                                                                                                                                                 <  1:656229/151‑1 (MQ=255) cagcagGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATa                                                                                                                                                 <  1:417521/151‑1 (MQ=255)  agcagGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAg                                                                                                                                                <  1:881597/151‑1 (MQ=255)  agcagGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAg                                                                                                                                                <  2:167479/151‑1 (MQ=255)   gcagGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGc                                                                                                                                               <  2:537725/151‑1 (MQ=255)             ggCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTg                                                                                                                                     >  1:842443/1‑151 (MQ=255)              gCTTCTGGTTCCACGTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTGCCCCATCGGGAAACGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCt                                                                                                                                      <  2:672226/149‑1 (MQ=255)                 tCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTa                                                                                                                                   <  1:238162/149‑1 (MQ=255)                    ggTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAAct                                                                                                                                <  1:778915/149‑1 (MQ=255)                    ggTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAAct                                                                                                                                <  2:262768/149‑1 (MQ=255)                       tCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGt                                                                                                                            <  2:180195/150‑1 (MQ=255)                         cACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATATAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTAtt                                                                                                                         <  2:95338/151‑1 (MQ=255)                         cACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGAGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTAtt                                                                                                                         >  2:246076/1‑151 (MQ=255)                             ttCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTg                                                                                                                       <  1:68058/149‑1 (MQ=255)                                   aTGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCg                                                                                                                >  1:429194/1‑150 (MQ=255)                                    tGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGa                                                                                                               >  2:226636/1‑150 (MQ=255)                                     ggTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGaaa                                                                                                             <  1:429518/151‑1 (MQ=255)                                              gTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGc                                                                                                    <  1:62425/151‑1 (MQ=255)                                              gTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGc                                                                                                    <  2:478250/151‑1 (MQ=255)                                                 gAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCcgc                                                                                                 >  1:566802/1‑151 (MQ=255)                                                 gAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCcgc                                                                                                 >  2:354627/1‑151 (MQ=255)                                                 gAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCcg                                                                                                  >  1:630136/1‑150 (MQ=255)                                                  aaGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCcgcg                                                                                                <  2:511448/151‑1 (MQ=255)                                                   aGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGt                                                                                               >  2:428190/1‑151 (MQ=255)                                                    gTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTa                                                                                              <  2:38496/151‑1 (MQ=255)                                                       ggCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTc                                                                                           <  1:356342/151‑1 (MQ=255)                                                            cGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACggg                                                                                      >  1:687442/1‑151 (MQ=255)                                                             gCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGggc                                                                                     <  1:379856/151‑1 (MQ=255)                                                                      aTCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCaa                                                                            <  1:871293/151‑1 (MQ=255)                                                                      aTCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCa                                                                             >  1:750083/1‑150 (MQ=255)                                                                            aGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGc                                                                       >  1:461058/1‑150 (MQ=255)                                                                               gAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCAttt                                                                   >  1:433852/1‑151 (MQ=255)                                                                                aggaGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGAGCAttt                                                                   >  2:268273/4‑150 (MQ=255)                                                                                          aTTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCg                                                        <  2:505566/151‑1 (MQ=255)                                                                                             cGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTa                                                      <  2:107545/150‑1 (MQ=255)                                                                                              ggTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTAc                                                     >  1:708504/1‑150 (MQ=255)                                                                                                 aTTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATcg                                                 <  1:392926/151‑1 (MQ=255)                                                                                                 aTTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATcg                                                 >  1:40110/1‑151 (MQ=255)                                                                                                 aTTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATcg                                                 >  1:575310/1‑151 (MQ=255)                                                                                                  ttGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATcg                                                 <  2:570888/150‑1 (MQ=255)                                                                                                   tGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATcgcg                                               >  1:440271/1‑151 (MQ=255)                                                                                                        ttCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCgcggcg                                            <  1:820989/149‑1 (MQ=255)                                                                                                        ttCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGa                                           >  1:527729/1‑150 (MQ=255)                                                                                                            ccATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATc                                      <  2:670360/151‑1 (MQ=255)                                                                                                              aTCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCa                                    <  1:343470/151‑1 (MQ=255)                                                                                                                cGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGc                                  <  2:711187/151‑1 (MQ=255)                                                                                                                cGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGc                                  >  1:424044/1‑151 (MQ=255)                                                                                                                     aaGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGCCGa                               >  1:474229/1‑149 (MQ=255)                                                                                                                       ggTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGCCGAAGAc                           >  1:453553/1‑151 (MQ=255)                                                                                                                        gTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGCCGAAGACg                          <  1:456878/151‑1 (MQ=255)                                                                                                                               gTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGCCGAAGACGCACGATa                   <  1:332977/151‑1 (MQ=255)                                                                                                                                tAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGCCGAAGACGCACGATa                   >  1:393129/1‑150 (MQ=255)                                                                                                                                 aGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGCCGAAGACGCACGATACg                 >  1:690535/1‑151 (MQ=255)                                                                                                                                     aaaGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGCCGAAGACGCACGATACGCTca             <  2:328746/151‑1 (MQ=255)                                                                                                                                     aaaGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGCCGAAGACGCACGATACGCTc              <  1:811286/150‑1 (MQ=255)                                                                                                                                     aaaGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGCCGAAGACGCACGATACGCTc              >  2:520740/1‑150 (MQ=255)                                                                                                                                      aaGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGCCGAAGACGCACGATACGCTcac            >  1:426674/1‑151 (MQ=255)                                                                                                                                      aaGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGCCGAAGACGCACGATACGCTcac            >  1:909840/1‑151 (MQ=255)                                                                                                                                                ttgGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGACCCCGTACATCGCGGCGAAGATCCAGCCGAAGACGCACGATACGCTCACTCCAATCaa   <  1:428540/150‑1 (MQ=255)                                                                                                                                                ttgGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGCCGAAGACGCACGATACGCTCACACCAATCAAc  >  1:403813/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGCCGAAGACGCACGATACGCTCACACCAATCAAC  >  NC_000913/4347177‑4347470

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