Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451 4,025,384 T→C E105E (GAA→GAG emrA ← multidrug efflux transport system ‑ membrane fusion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,025,3840TC98.2% 191.5 / NA 56E105E (GAA→GAGemrAmultidrug efflux transport system ‑ membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (26/29);  minor base A (1/0);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AGCTGCAGCCATGACATCGGGGTGCGAGGACACCACGGTGTCGTCGACCAGCGCGCTGCTGGTGCTGAGCTGTTGCCGCGCGCTGTTGACGGCAGCCTGGGCAACGCTCAGGTCGTCACGGGCGTGGGACAGTTCTTCTGCGGCAATCGCACCGCTGTCAGCCAGCACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCTTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTAGGCATCGTCGGTGCTTTCGTGCCAACGCCCGACAATGCTGTACCAGGCCCAGGTGCCGACCCCGGCAAGGATCAGCAGCAACAGCAGGCCAAGCAGCCAAGCCTTGCGCTTGCGCGACGGCTCAGGCGCGACGGAGGGAGTCG  >  AE015451/4025090‑4025676
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                      
aGCTGCAGCCATGACATCGGGGTGCGAGGACACCACGGTGTCGTCGACCAGCGCGCTGCTGGTGCTGAGCTGTTGCCGCGCGCTGTTGACGGCAGCCTGGGCAACGCTCAGGTCGTCACGGGCGTGGGACAGTTCTTCTGCGGCAATCGCACCGCTGTCAGCCAGCACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTcggcg                                                                                                                                                                                                                                                                                                  <  1:58181/301‑1 (MQ=255)
         cATGACATCGGGGTGCGAGGACACCACGGTGTCGTCGACCAGCGCGCGGCTGGTGCTGAGCTGTTGCCGCGCGCTGTTGACGGCAGCCTGGGCAACGCTCAGGTCGTCACGGGCGTGGGACAGTTCTTCTGCGGCAATCGCACCGCTGTCAGCCAGCACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAg                                                                                                                                                                                                                                                                                         <  2:42124/301‑1 (MQ=255)
           tGACATCGGGGTGCGAGGACACCACGGTGTCGTCGACCAGCGCGCTGCTGGTGCTGAGCTGTTGCCGCGCGCTGTTGACGGCAGCCTGGGCAACGCTCAGGTCGTCACGGGCGTGGGACAGTTCTTCTGCGGCAATCGCACCGCTGTCAGCCAGCACCTTGCGCCGGTTGAAGTCCTGTTTCCCCTTGCGCAGTTCGGCTTGGCGCGGTCCCGGTTGAGCCTTCAGGCAGTCGACGTTGCGGTAAAGCCCGCTCCCCTGCCGCCCGCTGCGTGCCGGATTGGCCTCGGCGCACTGCAAgac                                                                                                                                                                                                                                                                                       >  2:566054/1‑299 (MQ=255)
                          aGGACACCACGGTGTCGTCGACCAGCGCGCTGCTGGTGCTGAGCTGTTGCCGCGCGCTGTTGACGGCAGCCTGGGCAACGCTCCGGTCGTCACGGGCGTGGGACAGTTCTTCTGCGGCAATCGCACCGCTGTCAGCCAGCACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGc                                                                                                                                                                                                                                                                        >  2:238947/1‑301 (MQ=255)
                           ggcccaccgggggtcgtaGACGAGCGCGCTGCTGGTGCTGAGCTGGTGCCGCGCGCTGTTGACGGCAGCCTGGGAAACGCTCAGGTCGTCACGGGCGTGGGACAGTTCTTCTGCGGCAATCGCACCGCTGTCAGCCAGCACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGcc                                                                                                                                                                                                                                                                       <  2:606961/289‑1 (MQ=255)
                              caccacGGTGTCGTCGACCAGCGCGCTGCTGGTGCTGAGCTGTTGCCGCGCGCTGTTGACGGCAGCCTGGGCAACGCTCAGGTCGTCACGGGCGTTGGACAGTTCTTCTGCGGCAATCGCACCGCTGTCAGCCAGCACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGTCTTGCCGCGTTTCCAGTTGAGCCTTCCGCGCGTCGACGTTGCTGTATTGCCCTCGCACCTTCCGCACGTTGCGTGCCTGCTTGCCCTCGGCGGACTGTAAGGCCACTTcactatgcgccgag                                                                                                                                                                                                                                                                    >  2:457740/1‑288 (MQ=255)
                               cccacGGTGTCGTCGACCAGCGCGCTGCTGGTGCTGAGCTGTTGCCGCGCGCTGTTGACGGCAGCCTGGGCAACGCTCAGGTCGTCACGGGCGTGGGACAGTTCTTCTGCGGCAATCGCACCGCTGTCAGCCAGCACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGt                                                                                                                                                                                                                                                                   <  1:501274/300‑1 (MQ=255)
                                    ggTGTCGTCGACCAGCGCGCTGCTGGTGCTGAGCTGTTGCCGCGCGCTGTTGACGGCAGCCTGGGCAACGCTCAGGTCGTCACGGGCGTGGGACAGTTCTTCTGCGGCAATCGCACCGCTGTCAGCCAGCACCTTGCGCCGGTGGAAGTACTGTTGCGCCTTGCGCAGTTCGGCTTGGCTCGTTTCCAGTTTAGCCTTCAGCGAGTCGACGTTACTGTAAAGCCCGCTCACCTTCCGCACGCTGCGTGCCAGCTTGCCCTCTGCGGACTGCCAGCCCACTTCGCTGTCGGCCGGGTCGAAt                                                                                                                                                                                                                                                              >  2:681461/1‑301 (MQ=255)
                                     gtgtCGTCGACCAGCGCGCTGCTGGTGCTGAGCTGTTGCCGCGCGCTGTTGACGGCAGCCTGGGCAACGCTCAGGTCGTCACGGGCGTGGGACAGTTCTTCTGCGGCAATCGCACCGCTGTCAGCCAGCACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAAtt                                                                                                                                                                                                                                                             <  1:138554/301‑1 (MQ=255)
                                         ggcagACCCGCGCGCTGCTGGTGCTGAGCTGTTGCCGCGCGGTGTTGACGGCAGCCTGGGCAACGCTCAGGTCGTCACGGGCGTGGGACAGTTCTTCTGCGGCAATCGCACCGCTGTCAGCCAGCACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAACTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTgcag                                                                                                                                                                                                                                                        <  2:503745/297‑1 (MQ=255)
                                                 ggcgggcggcgggagctggtcTGTTGCCGGGCTCTGTGGACGGCGGCGTGCGGACCCCTCGGGTCGTCACGGGCGTGGGAGAGTTCTGCGACGTCAATCGCACCGTTGTCAGCCAGCACGTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGGTTGGCGGGTTTCCAGTTGAGCCATCAGCGAGTCGACGTTGCTGTAAAGGCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGACCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATagct                                                                                                                                                                                                                                                <  2:75848/281‑4 (MQ=255)
                                                   ggctgtgctgGTGCTGACCGGGTGCCACGCGCTGTATACGTCAGCCGTGTCACCGCTCCGGTCGGCACGGGCGTAGGACAGTTCTTCCGCGGCAATCGCACCGCTGTCAGCGAGCACCTTGCGCCGGTTGACGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCGTGGCCTCGGCGGGCTGTAAGGCCGCTTCGCTGTCGGGCGGGTCGAATTGCAGCAATACCTgg                                                                                                                                                                                                                                              <  2:714739/296‑1 (MQ=255)
                                                                       ggtgcggagggcggTTGACGGCAGCCTGGACAACGCTAAGGTCGTCACGGGCGTGGTACAGTTCTTCTGCGGCAATCGCACCGCTGTCATCCAGCACCTTGCGCCGGTTGACGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTGTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGACCTTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGcc                                                                                                                                                                                                                          <  2:672282/288‑1 (MQ=255)
                                                                       gTTGCCGCGCGCTGTTGACGGCAGCCTGGGCAACGCTCAGGTCGTCACGGGCGTGGGACAGTTCTTCTGCGGCAATCGCACCGCTGTCAGCCAGCACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGcc                                                                                                                                                                                                                          >  1:46532/1‑301 (MQ=255)
                                                                        ttGCCGCGCGCTGTTGACGGCAGCCTGGGCAACGCTCAGGTCGTCACGGGCGTGGGACAGTTCTTCTGCGGCAATCGCACCGCTGTCAGCCAGCACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCg                                                                                                                                                                                                                         <  1:310487/301‑1 (MQ=255)
                                                                                  ttGTTGACGGCAGCCGGGGCAACGCGCAGGTCGTCACGGGCGTGGGACAGTTCGTCTGCGGCAATCGCACCGCTGTCAGCCAGCACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGccc                                                                                                                                                                                                               <  1:164604/300‑1 (MQ=255)
                                                                                        gcGGCAGCCTGGGCAATGCTCAGGTCGTAACGGGCGTGGGACAGTTCTTCTGCGGCAATCGCACCGCTGTCAGCCAGCACCTTGCGCCGGTTGACGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGc                                                                                                                                                                                                         <  2:79955/300‑1 (MQ=255)
                                                                                         cGGCAGCCTGGGCAACGCTCAGGTCGTCACGGGCGTGGGACAGTTCTTCTGCGGCAATCGCACCGCTGTCAGCCAGCACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGGCATCGGCCCCGATGCt                                                                                                                                                                                                        >  1:724527/1‑301 (MQ=255)
                                                                                                cTGGGCAACGCTCAGGTCGTCACGGGCGTGGGACAGTTCTTCTGCGGCAATCGCACCGCTGTCAGCCAGCACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTTTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACg                                                                                                                                                                                                 >  2:112366/1‑301 (MQ=255)
                                                                                                   ggacaCGCTCAGGTCGTCACGTGCGTGGGACGGTTCTTCTGCGGCAATCGCACCGCTGTCAGCCAGCACCTTGCGGCAGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTg                                                                                                                                                                                              <  1:576742/297‑1 (MQ=255)
                                                                                                      aaCGCTCAGGTCGTCACGGGCGTGGGACAGTTCTTCTGCGGCAATCGCACCGCTGTCAGCCAGCACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCa                                                                                                                                                                                           <  1:593055/301‑1 (MQ=255)
                                                                                                       ccGCTCAGGGCGTCACGGGCGTGGGACAGGTCTTCTGCGGCAATCGCACCGCTGTCAGCCAGCACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTGGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAg                                                                                                                                                                                          <  2:72242/300‑1 (MQ=255)
                                                                                                          cTCAGGTCGTCACGGGCGTGGGACAGTTCTTCTGCGGCAATCGCACCGCTGTCAGCCCGCACCTTGCGCCGGTTGAAGCCCTGTTGCGCCTTGCGCAGTTCGGCGTGGCGCGTCTCCAGTTGAGCCTTCAGCGAGTCGACGTGGCTGTAAAGCCCGCGCACCTGCCGCACGCCGCGTGCCAGCTTGGCACCGGCGGACTGCAAGGCCACCTCCCTGTCTGCCGGGTCGCATGGCCCGCATACATGGCCGGCCTGGACCAGGTCGCCGGCTTCGGcccctctgcctgcacaggggtcacgcc                                                                                                                                                                                       >  2:500580/1‑278 (MQ=255)
                                                                                                                 cgtcACGGGCGTGGGACAGTTCTTCTGCGCCAATCGCACCGCTGTCCGCCAGCACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCTAGTCGACGTTTCTGTAAAGCCCGCGCACCTGCGCCACGCTGCGTGCGAGCTTGGCCCCGGCGGACTGCAGGGCCACTGCCCTGTCGGCCGGGTCGAATTGGAGCAATAGCTGGACGGCGTGGACCAGGGCGCCGTCATCGTCGCCGATGCGGATAACgtgaccgtcgaccggcgg                                                                                                                                                                                >  2:596877/1‑284 (MQ=255)
                                                                                                                     aCGGGCGTGGGACAGTTCTTCTGCGGCAATCGCACCGCTGTCAGCCAGCACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTg                                                                                                                                                                            <  1:414599/301‑1 (MQ=255)
                                                                                                                               gACAGTTCTTCTGCGGCAATCGCACCGCTGTCAGCCAGCACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCCGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAACGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCCTGGACCCGGCCGCCGTCATCCGCCCCGATGCCGGTAACGGTCCCAGCGACCAGCCGGGCCATTTccaccc                                                                                                                                                                  >  1:700020/1‑300 (MQ=255)
                                                                                                                                           gatctaaTCGCACCGCTGTCAGCCAGCACCTTGCGCCGGTTGAAGTCCGGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGa                                                                                                                                                      <  2:156743/296‑1 (MQ=255)
                                                                                                                                              gcgaTCGCGCCGGTGGCAGCCAGCACGTTGCGCCGGTTTAAGTCCTGTGGCGCTGTGCGCGGTTCGGCTTGGCGCGTTTCCGGTTGGGCCTTCGGCGAGGCGACGTTGCTGGAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGt                                                                                                                                                   <  2:680026/298‑1 (MQ=255)
                                                                                                                                                aaTCGCACCGCTGTCAGCCAGCACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGaagatc                                                                                                                                                 >  1:156813/1‑296 (MQ=255)
                                                                                                                                                    gCACCGCTGTCAGCCAGCACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCCTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGGCGGCCGCGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTTCCAGCGACCAGCGGGGTGCTTCCCACCACGTTGCCGTTGACGTATGCAt                                                                                                                                             >  1:608534/1‑301 (MQ=255)
                                                                                                                                                       ccGCTGTCAGCCAGCACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGAAACGGTGCCAGCGACCAGCGGGGTGATTTCACCCACGTTGCCGTTGACGTAGGCAtcct                                                                                                                                          >  1:282981/1‑299 (MQ=255)
                                                                                                                                                       ccGCGGTCCGCCAGCACCTTGCGCCGGTGGAAGTCCTGTTGGGCCTTGCGCAGTTCGGCTTGGCGCGGTTCCAGTTGAGCCTTCAGCGAGTGGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTAGGCAtcgt                                                                                                                                          <  1:77667/301‑1 (MQ=255)
                                                                                                                                                               agccagcACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTAGGCATCGTCGGTGCtt                                                                                                                                  <  1:112414/301‑1 (MQ=255)
                                                                                                                                                                gccagcACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTAGGCATCGTCGGTGCttt                                                                                                                                 <  1:681751/301‑1 (MQ=255)
                                                                                                                                                                  cagcACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTAGGCATCGTCGGTGCTTTCg                                                                                                                               >  1:573115/1‑301 (MQ=255)
                                                                                                                                                                           gggcCGGTTGAAGTCCTGTTGGGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTAGGCATCGTCGGTGCTTTCGTGCCAACGc                                                                                                                      <  1:667939/299‑1 (MQ=255)
                                                                                                                                                                                    gAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTAGGCATCGTCGGTGCTTTCGTGCCAACGCCCGACAAtt                                                                                                             >  1:486607/1‑300 (MQ=255)
                                                                                                                                                                                      aGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTCGGCATCGTCGGTGCTTTCGTGCCAACGCCCGACAATGCt                                                                                                           >  2:42299/1‑301 (MQ=255)
                                                                                                                                                                                         ccTGTTGCGCCTGGGGCAGTTCGGGTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGAACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTAGGCATCGTCGGTGCTTTCGTGCCAACGCCCGACAATGCTGTa                                                                                                        <  2:135379/301‑1 (MQ=255)
                                                                                                                                                                                                         cAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTAGGCATCGTCGGTGCTTTCGTGCCAAGGCCCGGCAATGCTGTACCAGGCCCAgggtccg                                                                                        >  1:697880/1‑296 (MQ=255)
                                                                                                                                                                                                          ggTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTAGGCATCGTCGGTGCTTTCGTGCCAACGCCCGACAATGCTGTACCAGGCCCAGGTGCCGa                                                                                       <  1:584257/300‑1 (MQ=255)
                                                                                                                                                                                                                                    ggCCTGCAGGGAGTCGGCGTTGCTGTAAAGCCGGCGCACCGGCCGAACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGAATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTAGGCATCGTCGGTGCTTTCGTGCCAACGCCCGACAATGCTGTACCAGGCCCAGGTGCCGACCCCGGCAAGGATCAGCAGCAAcagc                                                             <  2:210942/300‑1 (MQ=255)
                                                                                                                                                                                                                                       ttacagcaggtgggcGTTGCTGTAAAGTCCGCGCAGCTGCCGCACGCTGCGTGACAGCTTGGCCTCGGCGGAATGCAAGGCCACTTCGCTGTCAGCCGGGTCGAATTGCAGCAATACCTGGGCGGCATGGACCAGGTGGCCGTCATCGGCCGCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTAGGCATCGTCGGTGCTTTCGTGCCAACGCCCGACAATGCTGTACCAGGCCCAGGTGCCGACCCCGGCAAGGATCAGCAGCAACAGCAgg                                                          <  1:573857/287‑1 (MQ=255)
                                                                                                                                                                                                                                                       ttGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTAGGCATCGTCGGTGCTTTCGTGCCAACGCCCGACAATGCTGTACCAGGCCCAGGTGCCGACCCCGGCAAGGATCAGCAGCAACAGCAGGCCAAGCAGCCAAGccc                                          >  1:409241/1‑300 (MQ=255)
                                                                                                                                                                                                                                                         ttcagacgtgtgctcttccgatctgCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTAGGCATCGTCGGTGCTTTCGTGCCAACGCCCGACAATGCTGTACCAGGCCCAGGTGCCGACCCCGGCAAGGATCAGCAGCAACAGCAGGCCAAGCAGCCAAGCctt                                         <  1:507857/277‑1 (MQ=255)
                                                                                                                                                                                                                                                           tGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTAGGCATCGTCGGTGCTTTCGTGCCAACGCCCGACAATGCTGTACCAGGCCCAGGTGCCGACCCCGGCAAGGCTCAGCAGCAACAGCAGGCCAAGCAGCCAAGCctcgcc                                      >  2:357465/1‑297 (MQ=255)
                                                                                                                                                                                                                                                            gTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTAGGCATCGTCGGTGCTTTCGTGCCAACGCCCGACAATGCTGTACCAGGCCCAGGTGCCGACCCCGGCAAGGATCAGCAGCAACAGCAGGCCAAGCAGCCAAGCcttgcgc                                     >  1:97068/1‑301 (MQ=255)
                                                                                                                                                                                                                                                            gTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTAGGCATCGTCGGTGCTTTCGTGCCAACGCCCGACAATGCTGTACCAGGCCCAGGTGCCGACCCCGGCAAGGATCAGCAGCAACAGCAGGCCAAGCAGCCAAGCcttgcgc                                     >  1:635001/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                             tGCCGAACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTAGGCATCGTCGGTGCTTTCGTGCCAACGCCCGACAATGCTGTACCAGGCCCAGGTGCCGACCCCGGCAAGGATCAGCAGCAACAGCAGGCCAAGCAGCCAAGCCTTGCGCTTGCGCGACGGCTCAg                    <  1:42322/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                              gCCGCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTAGGCATCGTCGGTGCTTTCGTGCCAACGCCCGACAATGCTGTACCAGGCCCAGGTGCCGACCCCGGCAAGGATCAGCAGCAACAGCAGGCCAAGCAGCCAAGCCTTGCGCTTGCGCGACGGCTCAgg                   <  1:14333/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 gCACGCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTAGGCATCGTCGGTGCTTTCGTGCCAACGCCCGACAATGCTGTACCAGGCCCAGGTGCCGACCCCGGCAAGGATCAGCAGCAACAGCAGGCCAAGCAACCAAGccctagatcggaagagcgtcgtgtaggg                 >  2:507649/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                                     gCTGCGTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTAGGCATCGTCGGTGCTTTCGTGCCAACGCCCGACAATGCTGTACCAGGCCCAGGTGCCGACCCCGGCAAGGATCAGCAGCAACAGCAGGCCAAGCAGCCAAGCCTTGCGCTTGCGCGACGGCTCAGGCGCGACg            >  1:267375/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                          gTGCCAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTAGGCATCGTCGGTGCTTTCGTGCCAACGCCCGACAATGCTGTACCAGGCCCAGGTGCCGACCCCGGCAAGGATCAGCAGCAACAGCAGGCCAAGCAGCCAAGCCTTGCGCTTGCGCGACGGCTCAGGCGCGACggaggg       >  1:167783/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                             ccAGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTAGGCATCGTCGGTGCTTTCGTGCCAACGCCCGACAATGCTGTACCAGGCCCAGGTGCCGACCCCGGCAAGGATCAGCAGCAACAGCAGGCCAAGCAGCCAAGCCTTGCGCTTGCGCGACGGCTCAGGCGCGACGGAGGGAGt    >  1:159501/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                             ccAGCTTGGCCCCGGCTGACTGTAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGCCATCGGCCCCGATGCCGGTAACGCTGCCAGCGACCAGCGGGGTGATTTCCACCACGCTGCCGTTGCCGTCGCCATCGTCGGTGCCTTCGTGTCAACGCCCGACAATGCTGTACCCGGCCCAGGTGCCGACCCCGCCAAGGACCAGCAGCGACAGCCGACCAAGCAACAAAGCCTTGCGCTTGCCCGACGCGCCAGGCGCGGCgggggggtg    >  2:677340/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                                               aGCTTGGCCTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTAGGCATCGTCGGTGCTTTCGTGCCAACGCCCGACAATGCTGTACCAGGCCCAGGTGCCGACCCCGGCAAAGATAAGCAGCAACAGCAGACCAAGCAGCCAAGCCTTGCGCTTGCGCGCAGCCCCAGGCGCGACggagggcatcg  >  2:268267/1‑296 (MQ=255)
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                      
AGCTGCAGCCATGACATCGGGGTGCGAGGACACCACGGTGTCGTCGACCAGCGCGCTGCTGGTGCTGAGCTGTTGCCGCGCGCTGTTGACGGCAGCCTGGGCAACGCTCAGGTCGTCACGGGCGTGGGACAGTTCTTCTGCGGCAATCGCACCGCTGTCAGCCAGCACCTTGCGCCGGTTGAAGTCCTGTTGCGCCTTGCGCAGTTCGGCTTGGCGCGTTTCCAGTTGAGCCTTCAGCGAGTCGACGTTGCTGTAAAGCCCGCGCACCTGCCGCACGCTGCGTGCCAGCTTGGCTTCGGCGGACTGCAAGGCCACTTCGCTGTCGGCCGGGTCGAATTGCAGCAATACCTGGCCGGCATGGACCAGGTCGCCGTCATCGGCCCCGATGCTGGTAACGGTGCCAGCGACCAGCGGGGTGATTTCCACCACGTTGCCGTTGACGTAGGCATCGTCGGTGCTTTCGTGCCAACGCCCGACAATGCTGTACCAGGCCCAGGTGCCGACCCCGGCAAGGATCAGCAGCAACAGCAGGCCAAGCAGCCAAGCCTTGCGCTTGCGCGACGGCTCAGGCGCGACGGAGGGAGTCG  >  AE015451/4025090‑4025676

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