Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 195,424 G→A 11.5% V310V (GTG→GTA PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451195,4240GA11.5% 143.1 / 10.0 52V310V (GTG→GTAPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (25/21);  new base A (4/2);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

TCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCAC  >  AE015451/195219‑195636
                                                                                                                                                                                                             |                                                                                                                                                                                                                    
tCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGCCGCCGGTGACCGTGACCCTGTCCAACGGCCAGGtcat                                                                                                                                                                                         >  1:23882/1‑235 (MQ=255)
     ggTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGcc                                                                                                                                                                                                >  1:679767/1‑223 (MQ=255)
      gTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGgtgaccgtgacc                                                                                                                                                                                                             >  2:1124045/1‑209 (MQ=255)
       tGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGgtgac                                                                                                                                                                                                                    >  1:808271/1‑201 (MQ=255)
                           tGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAAc                                                                                                                                                                                                    <  1:800376/197‑1 (MQ=255)
                             ccAGGTAGGACCGGTGGGGGCGGGCTGACGGTGCCATGGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGCCCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACgg                                                                                                                                                                                                  <  1:231175/196‑1 (MQ=255)
                                        cGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGGCCAACGGCCAgg                                                                                                                                                                                             >  1:222786/1‑191 (MQ=255)
                                          ggggtggCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                              <  1:711610/234‑1 (MQ=255)
                                             gtggCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGAcc                                                                                                                                           >  1:336190/1‑236 (MQ=255)
                                                  gggCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAAccccg                                                                                                                                                                >  1:886306/1‑208 (MQ=255)
                                                     cTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAACTGGCCGATACCATCGACCCTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCg                                                                                                                                                                           <  1:893715/196‑1 (MQ=255)
                                                            tGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGcc                                                                                                                                                                                                >  1:652726/1‑168 (MQ=255)
                                                                                        ccAGGTGACCGATACAATCGACACTACAACCGCGACCCTGACGGCGTCCCCATCCGTCACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAAc                                                                                                                                                                                                    <  2:823132/136‑1 (MQ=11)
                                                                                              gACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAggg                                                                                                                                                              >  1:688591/1‑168 (MQ=14)
                                                                                               aCCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCa                                                                                                                                                                                               <  1:975384/134‑1 (MQ=255)
                                                                                                 cGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCCGGTCATCACTGTTGAAGCCGGc                                                                                                                                                                         >  1:505480/1‑154 (MQ=255)
                                                                                                 cGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGCCAc                                                                                               >  2:28594/1‑228 (MQ=255)
                                                                                                 cGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCAc                                                                                               >  1:99624/1‑228 (MQ=255)
                                                                                                   aTACAATGGACACGACAACCGCGACCCTGACCGCCTCGCCAGGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAAc                                                                                                                                                                                                    <  2:35227/125‑1 (MQ=255)
                                                                                                      cAATCGCCCCTACAACCGCGCCCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCg                                                                                                                                                                           <  1:235236/147‑1 (MQ=255)
                                                                                                            acacTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGcc                                                                                                                                                                                                >  2:587740/1‑120 (MQ=255)
                                                                                                                tACACCCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAACCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCa                                                                                                <  1:226184/212‑1 (MQ=255)
                                                                                                                          aCCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCg                                                                                                                                                                           <  1:612352/127‑1 (MQ=255)
                                                                                                                          aCCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCg                                                                                                                                                                           >  2:612351/1‑127 (MQ=255)
                                                                                                                           cccTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGCCCAAc                                                                                                                                                                                                    >  1:625983/1‑101 (MQ=255)
                                                                                                                           cccTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGGCTACAACAACGGCTCGACCg                                                                                                         >  1:154182/1‑192 (MQ=255)
                                                                                                                                  cGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACGCCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCg                                                                                                                                                                           <  2:901012/119‑1 (MQ=255)
                                                                                                                                  cGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCgg                                                                                                                                       <  2:336189/155‑1 (MQ=255)
                                                                                                                                   gCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGGCCAACGGCCAGGTCATCACTGTTGAAGCGGGCAAAACCCAGGGc                                                                                                                                                             >  1:348699/1‑132 (MQ=255)
                                                                                                                                                tGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCg                                        <  1:954571/236‑1 (MQ=255)
                                                                                                                                                     gAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCATGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCg                                                                                                         >  1:312485/1‑166 (MQ=38)
                                                                                                                                                       aGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGa                                 <  1:1124046/236‑1 (MQ=255)
                                                                                                                                                         gcggCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAAc                                                                                                                                                                                                    <  2:48444/69‑1 (MQ=9)
                                                                                                                                                              gTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCg                                              <  1:333672/216‑1 (MQ=255)
                                                                                                                                                                 aTCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCg                                        >  1:16880/1‑219 (MQ=255)
                                                                                                                                                                         cacCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCg                <  1:25632/235‑1 (MQ=255)
                                                                                                                                                                            cGTGACCCTGAGCAACCCTGCCCAGACGCCGgtga                                                                                                                                                                                                                     <  2:517298/35‑1 (MQ=37)
                                                                                                                                                                               gACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAACCCGAcgc                                          >  1:105785/1‑201 (MQ=37)
                                                                                                                                                                                aCCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCg                                              <  2:198957/198‑1 (MQ=255)
                                                                                                                                                                                 cccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCa                                                                                                                                                                 >  1:100942/1‑82 (MQ=38)
                                                                                                                                                                                 cccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCAAg                                                                                                                                              <  2:924233/101‑1 (MQ=38)
                                                                                                                                                                                 cccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCAAg                                                                                                                                              >  1:924234/1‑101 (MQ=25)
                                                                                                                                                                                   cTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACgg                                                                                                                                                                                                  >  1:314862/1‑47 (MQ=255)
                                                                                                                                                                                      aGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAtgccaccac  <  1:28594/236‑1 (MQ=255)
                                                                                                                                                                                       gCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGaa                                                  >  2:725885/1‑187 (MQ=255)
                                                                                                                                                                                       gCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGaa                                                  <  1:725886/187‑1 (MQ=255)
                                                                                                                                                                                         aaCCCTGCCCAGACGCCGGTGACCGTGACCCTGTc                                                                                                                                                                                                        >  1:517299/1‑35 (MQ=255)
                                                                                                                                                                                         aaCCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGAccc                                          >  1:998646/1‑193 (MQ=255)
                                                                                                                                                                                           cccTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCg                <  1:1060068/217‑1 (MQ=255)
                                                                                                                                                                                            ccTGCCCAGACGCCGGGGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGaa                                                  <  2:318409/182‑1 (MQ=39)
                                                                                                                                                                                              tGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCTGAccc                                                                                                                                          >  1:625422/1‑92 (MQ=18)
                                                                                                                                                                                              tGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGaa                                                                                      >  1:827253/1‑144 (MQ=17)
                                                                                                                                                                                                        ccGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAATTTCGAGCAACTGa                                                         >  1:409170/1‑163 (MQ=37)
                                                                                                                                                                                                        ccGGTGACCGTGACCCTGTCCAACGGCCAGACCATCACTGTTGAAGCCGTCAAAACCCAGAGCAGCGTCGATTTCCAGACCCCTGCCAATGACGTCTACAACAACga                                                                                                                 >  1:1123773/1‑106 (MQ=255)
                                                                                                                                                                                                        ccGGTGACCGTCACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGTCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAACTTCCAGCAACTGACCCCGa                                                   >  1:557698/1‑169 (MQ=255)
                                                                                                                                                                                                             gaccgtgaccCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAcacc       >  2:1033400/1‑208 (MQ=17)
                                                                                                                                                                                                             |                                                                                                                                                                                                                    
TCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCAC  >  AE015451/195219‑195636

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