Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 653,157 Δ1 bp 8.2% coding (1586/1962 nt) ctpL ← Methyl‑accepting chemotaxis protein CtpL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451653,1570G. 8.2% 162.1 / 8.7 49coding (1586/1962 nt)ctpLMethyl‑accepting chemotaxis protein CtpL
Reads supporting (aligned to +/- strand):  ref base G (21/24);  new base . (2/2);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.82e-02

CTGTGCCGTTACATCGGCAATCCGCACCGCTGTATCGGAAATGGTGCGAATGGCACAGACAATCTCATCCAGTGCCCCATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCTTCGATCGCGGCATTCAGTGCCAGCAGGTTGGTCTGCTCGGCAATCGCGCGGATCACCGTCAGTACCCCGCCAATGGTGGCCGACTCTTCCGCCAATTGTTCGATCATGCGGGCGTTGCCTTGCACTTCATCGACCAGCTCA  >  AE015451/652940‑653363
                                                                                                                                                                                                                          |                                                                                                                                                                                                               
ctgTGCCGTTACATCGGCAATCCGCACCGCTGTATCGGAAATGGTGCGAATGGCACAGACAATCTCATCCAGTGCCCCATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAgcg                                                                                                                                                                                                      <  1:639138/228‑1 (MQ=255)
       gTTACATCGGCAATCCGCACCGCTGTATCGGAAATGGTGCGAATGGCACAGACAATCTCATCCAGTTCCCCATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCg                                                                                                                                                                                         <  2:489668/234‑1 (MQ=255)
        ttACATCGGCAATCCGCACCGCTGTATCGGGAAATGGTGCGAATGGCACAGACAATCTCATCCAGTGCCCCATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATcctgtct                                                                                                                                                                                       >  1:905707/1‑230 (MQ=255)
            aTCGGCAATCCGCACCGCTGTATCGGAAATGGTGCGAATGGCACAGACAATCTCATCCAGTGCCCCATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACg                                                                                                                                                                                   <  1:721690/235‑1 (MQ=255)
                 cAATCCGCACCGCTGTATCGGAAATGGTGCGAATGGCACAGACAATCTCATCCAGTGCCCCATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAg                                                                                                                                                                             <  1:303220/236‑1 (MQ=255)
                   aTCCGCACCGCTGTATCGGAAATGGTGCGAATGGCACAGACAATCTCATCCAGTGCCCCATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGc                                                                                                                                                                            >  1:569611/1‑235 (MQ=255)
                                           ggTGCGAATGGCACAGACAATCTCATCCAGTGCCCCATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGCCGGCACggg                                                                                                                                                    >  1:115849/1‑236 (MQ=255)
                                            gTGCGAATGGCACAGACAATCTCATCCAGTGCCCCATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCT‑‑‑‑CACTG‑‑‑G‑‑‑CTCAGCGCACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCTTCGAt                                                                                                                                          >  2:317781/1‑233 (MQ=255)
                                              gCGAATGGCACAGACAATCTCATCCAGTGCCCCATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCCCCTCATCGGCGACCACgg                                                                                                                                                                                  >  1:596067/1‑203 (MQ=255)
                                              gCGAATGGCACAGACAATCTCATCCAGTGCCCCATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCt                                                                                                                                                                                                 >  2:597347/1‑188 (MQ=255)
                                              gCGAATGGCACAGACAATCTCATCCAGTGCCCCATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGcc                                                                                                                                                  >  1:814869/1‑235 (MQ=255)
                                                   tGGCACAGACAATCTCATCCAGTGCCCCATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACCGCCCATCTCGCCGGCACGGGCCGCTTc                                                                                                                                             >  1:407549/1‑235 (MQ=255)
                                                      cacaGACAATCTCATCCAGTGCCCCATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCTTCGATc                                                                                                                                         <  1:397610/236‑1 (MQ=255)
                                                                 ccATCCAGTGCCCCATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGt                                                                                                                                                                                                           <  2:905704/158‑1 (MQ=255)
                                                                  cATCCAGTGCCCCATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTCACTGGCTCAGCGCACCTCATCGGCGACCACGGCAAAGCCAc                                                                                                                                                                         <  2:323104/192‑1 (MQ=255)
                                                                  cATCCAGTGCCCCATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTCACTGGCTCAGCGCACCTCATCGGCGACCACGGCAAAGCCAc                                                                                                                                                                         >  1:323107/1‑192 (MQ=255)
                                                                   agccAGTGCCCCATCCGCCGCCTGGGCCTGGTGGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCtt                                                                                                                                              <  1:69013/216‑1 (MQ=255)
                                                                   aTCCAGTGCCCCATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCT‑‑‑‑CACTG‑‑‑G‑‑‑CTCAGCGCACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCTTCGATCGCGGCATTCAGTGCCAGCAGGTTg                                                                                                                 <  1:317784/235‑1 (MQ=255)
                                                                    tCCAGTGCCCCATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGc                                                                                                                                                                                          <  1:133220/173‑1 (MQ=255)
                                                                       aGTGCCCCATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCTTCGATCGCGGCATTCAGTGCcag                                                                                                                        >  1:41279/1‑236 (MQ=255)
                                                                        gtaccccATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCCCCGGTAGCCCCCGTGGTGCGCTGAGCCAGGGAGCGCACCTCATCGGCg                                                                                                                                                                                         >  1:91928/4‑170 (MQ=255)
                                                                           ccccATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCTTCGATCGCGGCATTCAGTg                                                                                                                            <  1:299624/228‑1 (MQ=255)
                                                                            cccATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGc                                                                                                                                                                                          <  2:901382/165‑1 (MQ=255)
                                                                            cccATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGc                                                                                                                                                                                          >  1:901385/1‑165 (MQ=255)
                                                                            cccATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCTTCGATCGCGGCATTCAGTg                                                                                                                            <  1:157510/227‑1 (MQ=255)
                                                                             ccATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCa                                                                                                                                                                          <  1:187284/180‑1 (MQ=255)
                                                                                        ctgggcctggTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCTTCGATCGCGGCATTCAGTGCCAGCAGGTTGGTCTGCTCg                                                                                                        >  1:917812/1‑235 (MQ=255)
                                                                                                                         cAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCTTCGATCGCGGcc                                                                                                                                   >  1:140044/1‑174 (MQ=255)
                                                                                                                             tGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGcc                                                                                                                                                            <  2:791073/146‑1 (MQ=255)
                                                                                                                              gcgcgcgcATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCTTCGATCGCGGCATTCAGTGCCAGCAGGTTGGTCTGCTCGGCAATCGCGCGGATCACCGTCAGTACCCCGCCAAtggtg                                                                 >  1:642487/1‑236 (MQ=255)
                                                                                                                                 cgcgcAGCTCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGAGGTCGATCAACGCCAGGCTCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCaa                                                                                                                                                                               <  2:747582/123‑1 (MQ=39)
                                                                                                                                   cgcATCCCCGCCACCGACCTTCTGGCGGCCTGCTGCAGGCGGTCGATCACCGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCAtc                                                                                                                                                                 <  1:173961/135‑1 (MQ=255)
                                                                                                                                                 cGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCTTCGATCGCGGCATTCAGTGCCAGCAGGTTGGTCTGCTCGGCAATCGCGCGGATCACCGTCAGTACCCCGCCAATGGTGGCCGACTCTTCCGCCAAtt                                              <  1:27136/236‑1 (MQ=255)
                                                                                                                                                     cTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCTTCGATCGCGGCAtt                                                                                                                                 >  2:445703/1‑149 (MQ=255)
                                                                                                                                                         tcgcggcCTGCTGCAGGCGGTCGATCACCGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCTTCGATCGCGGCATTCAGTGCCAGCAGGTTGGTCTGCTCGGCAATc                                                                                                  <  2:116290/174‑1 (MQ=255)
                                                                                                                                                                        ggCGGTCGATCAACGCCTGGATCTCGCCAGTAGCCCCCGTGGTGCGCTGAGCCAGt                                                                                                                                                                                                           <  2:479624/56‑1 (MQ=255)
                                                                                                                                                                            gTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCTTCGATCGCGGCATTCAGTGCCAGCAGGTTGGTCTGCTCGGCAATcgcgcg                                                                                             >  1:210681/1‑162 (MQ=255)
                                                                                                                                                                                aTCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCTTCGATCGCGGCATTCAGTGCCAGCAGGTTGGTCTGCTCGGCAATCGCGCGGATCACCGTCAGTACCCCGCCAATGGTGGCCGACTCTTCCGCCAATTGTTCGATCATGCGGGCGTTGCCTTGCACtt                >  1:636623/1‑235 (MQ=255)
                                                                                                                                                                                  cAACGCCTGGATCTCGCCGGTAGCGCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGAGCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCTTCGATCGCGGCATTCAGTGCCAGCAGGTTGGTCTGCTCGGCAAt                                                                                                   <  1:714837/150‑1 (MQ=255)
                                                                                                                                                                                      gCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCTTCGATCGCGGCATTCAGTGCCAGCAGGTTGGTCTGCTCGGCAATCGCGCGGATCACCGTCAGTACCCCGCCAATGGTGGCCGACTCTTCCGCCAATTGTTCGATCATGCGGGCGTTGCCTTGCACTTCATc            <  2:642484/233‑1 (MQ=255)
                                                                                                                                                                                             tctcGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCTTCGATCGCGGCATTCAGTGCCAGCAGGTTGGTCTGCTCGGCAATCGCGCGGATCACCGTCAGTACCCCGCCAATGGTGGCCGACTCTTCCGCCAATTGTTCGATCATGCGGGCGTTGCCTTGCACTTCATCGACCAGCTCa  >  1:425924/1‑236 (MQ=255)
                                                                                                                                                                                              ctcGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCTTCGATCGCGGCATTCAGTGCCAGCAGGTTGGTCTGCTCGGCAATCGCGCGGATCACCGTCAg                                                                                 >  2:387366/1‑156 (MQ=255)
                                                                                                                                                                                              ctcGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCTTCGATCGCGGCATTCAGTGCCAGCAGGTTGGTCTGCTCGGCAATCGCGCGGATCACCGTCAg                                                                                 <  1:387369/156‑1 (MQ=255)
                                                                                                                                                                                              ctcGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCTTCGATCGCGGCATTCAGTGCCAGCAGGTTGGTCTGCTCGGCAATCGCGCGGATCACCGTCAGTACCCCGCCAATGGTGGCCGACTCTTCCGCCAATTGTTCGATCATGCGGGCGTTGCCTTGCACTTCATc            <  1:557290/225‑1 (MQ=255)
                                                                                                                                                                                              ctcGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCTTCGATCGCGGCATTCAGTGCCAGCAGGTTGGTCTGCTCGGCAATCGCGCGGATCACCGTCAGTACCCCGCCAATGGTGGCCGACTCTTCCGCCAATTGTTCGATCATGCGGGCGTTGCCTTGCACTTCATc            >  2:557287/1‑225 (MQ=255)
                                                                                                                                                                                               tcGCCGGTAGCCCCCGTGGTTCGCTGAGCCAGTGAGCGCACCTCTTCGTCTACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCTTCGATCGCGGCATTCAGTGCcagc                                                                                                                       <  2:832504/117‑1 (MQ=38)
                                                                                                                                                                                               tcGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGa                                                                                                                                                                                                         <  2:146366/35‑1 (MQ=255)
                                                                                                                                                                                                         cccccGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCCTCGATCGCGGCATTCAGTGCCAGCAGGTTGGTCTGCTCGGCAATCGCGCGGATCAcc                                                                                      >  1:895668/1‑140 (MQ=255)
                                                                                                                                                                                                                       tGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCTTCGATCGCGGCATTCAGTGCCAGCAGGTTGGTCTGCTCGGCAATCGCGCGGATCACCGTCAGTACCCCGCCAATGGTGGCCGACTCTTCCGCCAATTGTTCGATCATGCGGGCGTTGCCTTGCACTTCATCGa          >  2:437709/1‑202 (MQ=255)
                                                                                                                                                                                                                          |                                                                                                                                                                                                               
CTGTGCCGTTACATCGGCAATCCGCACCGCTGTATCGGAAATGGTGCGAATGGCACAGACAATCTCATCCAGTGCCCCATCCGCCGCCTGGGCCTGGTTGGCCGTGGCCTCGGCATGTTCCAGTTGCGCGCGCATCCCCGCCACCGAACTTCTGGCGGCCTGCTGCAGGCGGTCGATCAACGCCTGGATCTCGCCGGTAGCCCCCGTGGTGCGCTGAGCCAGTGAGCGCACCTCATCGGCGACCACGGCAAAGCCACGGCCCATCTCGCCGGCACGGGCCGCTTCGATCGCGGCATTCAGTGCCAGCAGGTTGGTCTGCTCGGCAATCGCGCGGATCACCGTCAGTACCCCGCCAATGGTGGCCGACTCTTCCGCCAATTGTTCGATCATGCGGGCGTTGCCTTGCACTTCATCGACCAGCTCA  >  AE015451/652940‑653363

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