Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,096,230	T→C	T407A (ACG→GCG)	rhaB ←	rhamnulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,096,230	0	T	C	98.7%	269.4 / ‑3.8	77	T407A (ACG→GCG)	rhaB	rhamnulokinase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (32/44);  total (33/44)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

GTGCTGACGACCTGACGGAAATCATCCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGTGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCG  >  NC_000913/4096082‑4096373                                                                                                                                                     |                                                                                                                                                gtgCTGACGACCTGACGGAAATCATCCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCgcgt                                                                                                                                               <  2:837628/151‑1 (MQ=255)  tgCTGACGACCTGACGGAAATCATCCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGtt                                                                                                                                              >  1:662232/1‑151 (MQ=255)   gCTGACGACCTGACGGAAATCATCCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTc                                                                                                                                             <  2:107407/151‑1 (MQ=255)      gacgacCTGACGGAAATCATCCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTgg                                                                                                                                          <  2:450971/151‑1 (MQ=255)      gacgacCTGACGGAAATCATCCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTgg                                                                                                                                          >  2:822348/1‑151 (MQ=255)      gacgacCTGACGGAAATCATCCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTgg                                                                                                                                          <  1:121719/151‑1 (MQ=255)       acgacCTGACGGAAATCATCCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGc                                                                                                                                         <  1:8343/151‑1 (MQ=255)           ccTGACGGAAATCATCCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGcc                                                                                                                                     >  1:288212/1‑151 (MQ=255)              gACGGAAATCATCCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGGTATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTcc                                                                                                                                  <  2:384608/151‑1 (MQ=255)              gACGGAAATCATCCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGAAGCCt                                                                                                                                    <  1:1063579/149‑1 (MQ=255)               aCGGAAATCATCCACATTGTTGAGTTCATCCAGCGGCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTAGCAGCCTcc                                                                                                                                  >  1:1081979/1‑150 (MQ=255)                cGGAAATCATCCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTccgc                                                                                                                                <  2:501803/151‑1 (MQ=255)                  gAAATCATCCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTccgccg                                                                                                                              <  1:177845/151‑1 (MQ=255)                      tcatcCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACa                                                                                                                           >  1:851673/1‑150 (MQ=255)                      tcatcCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACa                                                                                                                           >  2:392617/1‑150 (MQ=255)                       catcCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACaa                                                                                                                          <  1:938750/150‑1 (MQ=255)                       catcCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAat                                                                                                                         <  1:620550/151‑1 (MQ=255)                       catcCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAat                                                                                                                         <  2:239300/151‑1 (MQ=255)                              aTTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGt                                                                                                                  >  1:916368/1‑151 (MQ=255)                               ttgttgAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGtt                                                                                                                 >  1:24730/1‑151 (MQ=255)                                 gttgAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTg                                                                                                                <  1:548919/150‑1 (MQ=255)                                   tgAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCg                                                                                                              <  1:248231/150‑1 (MQ=255)                                        tCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGaaa                                                                                                         <  1:1272875/150‑1 (MQ=255)                                          aTCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCt                                                                                                      >  2:128624/1‑151 (MQ=255)                                          aTCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCt                                                                                                      <  1:1044558/151‑1 (MQ=255)                                           tCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCtt                                                                                                     >  1:578363/1‑151 (MQ=255)                                            ccAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTc                                                                                                    >  1:649331/1‑151 (MQ=255)                                             cAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTc                                                                                                    >  1:1013481/1‑150 (MQ=255)                                              tgCGTCATTAACTGGATGCCGATATTGACGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGATGAGCAGCGCGTTCTGGCAGCCTCCGCCGAAAATATGCAGTTGCGAGAAATCTTCAc                                                                                                  >  1:320841/2‑151 (MQ=255)                                                 gTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACcgc                                                                                               >  2:1253490/1‑151 (MQ=255)                                                   cATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCgcgc                                                                                             <  2:1086932/151‑1 (MQ=255)                                                   cATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCgcgc                                                                                             >  2:375996/1‑151 (MQ=255)                                                    aTTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCgcgc                                                                                             >  1:695179/1‑150 (MQ=255)                                                       aaCTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCt                                                                                         >  2:685610/1‑151 (MQ=255)                                                        aCTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTg                                                                                        <  2:18410/151‑1 (MQ=255)                                                        aCTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTg                                                                                        <  1:836859/151‑1 (MQ=255)                                                           ggATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTgcg                                                                                      <  2:1269408/150‑1 (MQ=255)                                                             aTGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCa                                                                                   <  1:991336/151‑1 (MQ=255)                                                                 cGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTc                                                                               <  1:1237640/151‑1 (MQ=255)                                                                 cGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTc                                                                               >  2:734388/1‑151 (MQ=255)                                                                     atTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGc                                                                           <  2:828031/151‑1 (MQ=255)                                                                        gCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAac                                                                        >  2:280341/1‑151 (MQ=255)                                                                          cGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCGAGCTCATGCAacac                                                                      <  1:937306/151‑1 (MQ=255)                                                                          cGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAacac                                                                      >  2:858533/1‑151 (MQ=255)                                                                               gTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCgg                                                                 <  2:700724/151‑1 (MQ=255)                                                                                 cGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCa                                                               <  2:1106603/151‑1 (MQ=255)                                                                                  gAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCAt                                                              <  2:315291/151‑1 (MQ=255)                                                                                   aGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCAt                                                              >  2:1145235/1‑150 (MQ=255)                                                                                       ttCAACAGGCCCTGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGCCAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATAcagc                                                         >  1:218183/1‑151 (MQ=255)                                                                                        tCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATAcag                                                          <  2:1152388/149‑1 (MQ=255)                                                                                         cAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATAcagca                                                        <  1:271331/150‑1 (MQ=255)                                                                                          accAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATAcagca                                                        <  1:550843/147‑1 (MQ=255)                                                                                          aaCAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACagcagc                                                      <  1:1208714/151‑1 (MQ=255)                                                                                           aCAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAgcg                                                     >  2:31436/1‑151 (MQ=255)                                                                                                cccGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGa                                                <  2:422707/151‑1 (MQ=255)                                                                                                 ccGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGAc                                               <  2:450450/151‑1 (MQ=255)                                                                                                 ccGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGAc                                               >  1:483866/1‑151 (MQ=255)                                                                                                     cGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTc                                           >  1:127766/1‑151 (MQ=255)                                                                                                           cccGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGaaaa                                      >  2:171256/1‑150 (MQ=255)                                                                                                           cccGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAAt                                     >  2:605325/1‑151 (MQ=255)                                                                                                             cGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGc                                   <  1:419576/151‑1 (MQ=255)                                                                                                              gAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCa                                  <  2:935070/151‑1 (MQ=255)                                                                                                              gAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTGCTGGCAGCCTCCGCCGACAATATGCAGTTGCGATATATCTTCACCGCGCAGCTGCGCCAGCTCATACAACACACCGCCATACAGCAGCGCCAGACTGTCGAAAATGCa                                  <  1:970244/151‑1 (MQ=255)                                                                                                               aaTACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGTGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGc                                   >  1:1077796/1‑149 (MQ=255)                                                                                                                  aCCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAgcgc                              <  2:1232838/151‑1 (MQ=255)                                                                                                                   ccGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAgcgcg                             >  2:673951/1‑151 (MQ=255)                                                                                                                         gCATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGtt                       >  1:1036413/1‑151 (MQ=255)                                                                                                                          cATCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTc                      <  1:546687/151‑1 (MQ=255)                                                                                                                           aTCGGCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCa                     <  1:509099/151‑1 (MQ=255)                                                                                                                              ggCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCa                  <  1:304988/151‑1 (MQ=255)                                                                                                                              ggCGCATAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGGCGAAAATGCAGCGCGCCAGTTCAGc                   >  1:813988/1‑150 (MQ=255)                                                                                                                                   aTAGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACt             >  2:1071081/1‑151 (MQ=255)                                                                                                                                     aGCTGGTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCGCCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCTGACTGTCGAAAATGCTGCGCGCCAGTTCAGCATCTCttt           <  1:916607/151‑1 (MQ=255)                                                                                                                                         ggTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCgg       <  1:596247/151‑1 (MQ=255)                                                                                                                                          gTTGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGTCATACATCTTCTCCTGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTcc         <  1:1044713/148‑1 (MQ=255)                                                                                                                                            tGAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGa     <  2:518237/150‑1 (MQ=255)                                                                                                                                             gAGCAGCGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATc   <  2:288227/151‑1 (MQ=255)                                                                                                                                              agcagcGCGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCg  >  1:85362/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                GTGCTGACGACCTGACGGAAATCATCCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAATACCGCAGGCATCGGCGCATAGCTGGTTGAGCAGCGTGTTCTGGCAGCCTCCGCCGACAATATGCAGTTGCGAGAAATCTTCACCGCGCAGCTGCGCCAGCTCATGCAACACATCGGCATACAGCAGCGCCAGACTGTCGAAAATGCAGCGCGCCAGTTCAGCATCACTTTCCGGGATCG  >  NC_000913/4096082‑4096373

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