Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 1,397,918 C→A 10.2% A289D (GCC→GAC) ‡ 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,9180CA10.2% 154.6 / 9.1 49A289D (GCC→GAC) ‡tolAcolicin S4 and filamentous phage transport system
Reads supporting (aligned to +/- strand):  ref base C (14/30);  new base A (3/2);  total (17/32)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

GCAGCTGAAGAGGCCAAGAAAAAGGCCGCTGCAGATGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTGGCTGCGGTGAAGAACATTGGTCGTTTGACCGAGATGCAGGGTATGAAGCCGAGCGATTTCAACC  >  AE015451/1397713‑1398121
                                                                                                                                                                                                             |                                                                                                                                                                                                           
gCAGCTGAAGAGGCCAAGAAAAAGGCCGCTGCAGATGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCcgcgcg                                                                                                                                                                               <  1:451834/236‑1 (MQ=255)
 cAGCTGAAGAGGCCAAGAAAAAGGCCGCTGCAGATGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGAATGCGCGCGGCCGAgg                                                                                                                                                                                                        >  1:275163/1‑210 (MQ=255)
      gAAGAGGCCAAGAAAAAGGCCGCTGCAGATGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAg                                                                                                                                                                           <  1:725464/234‑1 (MQ=255)
       aaGAGGCCAAGAAAAAGGCCGCTGCAGATGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAggg                                                                                                                                                                         <  1:673754/235‑1 (MQ=255)
        agagGCCAAGAAAAAGGCCGCTGCAGATGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGAGCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTg                                                                                                                                                                                                    >  1:789357/1‑207 (MQ=255)
                         ccGCTGCAGATGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGGGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGGGGGCTggtc                                                                                                                                                                                                 >  1:240548/1‑191 (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)
                                                                                                                                                                                                             |                                                                                                                                                                                                           
GCAGCTGAAGAGGCCAAGAAAAAGGCCGCTGCAGATGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTGGCTGCGGTGAAGAACATTGGTCGTTTGACCGAGATGCAGGGTATGAAGCCGAGCGATTTCAACC  >  AE015451/1397713‑1398121

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