Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,096,388	C→A	5.7%	R354L (CGG→CTG)	rhaB ←	rhamnulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,096,388	0	C	A	5.7%	226.7 / 2.7	65	R354L (CGG→CTG)	rhaB	rhamnulokinase
Reads supporting (aligned to +/- strand):  ref base C (33/28);  new base A (2/2);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

CTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTTGCTGCTCCTGAAGCACTCG  >  NC_000913/4096242‑4096536                                                                                                                                                   |                                                                                                                                                     cTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGa                                                                                                                                                    >  1:276830/1‑149 (MQ=255)             tatGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGACAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGa                                                                                                                                        >  1:431831/1‑149 (MQ=255)              atGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGaa                                                                                                                                       <  4:146282/149‑1 (MQ=255)                gCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAAtt                                                                                                                                     >  1:72010/1‑149 (MQ=255)                    ttGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGc                                                                                                                                 <  2:201721/149‑1 (MQ=255)                    ttGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCACCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGc                                                                                                                                 <  4:387659/149‑1 (MQ=255)                       cGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGc                                                                                                                              <  4:527540/149‑1 (MQ=255)                        gagaAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGca                                                                                                                             >  1:397515/1‑149 (MQ=255)                         agaAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGcac                                                                                                                            <  1:26554/149‑1 (MQ=255)                         agaAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGcac                                                                                                                            <  1:26688/149‑1 (MQ=255)                          gaAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGcaca                                                                                                                           >  4:384800/1‑149 (MQ=255)                          gaAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGcaca                                                                                                                           >  1:328616/1‑149 (MQ=255)                           aaaTCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACAt                                                                                                                          >  2:535454/1‑149 (MQ=255)                            aaTCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATc                                                                                                                         <  1:163182/149‑1 (MQ=255)                             aTCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCg                                                                                                                        >  1:16654/1‑149 (MQ=255)                              tCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGt                                                                                                                       <  2:328616/149‑1 (MQ=255)                               cTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGtc                                                                                                                      >  1:251730/1‑149 (MQ=255)                                  cACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCa                                                                                                                   >  4:411242/1‑149 (MQ=255)                                          aGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCAGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATtaat                                                                                                           >  4:371602/1‑149 (MQ=255)                                          aGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTACGCCGTTTCCAGACACACAGCCTGAATTTCGCTGCACATCGTCTCAGGATtaat                                                                                                           >  2:124277/1‑149 (MQ=255)                                           gCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATtaata                                                                                                          >  3:391299/1‑149 (MQ=255)                                                  cAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGAt                                                                                                   <  3:384800/149‑1 (MQ=255)                                                  cAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCAGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGAt                                                                                                   <  3:371602/149‑1 (MQ=255)                                                  cAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCAGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGAt                                                                                                   <  1:124277/149‑1 (MQ=255)                                                         tGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCAtt                                                                                            >  4:505380/1‑149 (MQ=255)                                                          gCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTg                                                                                           >  3:153526/1‑149 (MQ=255)                                                          gCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTg                                                                                           <  1:513690/149‑1 (MQ=255)                                                            aaCACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTggg                                                                                         >  3:313517/1‑149 (MQ=255)                                                             acacATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGa                                                                                        <  4:299234/149‑1 (MQ=255)                                                                 aTCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGa                                                                                    >  1:5703/1‑149 (MQ=255)                                                                     gCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAAt                                                                                <  3:485018/149‑1 (MQ=255)                                                                        tACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGaa                                                                             <  2:468760/149‑1 (MQ=255)                                                                         aCAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAg                                                                            >  3:332453/1‑149 (MQ=255)                                                                          cagcagCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGc                                                                           <  3:94216/149‑1 (MQ=255)                                                                           agcagcGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCg                                                                          <  1:157308/149‑1 (MQ=255)                                                                           agcagcGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCg                                                                          >  4:453666/1‑149 (MQ=255)                                                                                  ccAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGcc                                                                   >  3:269521/1‑149 (MQ=255)                                                                                     gACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGa                                                                >  2:84197/1‑149 (MQ=255)                                                                                          tCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGc                                                           >  4:414780/1‑149 (MQ=255)                                                                                             aaaaTGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCtg                                                        <  2:276830/149‑1 (MQ=255)                                                                                              aaaTGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCtgt                                                       >  1:153966/1‑149 (MQ=255)                                                                                                aTGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTgtgt                                                     >  1:537730/1‑149 (MQ=255)                                                                                                    aGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGAtt                                                                                                              >  2:556410/1‑88 (MQ=255)                                                                                                    aGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGAtt                                                                                                              <  1:556410/88‑1 (MQ=255)                                                                                                    aGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTcgcg                                                 >  3:199319/1‑149 (MQ=255)                                                                                                     gcgcgcCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCgg                                                <  3:357037/149‑1 (MQ=255)                                                                                                         gcCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGAt                                            <  2:289101/149‑1 (MQ=255)                                                                                                             gTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAgc                                        <  2:397515/149‑1 (MQ=255)                                                                                                                  gCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGAt                                            >  3:39521/1‑140 (MQ=255)                                                                                                                  gCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGAGTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGAt                                            <  4:39521/140‑1 (MQ=255)                                                                                                                        cTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATc                             <  3:411242/149‑1 (MQ=255)                                                                                                                         tttCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATc                             >  2:552434/1‑148 (MQ=255)                                                                                                                         tttCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATc                             <  1:552434/148‑1 (MQ=255)                                                                                                                           tCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGtt                          <  2:90238/149‑1 (MQ=255)                                                                                                                           tCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGtt                          >  4:135241/1‑149 (MQ=255)                                                                                                                              gggATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGAt                       >  1:585656/1‑149 (MQ=255)                                                                                                                              gggATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGAt                       >  2:112449/1‑149 (MQ=255)                                                                                                                              gggATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGAt                       >  3:402525/1‑149 (MQ=255)                                                                                                                               ggATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGAtt                      <  1:398121/149‑1 (MQ=255)                                                                                                                                    ggTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTtgctg                 >  2:180566/1‑149 (MQ=255)                                                                                                                                         cgcCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTTGCTGCTCCt            >  3:524485/1‑149 (MQ=255)                                                                                                                                          gcCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGAt                                                                                   >  1:99388/1‑77 (MQ=255)                                                                                                                                          gcCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGAt                                                                                   <  2:99388/77‑1 (MQ=255)                                                                                                                                            cGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTTGCTGCTCCTGaa         <  2:5703/149‑1 (MQ=255)                                                                                                                                             gTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTTGCTGCTCCTGAAg        <  2:11463/149‑1 (MQ=255)                                                                                                                                                 cccGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTTGCTGCTCCTGAAGCACt    >  1:306389/1‑149 (MQ=255)                                                                                                                                                 cccGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTTGCTGCTCCTGAAGCACt    >  3:196448/1‑149 (MQ=255)                                                                                                                                                  ccGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTTGCTGCTCCTGAAGCACTc   >  1:559371/1‑149 (MQ=255)                                                                                                                                                   cGACGCGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTTGCTGCTCCTGAAGCACTCg  >  1:315391/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCGGTTGCGCCGTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTTGCTGCTCCTGAAGCACTCG  >  NC_000913/4096242‑4096536

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