Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,428,720	T→C	intergenic (‑3/‑215)	ytfB ← / → fklB	OapA family protein/FKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,428,720	0	T	C	100.0%	248.0 / NA	70	intergenic (‑3/‑215)	ytfB/fklB	OapA family protein/FKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/41);  total (29/41)

GGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTGGGATTTGTCATCATTCCCGC  >  NC_000913/4428572‑4428863                                                                                                                                                     |                                                                                                                                                ggAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGctt                                                                                                                                               >  1:265135/1‑151 (MQ=255)   aaGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTtcc                                                                                                                                             >  2:13249/1‑151 (MQ=255)    aGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTtcct                                                                                                                                            <  1:933388/151‑1 (MQ=255)      cAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTtcctcc                                                                                                                                          >  1:1185820/1‑151 (MQ=255)       aGAAATCC‑ACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCCTGAAACGAAAAATTTCCGGCGCGGGCCCGGCTTTTTCCCGGGGTGGGTTTTGGTTAAAACGCCCGGGCCAGGCTtcctcc                                                                                                                                          >  1:188056/1‑149 (MQ=255)        gAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTcctcct                                                                                                                                         <  1:1017703/150‑1 (MQ=255)        gAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTg                                                                                                                                        >  2:393382/1‑151 (MQ=255)          aaTCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCt                                                                                                                                      <  1:1123252/151‑1 (MQ=255)                accaccAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGt                                                                                                                                <  1:1085127/151‑1 (MQ=255)                     caATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCCGCGGGGCCCTGAAAGCAAAATTGTCCGGCGCGTGCCAGCCTTTTTCCAGGGGTGGTTTTAGGTCAAAGCGCCCGGGCATGGCCTCCCTCTGCTCAGCGTCTTg                                                                                                                            >  1:712111/1‑150 (MQ=255)                       aTACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTc                                                                                                                         <  2:899537/151‑1 (MQ=255)                       aTACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTc                                                                                                                         <  1:300447/151‑1 (MQ=255)                         aCAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTcc                                                                                                                        <  1:1005423/150‑1 (MQ=255)                         aCAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCt                                                                                                                       <  2:631780/151‑1 (MQ=255)                          cAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTc                                                                                                                      >  2:890712/1‑151 (MQ=255)                               gcggcAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATc                                                                                                                  >  2:1150157/1‑150 (MQ=255)                               gcggcAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCa                                                                                                                 <  1:13249/151‑1 (MQ=255)                                   cAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGc                                                                                                             <  2:360793/151‑1 (MQ=255)                                       aatGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTa                                                                                                          <  2:1121392/150‑1 (MQ=255)                                       aatGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTAt                                                                                                         <  1:1018292/151‑1 (MQ=255)                                             cccccGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATa                                                                                                    >  1:123720/1‑150 (MQ=255)                                             cccccGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATa                                                                                                    >  1:775607/1‑150 (MQ=255)                                                ccGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCt                                                                                                <  2:227432/151‑1 (MQ=255)                                                        gCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTat                                                                                        <  2:185472/151‑1 (MQ=255)                                                        gCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTat                                                                                        >  2:1045093/1‑151 (MQ=255)                                                         cATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTatg                                                                                       <  2:570692/151‑1 (MQ=255)                                                          aTCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAATTCAAAGCGCCCGGGCATGGCCTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTat                                                                                        >  1:759653/1‑149 (MQ=255)                                                            cggcggCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTatgatg                                                                                    <  1:176978/151‑1 (MQ=255)                                                             ggcggcAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGt                                                                                   >  2:445055/1‑151 (MQ=255)                                                                 gcAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGt                                                                               <  1:887232/151‑1 (MQ=255)                                                                     cGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTatgatg                                                                                    <  2:156122/142‑1 (MQ=255)                                                                     cGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTatgatg                                                                                    >  1:156120/1‑142 (MQ=255)                                                                     cGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGt                                                                            >  1:135010/1‑150 (MQ=255)                                                                        gTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAAc                                                                        >  1:234154/1‑151 (MQ=255)                                                                         tCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAAcc                                                                       >  1:151634/1‑151 (MQ=255)                                                                          ccATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAAccc                                                                      <  1:300503/151‑1 (MQ=255)                                                                           cATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCg                                                                     <  2:403914/151‑1 (MQ=255)                                                                           cATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCg                                                                     <  2:188058/151‑1 (MQ=255)                                                                              gAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCt                                                                  <  2:293685/151‑1 (MQ=255)                                                                                    aaaaTTGTCCGGCGCGTGCCAGATTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTCTGATGTCCGTGGTAAACCCGCCTttatt                                                             <  2:322118/150‑1 (MQ=255)                                                                                    aaaaTTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTttat                                                              <  1:937293/149‑1 (MQ=255)                                                                                     aaaTTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAGCCCGCCTTtattat                                                           <  1:855195/151‑1 (MQ=255)                                                                                       aTTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTatat                                                         <  1:561786/151‑1 (MQ=255)                                                                                         tGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATAtt                                                        <  2:866488/150‑1 (MQ=255)                                                                                          gTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCGGGGCATGGCTTCCTCCTGCTCCGCGTTTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTatat                                                         <  1:911262/148‑1 (MQ=255)                                                                                          gTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATAtt                                                        <  2:522513/149‑1 (MQ=255)                                                                                          gTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTc                                                       >  1:826104/1‑150 (MQ=255)                                                                                                       cAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGc                                         <  2:196723/151‑1 (MQ=255)                                                                                                            tttCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTg                                    >  2:437887/1‑151 (MQ=255)                                                                                                             ttCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGt                                   >  2:186609/1‑151 (MQ=255)                                                                                                                 ccAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGTTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGtt                                  <  1:692832/148‑1 (MQ=255)                                                                                                                   aGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCttt                             >  2:389063/1‑151 (MQ=255)                                                                                                                      gttggttTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTc                          >  1:832726/1‑151 (MQ=255)                                                                                                                       ttgtttTTAGTTCAAAGCGCCCGGGAATGGCTCCCTCTTGCTCCGCGTTTTGTTCCTAAATCATAGCCTATGAATAAGCTAACGTTTTGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCttttt                           <  2:770281/149‑1 (MQ=255)                                                                                                                        tggttTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTTTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTcc                         <  2:1061687/150‑1 (MQ=255)                                                                                                                         ggttTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTccc                        <  2:267014/150‑1 (MQ=255)                                                                                                                          gttTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTg                      >  1:24554/1‑151 (MQ=255)                                                                                                                           ttttAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTgg                     >  1:582983/1‑151 (MQ=255)                                                                                                                            tttAGTTCAAGGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTgg                     <  1:770195/150‑1 (MQ=255)                                                                                                                            tttAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTggg                    >  1:220031/1‑151 (MQ=255)                                                                                                                              tAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTGGGAt                  >  2:384687/1‑151 (MQ=255)                                                                                                                                 ttCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTGTTGTTCCTCAATCATAGCCTATGAATAGGCTAACCCTATGATGTCCTTGTTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTGGGAttt                <  1:1145386/150‑1 (MQ=255)                                                                                                                                  tCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTGGGATTTGt              >  1:831055/1‑151 (MQ=255)                                                                                                                                    aaaGCTCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTGGGATTTGt              <  1:139875/149‑1 (MQ=255)                                                                                                                                     aaGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTGGGATTTGtcat           >  2:658161/1‑151 (MQ=255)                                                                                                                                         gcCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTGGGATTTGTCATCAtt       >  2:830555/1‑151 (MQ=255)                                                                                                                                          cccGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTGGGATTTGTCATCAtt       <  2:1106171/150‑1 (MQ=255)                                                                                                                                           ccGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTGGGATTTGTCATCAtt       <  2:1175617/149‑1 (MQ=255)                                                                                                                                           ccGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTGGGATTTGTCATCATTcc     >  1:694431/1‑151 (MQ=255)                                                                                                                                             gggCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTGGGATTTGTCATCATTCCcg   <  2:977596/151‑1 (MQ=255)                                                                                                                                              ggCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCTTATGAATAAGCTAACGCTATGATGTCCGTGTTAACTCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTGGGATTTGTCATCATTCCcgc  <  1:625032/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                GGAAGCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTGGGATTTGTCATCATTCCCGC  >  NC_000913/4428572‑4428863

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