Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 1,397,940 T→G 9.6% P296P (CCT→CCG) ‡ tolA → colicin S4 and filamentous phage transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,397,9400TG 9.6% 164.0 / 7.5 52P296P (CCT→CCG) ‡tolAcolicin S4 and filamentous phage transport system
Reads supporting (aligned to +/- strand):  ref base T (15/32);  new base G (3/2);  total (18/34)
Fisher's exact test for biased strand distribution p-value = 3.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.28e-01

GCAGCTGAAGAGGCCAAGAAAAAGGCCGCTGCAGATGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTGGCTGCGGTGAAGAACATTGGTCGTTTGACCGAGATGCAGGGTATGAAGCCGAGCGATTTCAACCAATATCGTTCGTTCAAGATGACATTTACACCTGAGGATCTAGCGT  >  AE015451/1397713‑1398166
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                  
gCAGCTGAAGAGGCCAAGAAAAAGGCCGCTGCAGATGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCcgcgcg                                                                                                                                                                                                                            <  1:451834/236‑1 (MQ=255)
      gAAGAGGCCAAGAAAAAGGCCGCTGCAGATGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAg                                                                                                                                                                                                                        <  1:725464/234‑1 (MQ=255)
       aaGAGGCCAAGAAAAAGGCCGCTGCAGATGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAggg                                                                                                                                                                                                                      <  1:673754/235‑1 (MQ=255)
                               cAGATGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCAt                                                                                                                                                                                                                   >  2:1089195/1‑214 (MQ=255)
                                   tGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCcgcg                                                                                                                                                                                                                              <  2:749072/199‑1 (MQ=255)
                                   tGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATg                                                                                                                                                                                          <  1:848239/235‑1 (MQ=255)
                                   ggCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACggtg                                                                                                                                                                                                            <  2:486993/216‑1 (MQ=255)
                                    gCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCatgcg                                                                                                                                                                                                                      <  1:152419/206‑5 (MQ=255)
                                      ccAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGAAGGCGGACGTGCCCAGCCCTCGGCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTg                                                                                                                                                                                       <  1:407518/235‑1 (MQ=255)
                                       cAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGAAGGCGGACGTGCCCAGCCCTCGGCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGt                                                                                                                                                                                         >  2:407517/1‑232 (MQ=255)
                                         gaagaagaaGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCg                                                                                                                                                                                    <  1:849978/235‑1 (MQ=255)
                                             aagaagGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCcgcgc                                                                                                                                                                                                                             <  2:446118/190‑1 (MQ=255)
                                                aagGCCCAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTcc                                                                                                                                                                                                                                 <  2:106850/183‑1 (MQ=255)
                                                     acaGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGTTCCTGCAGATCAACATGTTGCCGGa                                                                                                                                                                                  >  1:378121/1‑225 (MQ=255)
                                                       aGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGc                                                                                                                                                                                                     >  2:1101435/1‑204 (MQ=255)
                                                       aGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCAc                                                                                                                                                                      >  1:11555/1‑235 (MQ=255)
                                                         gAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCatgcgg                                                                                                                                                                                                                     <  1:604029/186‑6 (MQ=255)
                                                         gAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCa                                                                                                                                                                    <  1:549520/235‑1 (MQ=255)
                                                                cccGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGAc                                                                                                                                                                                 <  1:814703/215‑1 (MQ=255)
                                                                cccGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCa                                                                                                                                                              <  1:477611/234‑1 (MQ=255)
                                                                   gCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTg                                                                                                                                                         <  1:12098/236‑1 (MQ=255)
                                                                    cAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTg                                                                                                                                                         <  1:522963/235‑1 (MQ=255)
                                                                           gcAGAAGCCAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGa                                                                                                                                                                                  <  1:864911/203‑1 (MQ=255)
                                                                             agaagaCAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCACCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGGGGGCTGGGCACGTCCGCCTTcc                                                                                                                                                                                                                                 >  2:690357/1‑154 (MQ=255)
                                                                                        caGCCCAGGCACTGGCCGACGTGTCGTCCGATGCCACCGAGCGGCAGCCGGCGCTGGCCGACGAGCAGGGTGACGAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGcc                                                                                                                                                                                     <  2:125749/186‑1 (MQ=255)
                                                                                          gCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGc                                                                                                                                    <  1:1034221/234‑1 (MQ=255)
                                                                                                        cGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGAAGGCGGACGTGCCCAGCCCTCGGCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTc                                                                                                                      >  1:994573/1‑234 (MQ=255)
                                                                                                         gAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTcggc                                                                                                                   >  1:674307/1‑236 (MQ=255)
                                                                                                                            ccGAGCGGCAGCGGGCGCTGGCCGACGAGCAGGGGGACCAGGTGGCCGGCGACTTCGACGACCTGATCAGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAggg                                                                                                                                                                                                                      <  2:53287/118‑1 (MQ=255)
                                                                                                                             cGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGAc                                                                                                                                                                                 <  1:247581/154‑1 (MQ=255)
                                                                                                                                          gcggTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGAc                                                                                                                                       <  1:150836/183‑1 (MQ=255)
                                                                                                                                             cTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTGGCTGCGGTGAAGAACATTGGTCGTTTg                                                                                    <  1:1067374/231‑1 (MQ=255)
                                                                                                                                              tGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGAAGGCGGACGTGCCCAGCCCTCGGCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTgg                                                                                                              <  2:994572/204‑1 (MQ=255)
                                                                                                                                                 ccGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGc                                                                                                                                                       <  2:65228/160‑1 (MQ=255)
                                                                                                                                                               gACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTGGCTGCGGTGAAGAACATTGGTCGTTTGACCGAGATGCAGGGTATGAAGcc                                                             <  1:681944/236‑1 (MQ=255)
                                                                                                                                                                     gTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAgtggtg                                                                                                                                         <  2:1075121/154‑1 (MQ=255)
                                                                                                                                                                     gTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGAGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTGGCTGCGGTGAAGAACATTGGTCGTTTGACCGAGATGCAGGGTATGAAGCCGAGCg                                                        <  1:700217/235‑1 (MQ=255)
                                                                                                                                                                             cGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGt                                                                                                                                >  2:247885/1‑155 (MQ=255)
                                                                                                                                                                             cGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACa                                                                                                                          >  2:914803/1‑161 (MQ=255)
                                                                                                                                                                             cGACTTCGACGACCTGATCCGCATGCGCGCGGAAGGCGGACGTGCCCAGCCCTCGGCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTGGCTGCGGTGAAGAACATTGGTCGTTTGACCGAGATGCAGGGTATGAAGCCGAGCGATTTCAAcc                                               >  1:934950/1‑236 (MQ=255)
                                                                                                                                                                               aCTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTg                                                                                                               >  2:324759/1‑170 (MQ=255)
                                                                                                                                                                                 ctCGGCGACCTGCTCCGCAGGCGCGCGGCGGAGGGCTGGGCACGTGCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAg                                                                                                                                                                                                   <  2:157734/83‑1 (MQ=255)
                                                                                                                                                                                  tCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGc                                                                                                                                    >  1:269021/1‑146 (MQ=255)
                                                                                                                                                                                     gcgacCTGAGCCGCAGGCGGGCGGCCGAGGGCTGGGCACGGCGGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTg                                                                                                                                                         <  2:23208/121‑1 (MQ=255)
                                                                                                                                                                                                    tGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTgg                                                                                                              >  2:26512/1‑150 (MQ=255)
                                                                                                                                                                                                      cgcgcgGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTGGCTGCGGTGAAGAACATTGGTCGTTTGa                                                                                   <  2:892630/175‑1 (MQ=255)
                                                                                                                                                                                                                       ggCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTGGCTGCGGTGAAGAACATTGGTCGTTTGACCGAGATGCAGGGTATGAAGCCGAGCGATTTCAACCAATATCGTTCGTTCAAGATGACATTTACACCTGAGGATCTa      <  1:492489/235‑1 (MQ=255)
                                                                                                                                                                                                                        gCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTGGCTGCGGTGAAGAACATTGGTCGTTTGACCGAGATGCAGGGTATGAAGCCGAGCGATTTCAACCAATATCGTTCGTTCAAGATGACATTTacac                 <  2:674306/223‑1 (MQ=255)
                                                                                                                                                                                                                        ccgtcatgcccTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTGGCTGCGGGGAAGAACATTGGTCGTTTGACCGAGATGCAGGGTATGAAGCCGAGCGATTTCAACCAatat                                          >  1:295218/10‑198 (MQ=255)
                                                                                                                                                                                                                          aCGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTGGCTGCGGTGAAGAACATTGGTCGTTTGACCGAGATGCAGGGTATGAAGCCGAGCGATTTCAAc                                                >  2:1127999/1‑190 (MQ=255)
                                                                                                                                                                                                                          aCGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTGGCTGCGGTGAAGAACATTGGTCGTTTGACCGAGATGCAGGGTATGAAGCCGAGCGATTTCAACCAATATCGTTCGTTCAAGATGACATTTACACCTGAGGATCTa      >  2:787091/1‑232 (MQ=255)
                                                                                                                                                                                                                          aCGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTGGCTGCGGTGAAGAACATTGGTCGTTTGACCGAGATGCAGGGTATGAAGCCGAGCGATTTCAACCAATATCGTTCGTTCAAGATGACATTTACACCTGAGGATCTAGCGt  >  1:511066/1‑236 (MQ=255)
                                                                                                                                                                                                                                   ttCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTGGCTGCGGTGAAGAACATTGGTCGTTTGACCGAGATGCAGGGTATGAAGCCGAGCGATTTCAACCAATATCGTTCGTTCAAGATGACATTTACACCTGAGGATCTAGc    <  1:344926/225‑1 (MQ=255)
                                                                                                                                                                                                                                   ttCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTGGCTGCGGTGAAGAACATTGGTCGTTTGACCGAGATGCAGGGTATGAAGCCGAGCGATTTCAACCAATATCGTTCGTTCAAGATGACATTTACACCTGAGGATCTAGCGt  <  1:449744/227‑1 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                  
GCAGCTGAAGAGGCCAAGAAAAAGGCCGCTGCAGATGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTGGCTGCGGTGAAGAACATTGGTCGTTTGACCGAGATGCAGGGTATGAAGCCGAGCGATTTCAACCAATATCGTTCGTTCAAGATGACATTTACACCTGAGGATCTAGCGT  >  AE015451/1397713‑1398166

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