Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 3,029,299 A→G 6.9% T318A (ACC→GCC) ‡ pcaY → aromatic acid chemoreceptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154513,029,2990AG 6.9% 187.3 / 3.8 58T318A (ACC→GCC) ‡pcaYaromatic acid chemoreceptor
Reads supporting (aligned to +/- strand):  ref base A (24/30);  new base G (2/2);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.28e-01

GCGACCTGAGTGTGCCGGTGCAAGGGCAGGGTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGACCAGCAGTGCGCACCTGCAGCAGTTGGTGGGGCAGGTGCGGGACATCGGCAAGGTGCTGGAAGTGATCCGCTCGGTGTCCGAGC  >  AE015451/3029066‑3029524
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                 
gCGACCGGAGTGTGCCGGTGCAAGGGCAGGGTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGAc                                                                                                                                                                                                                                  <  1:932009/235‑1 (MQ=255)
     cTGAGTGTGCCGTTGCAAGGGCAGGGTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGAcc                                                                                                                                                                                                                                 <  1:200866/231‑1 (MQ=255)
          tgtgCCGGTGCAAGGGCAGGGTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTgg                                                                                                                                                                                                                       <  1:781479/236‑1 (MQ=255)
                       gggCAGGGTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGGGGAGGAAGTTgcgcg                                                                                                                                                                                                          >  1:173084/1‑236 (MQ=255)
                        ggCAGGGTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTgcgcgc                                                                                                                                                                                                         <  1:415112/236‑1 (MQ=255)
                              gTAACAACGGAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAAcgc                                                                                                                                                                                                    <  1:846333/235‑1 (MQ=255)
                               tAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAg                                                                                                                                                                                                                 >  2:630566/1‑221 (MQ=255)
                               tAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCa                                                                                                                                                                                                        >  1:610079/1‑230 (MQ=255)
                                      gAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATcaccac                                                                                                                                                                                           >  1:309325/1‑236 (MQ=255)
                                                               tCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCa                                                                                                                                                                    >  2:916376/1‑234 (MQ=255)
                                                                  aGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGGTCATTTCGGTCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGGGAATCCAAcc                                                                                                                                                                 >  2:499308/1‑234 (MQ=255)
                                                                  aGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAAc                                                                                                                                                                  >  2:595894/1‑233 (MQ=255)
                                                                    cgcATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCt                                                                                                                                                                                         >  2:366382/1‑208 (MQ=255)
                                                                           agcagAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAaccac                                                                                                                                                                               >  1:367386/1‑211 (MQ=255)
                                                                           agcagAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAaccac                                                                                                                                                                               <  2:367383/211‑1 (MQ=255)
                                                                                        cgACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAgccgcc                                                                                                                                            <  2:543407/233‑1 (MQ=255)
                                                                                          acacTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCa                                                                                                                                        >  1:476860/1‑235 (MQ=255)
                                                                                                   ggCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGa                                                                                                                               >  1:362997/1‑235 (MQ=255)
                                                                                                           aacaacTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTgcgc                                                                                                                                                                                                           >  2:547843/1‑151 (MQ=255)
                                                                                                           aacaacTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCACGTCAGCGAAAACATCGa                                                                                                                      >  1:54102/1‑236 (MQ=255)
                                                                                                           aacaacTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGa                                                                                                                      >  1:287632/1‑236 (MQ=255)
                                                                                                             caacTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGa                                                                                                                      <  2:389415/234‑1 (MQ=255)
                                                                                                                cTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGcc                                                                                                                 >  1:539865/1‑235 (MQ=255)
                                                                                                                     cccGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCAcc                                                                                                               >  2:494366/1‑233 (MQ=255)
                                                                                                                      ccGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAg                                                                                                           >  1:264743/1‑236 (MQ=255)
                                                                                                                           ggCTGGGCGCCGCCGCCACGGCGTTGGACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATca                                                                                                                                                                                               <  2:174238/146‑1 (MQ=255)
                                                                                                                                    ccgccgccACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTc                                                                                                                                                                                        <  2:277830/145‑1 (MQ=255)
                                                                                                                                       ccgccACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGGTCATTTCGGTCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAacctcacc                                                                                                                                                                            <  2:431358/154‑6 (MQ=255)
                                                                                                                                       ccgccACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGGTCATTTCGGTCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAacctcacc                                                                                                                                                                            >  1:431361/1‑149 (MQ=255)
                                                                                                                                              ggCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGAcc                                                                                    <  1:630569/235‑1 (MQ=255)
                                                                                                                                                    gAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGACCAGCAGTg                                                                             <  1:502841/236‑1 (MQ=255)
                                                                                                                                                     aaCGCCGTCACCGAGGAAAGCGCCAACAACCTGCGGCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCaccacc                                                                                                                                                                                          <  2:161675/126‑1 (MQ=255)
                                                                                                                                                      aCGCCGTCCCCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGCGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGaa                                                                                             <  1:171906/218‑1 (MQ=255)
                                                                                                                                                      aCGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGACCAGCAGTgcg                                                                           <  1:839383/236‑1 (MQ=255)
                                                                                                                                                             cACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGACCAGCAGTGCGCAcc                                                                       <  2:362994/233‑1 (MQ=255)
                                                                                                                                                               ccGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGACCAGCAGTGCGCACCTgcag                                                                  >  1:322015/1‑236 (MQ=255)
                                                                                                                                                                     aaaGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGACCAGCAGTGCGCACCTGCAGCAGt                                                              <  2:476857/234‑1 (MQ=255)
                                                                                                                                                                                 acCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGGTCATTTCGGTCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAATATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGACCAGCAGTGCGCACCTGCAGCAGTTGGTGGGGCAGGt                                                 <  1:499311/235‑1 (MQ=255)
                                                                                                                                                                                         ccgcagGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGACCAGCAGTg                                                                             <  2:399695/197‑1 (MQ=255)
                                                                                                                                                                                            cagGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGACCAGCAGTGCGCACCTGCAGCAGTTGGTGGGGCAGGTGCGGGACATCgg                                     <  1:628789/236‑1 (MQ=255)
                                                                                                                                                                                            cagGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGACCAGCAGTGCGCACCTGCAGCAGTTGGTGGGGCAGGTGCGGGACATCgg                                     <  1:595897/236‑1 (MQ=255)
                                                                                                                                                                                            cagGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGACCAGCAGTGCGCACCTGCAGCAGTTGGTGGGGCAGGTGCGGGACATCg                                      <  1:75992/235‑1 (MQ=255)
                                                                                                                                                                                               gggCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGACCAGCAGTGCGCACCTGCAGCAGTTGGTGGGGCAGGTGCGGGACATCGGCaa                                  >  1:496442/1‑236 (MQ=255)
                                                                                                                                                                                                  cAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGACCAGCAGTGCGCACCTGCAGCAGTTGGTGGGGCAgg                                                  >  2:397155/1‑217 (MQ=255)
                                                                                                                                                                                                         tGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGACCAGCAGTGCGCACCTGCAGCAGTtggtgg                                                        >  2:417151/1‑204 (MQ=255)
                                                                                                                                                                                                         tGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGACCAGCAGTGCGCACCTGCAGCAGTtggtgg                                                        <  1:417154/204‑1 (MQ=255)
                                                                                                                                                                                                             gCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGACCAGCAGTGCGCACCTGCAGCAGTTGGTGGGGCAGGTGCGGGACATCGGCAAGGTGCTGGAAGTg                     <  1:481450/235‑1 (MQ=255)
                                                                                                                                                                                                              cAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAg                                                                                                                                  >  2:897871/1‑125 (MQ=255)
                                                                                                                                                                                                                ggCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGACCAGCAGTGCGCACCTGCAGCAGTTGGTGGGGCAGGTGCGGGACATCGGCAAGGTGCTGGAAGTGATcc                 <  1:497943/236‑1 (MQ=255)
                                                                                                                                                                                                                 gccgccACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGACCAGCAGTGCGCACCTGCAGCAGTTGGTGGGGCAGGTGCGGGACATCGGCAAGGTGCTGGAAGTGATCCg                <  1:293004/236‑1 (MQ=255)
                                                                                                                                                                                                                   cgccACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGACCAGCAGTGCGCACCTGCAGCAGTTGGGGGGGCAGGTGCGGGACATCGGCAAGGTGCTGGAAGTGATCCGCt              >  1:35536/1‑236 (MQ=255)
                                                                                                                                                                                                                      cACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGAc                                                                                     <  2:221124/162‑1 (MQ=255)
                                                                                                                                                                                                                       aCTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACAAGTGAATCAAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATg                                                                                                      >  2:811846/1‑144 (MQ=255)
                                                                                                                                                                                                                       aCTGCCGTGACCGAAATGACCACGGCAGAGGAGGAAGTTGCGCGCAACGCGAGCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGACCAGCAGTGCGCACCTGCAGCAGTTGGTGGGGCAGGTGCGGGACATCgg                                     <  2:38422/209‑1 (MQ=255)
                                                                                                                                                                                                                                 ccGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGACCAGCAGTGCGCACCTGCAGCAGTTGGTGGGGCAGGTGCGGGACATCGGCAAGGTGCTGGAAGTGATCCGCTCGGTGTCCGAGc  <  2:496439/234‑1 (MQ=255)
                                                                                                                                                                                                                                   gAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGAc                                                                                     <  2:528679/149‑1 (MQ=255)
                                                                                                                                                                                                                                   gAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGAc                                                                                     >  1:528682/1‑149 (MQ=255)
                                                                                                                                                                                                                                        ggACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGACCAGCAGTGCGCACCTGCAGCAGTTGGTGGGGCAGGTGCGGGACATCgg                                     <  2:173082/192‑1 (MQ=255)
                                                                                                                                                                                                                                          aCCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCGGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGACCAGCAGTGCGCACCTGCAGCAGTTGGTGGGGCAGGTGCGGGACATCg                                      <  2:127425/190‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                 
GCGACCTGAGTGTGCCGGTGCAAGGGCAGGGTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGCACCGAAGCCATGACCCGCGAAATCCAGACCAGCAGTGCGCACCTGCAGCAGTTGGTGGGGCAGGTGCGGGACATCGGCAAGGTGCTGGAAGTGATCCGCTCGGTGTCCGAGC  >  AE015451/3029066‑3029524

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