Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 195,706 C→T 7.4% N404N (AAC→AAT PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451195,7060CT 7.4% 181.7 / 4.6 67N404N (AAC→AATPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (26/35);  new base T (3/2);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTC  >  AE015451/195484‑195921
                                                                                                                                                                                                                              |                                                                                                                                                                                                                       
cGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACAGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCCAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCGCCGCGACCCTGACGGCGAGCGGGTCGGTCACCTAAGGCGGCGAGATCACTTACACCGTGACCCTGAGCAACCCt                                                                                                                                                                                                                      <  1:547062/226‑1 (MQ=255)
   cGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCg                                                                                                                                                                                                          >  1:1040475/1‑235 (MQ=25)
   cGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGctgt                                                                                                                                                                                                         >  1:502198/1‑233 (MQ=255)
                     gaaTGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                                                                                           >  2:799146/2‑216 (MQ=9)
                              cTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGtgacag                                                                                                                                                                                                   <  2:903792/215‑6 (MQ=2)
                               tACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAAc                                                                                                                                                                             <  1:317725/236‑1 (MQ=12)
                                                 gTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                            >  1:627840/1‑235 (MQ=255)
                                                  acacagTCACCATCGAGAACGCCACAGGCGCTAATTTCGAGCAGCTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACATCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCTACGGCCAAACGATTACCGTTGAAGccg                                                                                                                                                            <  1:1068851/229‑3 (MQ=11)
                                                        tCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGccg                                                                                                                                                            <  2:582209/228‑3 (MQ=2)
                                                               cGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCt                                                                                                                                                                                            >  1:19047/1‑189 (MQ=25)
                                                                      gCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAAc                                                                                                                                                                                                     >  1:827003/1‑173 (MQ=12)
                                                                          ccGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTg                                                                                                                                                                                                 <  1:538606/173‑1 (MQ=255)
                                                                           aggcggTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGCCCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                            <  1:349008/208‑1 (MQ=255)
                                                                             gcggTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                               <  1:1091451/236‑1 (MQ=255)
                                                                                   aTTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGcc                                                                                                                                                                                                                   <  2:563807/146‑1 (MQ=17)
                                                                                   aTTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGcc                                                                                                                                                                                                                   >  1:563808/1‑146 (MQ=17)
                                                                                    tctCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACg                                                                                                                                                                                                             <  1:92690/149‑1 (MQ=9)
                                                                                          gCAGCTGACCCCGACTCCGACCCCGGCGCAGACCCCGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGcc                                                                                                                                                                                                           <  1:744391/147‑1 (MQ=255)
                                                                                          gCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGTTCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCACCCCTGCCCAGACGCCGGTAACCGTg                                                                                                                                                                                                 >  1:954116/1‑157 (MQ=255)
                                                                                                    ccGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGcc                                                                                                                                                                                                           <  1:327486/137‑1 (MQ=17)
                                                                                                       aaTCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGcc                                                                                                                                                                                 >  1:1078198/1‑160 (MQ=38)
                                                                                                           cGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACAATTACCGTTGa                                                                                                                                                                 <  1:712047/172‑1 (MQ=255)
                                                                                                           cGACCCCGGCTCAGACCACCATGAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGTGCCCGTCCGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTc                                                                                                                                                                                         <  1:517154/148‑1 (MQ=255)
                                                                                                             aCCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                            <  2:834135/175‑1 (MQ=255)
                                                                                                              ccccGGCTCAGACCCCGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAAc                                                                                                                                                                             <  1:439296/157‑1 (MQ=255)
                                                                                                              ccccGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACgg                                                                                                                                                                                   <  1:246302/151‑1 (MQ=38)
                                                                                                                                             cGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTc                                                                                                             >  2:501623/1‑190 (MQ=38)
                                                                                                                                             cGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCAc                                                                                >  1:105553/1‑219 (MQ=39)
                                                                                                                                               atgccaccaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGcc                                                                                                                                                                                                           <  2:1040762/91‑1 (MQ=21)
                                                                                                                                                             aCCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACACCGGCCAATGAc                                                                                                                >  1:94591/1‑171 (MQ=38)
                                                                                                                                                              cccTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGccg                                                                                                                                                            <  2:870510/126‑3 (MQ=2)
                                                                                                                                                              cccTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGccg                                                                                                                                                            >  1:870511/1‑124 (MQ=2)
                                                                                                                                                                    cGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGt                                                                                                                                                                                          <  2:264333/90‑1 (MQ=255)
                                                                                                                                                                           cccGTGGGTCACCGAAGGCGGCGTGATCACATACACCGTGACCCTGAGCAACCCTGCCCAGACGc                                                                                                                                                                                                            <  2:95239/65‑1 (MQ=255)
                                                                                                                                                                           cccGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCa                                                                                                                                                                                                                 <  1:725726/60‑1 (MQ=38)
                                                                                                                                                                             cgacGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAAc                                                                                                                                                                             <  2:557762/91‑1 (MQ=255)
                                                                                                                                                                             cGTGGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                 <  2:348698/106‑1 (MQ=255)
                                                                                                                                                                                cGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGtt                                                                                                                                                                   <  2:652725/101‑1 (MQ=255)
                                                                                                                                                                                 ggTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCGGACGCCGGTAACCGTGAGCCTGTCCAACGGCCa                                                                                                                                                                                <  2:225819/87‑1 (MQ=255)
                                                                                                                                                                                  gtgaCCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAAt                                                                                                                   <  2:48221/144‑1 (MQ=21)
                                                                                                                                                                                   tCACCGAAGGTGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGc                                                                                                                                                                                                            <  2:240167/57‑1 (MQ=16)
                                                                                                                                                                                         aaGGCGGCGTGATCACTTACACCGTGAGCCTGTGCAAATCTGCCAAGACGCCGGTAACCGTGAAACTGTGCAAGGGCTAAACGATTAGAGATGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGATCGACGGTCAGCGTCa                                                                                 >  1:769229/1‑174 (MQ=11)
                                                                                                                                                                                            gcggcgTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAAc                                                                                                                                                                                     >  1:48444/1‑71 (MQ=255)
                                                                                                                                                                                             aggcgTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGcc                                                                                                                                                                                                           <  2:13454/47‑1 (MQ=38)
                                                                                                                                                                                                cgTAATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTc                                                                                                                                                                                         <  2:563188/63‑1 (MQ=255)
                                                                                                                                                                                                     ccACGTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTggcg                                                            <  1:980770/182‑1 (MQ=255)
                                                                                                                                                                                                       aCTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACg      <  1:501624/235‑1 (MQ=255)
                                                                                                                                                                                                             acCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTggcgtca                                                         <  2:1080577/178‑4 (MQ=255)
                                                                                                                                                                                                                  gACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAAc                                                                                                                                                                             >  2:1028684/1‑57 (MQ=255)
                                                                                                                                                                                                                    cccTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGAtt                                                                                                                                    >  2:1037477/1‑96 (MQ=255)
                                                                                                                                                                                                                    cccTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCAcc                                                                               <  2:177855/149‑1 (MQ=37)
                                                                                                                                                                                                                    cccTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCAcc                                                                               >  1:177855/1‑149 (MQ=37)
                                                                                                                                                                                                                    cccTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAACGCCACCggc                                                             >  2:895494/1‑167 (MQ=255)
                                                                                                                                                                                                                    cccTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACAATCGAGAGCGCCACTGGCGGCAACTTCGAGCAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTc  >  1:987862/1‑226 (MQ=255)
                                                                                                                                                                                                                    cccTGAGCAACCCTGCCCAGACGCAGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCATTGACGACTACaacaac                                                                                                    >  1:302452/1‑128 (MQ=255)
                                                                                                                                                                                                                      cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTc                                                                                                                                                                                         >  2:913500/1‑41 (MQ=255)
                                                                                                                                                                                                                      cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTc                                                                                                                                                                                         <  1:913501/41‑1 (MQ=255)
                                                                                                                                                                                                                      cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCTACGGCAGCGTCGATTTCCAGACCCCGGCCAATTACACCTACAACAGCgc                                                                                                  >  1:1076058/1‑127 (MQ=11)
                                                                                                                                                                                                                      cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCgg                        >  2:1123014/1‑202 (MQ=255)
                                                                                                                                                                                                                      cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAACGCCACGggcggc                                                          >  2:468264/1‑168 (MQ=255)
                                                                                                                                                                                                                       cgAGCAACCCTGCCCAGACGCCGGTAACCGTGAGCCGGTCCAACGGCCAAACGATTTCCGtt                                                                                                                                                                   <  2:380739/61‑1 (MQ=14)
                                                                                                                                                                                                                          gCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGAccc                                                                                                                           >  1:397553/1‑99 (MQ=255)
                                                                                                                                                                                                                            aaCCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTacaaca                                                                                                      <  2:641975/118‑1 (MQ=12)
                                                                                                                                                                                                                            aaCCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTacaaca                                                                                                      >  1:641976/1‑118 (MQ=255)
                                                                                                                                                                                                                            aaCCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCaa                                                                                                                                                                               >  1:13454/1‑45 (MQ=255)
                                                                                                                                                                                                                            aaCCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGGCAGCGTCACCATCGAGAgcgc                                                                    >  1:343465/1‑152 (MQ=38)
                                                                                                                                                                                                                            aaCCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGCAAAACCCAGGGCAGCGTCGAttt                                                                                                                                   >  1:890804/1‑89 (MQ=255)
                                                                                                                                                                                                                              cccTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTgg                                                              <  2:409169/156‑1 (MQ=40)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                       
CGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTC  >  AE015451/195484‑195921

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