Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,964,454 C→T 23.3% D302N (GAC→AAC)  fleQ ← Transcriptional regulator FleQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,964,4540CT23.3% 77.9 / 25.7 43D302N (GAC→AAC) fleQTranscriptional regulator FleQ
Reads supporting (aligned to +/- strand):  ref base C (15/18);  new base T (5/5);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.41e-01

ACGGGTGCATGATCGCCATGCGCTCGACCAGGTTGGCCAGCTCGCGTACGTTGCCCGGCCAGCCGTGACGACACAGCGACATGATCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGAT  >  AE015451/4964226‑4964679
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                 
aCGGGTGCATGATCGCCATGCGCTCGACCAGGTTGGCCAGCTCGCGTACGTTGCCCGGCCAGCCGTGACGACACAGCGACATGATCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTTTTcgcg                                                                                                                                                                                                                             >  1:560246/1‑235 (MQ=255)
   ggTGCATGATCGCCATGCGCTCGACCAGGTTGGCCAGCTCGCGTACGTTGCCCGGCCAGCCGTGACGACACAGCGACATGATCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAg                                                                                                                                                                                                                         <  1:664845/236‑1 (MQ=255)
   ggTGCATGATCGCCATGCGCTCGACCAGGTTGGCCAGCTCGCGTACGTTGCCCGGCCAGCCGTGACGACACAGCGACATGATCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAg                                                                                                                                                                                                                         <  1:63738/236‑1 (MQ=255)
        aTGATCGCCATGCGCTCGACCAGGTTGGCCAGCTCGCGTACGTTGCCCGGCCAGCCGTGACGACACAGCGACATGATCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCt                                                                                                                                                                                                                                  >  2:646536/1‑222 (MQ=255)
        aTGATCGCCATGCGCTCGACCAGGTTGGCCAGCTCGCGTACGTTGCCCGGCCAGCCGTGACGACACAGCGACATGATCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCt                                                                                                                                                                                                                                  <  1:646537/222‑1 (MQ=255)
        aTGATCGCCATGCGCTCGACCAGGTTGGCCAGCTCGCGTACGTTGCCCGGCCAGCCGTGACGACACAGCGACATGATCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGcc                                                                                                                                                                                                                    <  1:190819/236‑1 (MQ=255)
                                           gcgTACGTTGCCCGGCCAGCCGTGACGACACAGCGACATGATCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTTTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTtgtg                                                                                                                                                                                      <  2:295434/231‑1 (MQ=255)
                                                     cccGGCCAGCCGTGACGACACAGCGACATGATCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTGCGAGGTTCTTGTGCGTCg                                                                                                                                                                                 >  2:501232/1‑226 (MQ=255)
                                                           cAGCCGTGACGACACAGCGACATGATCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTTTTCGCGGAAg                                                                                                                                                                                                                         <  2:550519/180‑1 (MQ=255)
                                                           cAGCCGTGACGACACAGCGACATGATCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTTTTCGCGGAAg                                                                                                                                                                                                                         >  1:550520/1‑180 (MQ=255)
                                                            aGCCGTGACGACACAGCGACATGATCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACg                                                                                                                                                                <  1:141146/236‑1 (MQ=255)
                                                                        acaGCGACATGATCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTgg                                                                                                                                                    <  1:172881/236‑1 (MQ=255)
                                                                                   aTCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCtt                                                                                                                                                                                         >  2:426569/1‑188 (MQ=255)
                                                                                      gAGGCGGAATTGAAGCGAATGTAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTTTTCGCGGAAGGTGcc                                                                                                                                                                                                                    <  2:304809/158‑1 (MQ=255)
                                                                                                          ggAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTggg                                                                                                                                                   >  1:460571/1‑203 (MQ=255)
                                                                                                            acccGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCt                                                                                                                                                <  2:227391/202‑1 (MQ=255)
                                                                                                                    ttCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGc                                                                                                                                          >  1:743513/1‑202 (MQ=255)
                                                                                                                      cgcGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCacac                                                                                                      <  1:108242/234‑1 (MQ=255)
                                                                                                                         gTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGACGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGAGCACGTCGATGCTTTGGGTCTTGTTGCTACCCaa                                                                                                                                   >  1:465959/1‑203 (MQ=255)
                                                                                                                                 tGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTTTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCagcagc                                                                                              >  2:252368/1‑233 (MQ=255)
                                                                                                                                    gcgAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGAc                                                                                         >  1:79227/1‑235 (MQ=255)
                                                                                                                                      gAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCacac                                                                                                      <  1:233561/220‑1 (MQ=255)
                                                                                                                                               tCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACGCGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCt                                                                                                                                                <  2:208570/169‑1 (MQ=255)
                                                                                                                                                 gTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCTTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTAc                                                                                                                                       >  2:438804/1‑176 (MQ=255)
                                                                                                                                                     agcaGCAACGGTATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAAc                                                                            <  1:181219/229‑1 (MQ=255)
                                                                                                                                                                  tGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTtcg                                                          <  1:294684/236‑1 (MQ=255)
                                                                                                                                                                               cTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACg                                                                                                                    <  1:28147/165‑1 (MQ=255)
                                                                                                                                                                                tCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTTTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAggg                                             <  1:252369/235‑1 (MQ=255)
                                                                                                                                                                                 cACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTTTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAgg                                              >  2:648895/1‑233 (MQ=255)
                                                                                                                                                                                 cACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCGCACGCAGCAGCTTGACCTGCATCGGCAAc                                                                            >  1:157639/1‑203 (MQ=255)
                                                                                                                                                                                 cACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTa                                           >  1:532678/1‑236 (MQ=255)
                                                                                                                                                                                 cACGCAGCGCCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCg                                                                                                                                                                                                 <  2:581181/86‑1 (MQ=255)
                                                                                                                                                                                             ccATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGc                                >  1:166102/1‑235 (MQ=255)
                                                                                                                                                                                                 ttCGATAGGGAACACGTTCAGGCGGTAGTACAGGATTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCgatga                                                                                                                                                                        >  2:590641/1‑95 (MQ=255)
                                                                                                                                                                                                     aTAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGtt                                   <  1:561956/224‑1 (MQ=255)
                                                                                                                                                                                                              acGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCa                                                                                                                                                                               <  2:730317/75‑1 (MQ=255)
                                                                                                                                                                                                              acGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCa                                                                                                                                                                               >  1:730318/1‑75 (MQ=255)
                                                                                                                                                                                                                gTTCAGGCGTTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGGATCGGCAACGGCATGTCGCCGATTTCGTCGAg                                                     <  2:718507/195‑1 (MQ=255)
                                                                                                                                                                                                                    aGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAgg                                              <  2:631237/198‑1 (MQ=255)
                                                                                                                                                                                                                      gCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCt        >  1:58208/1‑234 (MQ=255)
                                                                                                                                                                                                                      gCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCACAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCt        >  2:404884/1‑234 (MQ=255)
                                                                                                                                                                                                                        ggTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCgg               >  1:653685/1‑225 (MQ=255)
                                                                                                                                                                                                                          tagtaCAGGTTTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGAt  <  1:648896/236‑1 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                 
ACGGGTGCATGATCGCCATGCGCTCGACCAGGTTGGCCAGCTCGCGTACGTTGCCCGGCCAGCCGTGACGACACAGCGACATGATCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGAT  >  AE015451/4964226‑4964679

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