Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 210,721 T→C 6.9% N5409N (AAT→AAC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451210,7210TC 6.9% 251.1 / 2.4 82N5409N (AAT→AACPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (38/38);  new base C (2/4);  total (40/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCG  >  AE015451/210501‑210956
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           
gCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAg                                                                                                                                                                                                                                       <  2:893282/227‑1 (MQ=255)
   gggCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGTCGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCa                                                                                                                                                                                                                             <  2:674557/232‑1 (MQ=255)
       ggACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCCACAATg                                                                                                                                                                                                                                            >  2:848320/1‑215 (MQ=255)
       ggACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCg                                                                                                                                                                                                                       >  1:116466/1‑236 (MQ=40)
        gACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCgg                                                                                                                                                                                                                      >  1:408179/1‑236 (MQ=255)
                 aaGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCga                                                                                                                                                                                                                   >  2:761285/1‑230 (MQ=40)
                   gTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAgcg                                                                                                                                                                                                           <  1:1271884/236‑1 (MQ=40)
                     tCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCa                                                                                                                                                                                                                             <  2:1255480/216‑1 (MQ=40)
                           aaCTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTa                                                                                                                                                                                                    >  1:486649/1‑235 (MQ=255)
                                aCGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACAccc                                                                                                                                                                                               <  1:921404/235‑1 (MQ=255)
                                      gACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCt                                                                                                                                                                                          <  1:605252/234‑1 (MQ=40)
                                       aCAAAACCCCGGTCAGCACGTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCa                                                                                                                                                                                                              <  1:1177663/213‑1 (MQ=18)
                                          aGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAg                                                                                                                                                                                    >  1:662859/1‑236 (MQ=40)
                                              gCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGACAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGTCCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGc                                                                                                                                                                                 <  1:798344/235‑1 (MQ=37)
                                                   tCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTgcgcctgtctcttat                                                                                                                                                                           >  1:910114/1‑225 (MQ=255)
                                                       cACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTa                                                                                                                                                                                                    <  2:910109/207‑1 (MQ=38)
                                                           tCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCCAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCg                                                                                                                                                                                <  2:773314/223‑1 (MQ=37)
                                                            cGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGAt                                                                                                                                                                              <  2:1096868/224‑1 (MQ=25)
                                                              gTCACCGACGAGCCAGGTACCCCAGGCAACGACGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCa                                                                                                                                                                                                              <  2:68663/190‑1 (MQ=37)
                                                                cACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGacg                                                                                                                                                              >  1:354491/1‑234 (MQ=18)
                                                                 aCCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGc                                                                                                                                                             <  1:488638/236‑1 (MQ=37)
                                                                         gCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGc                                                                                                                                                                                 <  1:741770/208‑1 (MQ=35)
                                                                          ccAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAg                                                                                                                                                                                                                                       <  2:731715/153‑1 (MQ=37)
                                                                                 ccccAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGc                                                                                                                                                >  2:819390/1‑233 (MQ=35)
                                                                                  cccGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGc                                                                                                                                                                                 <  2:739483/199‑1 (MQ=17)
                                                                                    cAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTa                                                                                                                                                                      <  1:364852/208‑1 (MQ=33)
                                                                                      ggCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTgcgc                                                                                                                                                                                      >  2:741765/1‑190 (MQ=33)
                                                                                          aCGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTg                                                                                                                                                                                         >  2:345928/1‑183 (MQ=18)
                                                                                           cGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACt                                                                                                                                   <  1:189129/236‑1 (MQ=33)
                                                                                            gAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTTCACGTCAACCAGCCGCAGGCTCACGCCCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCgg                                                                                                                                                                                                     >  1:401101/1‑168 (MQ=255)
                                                                                             aaGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAg                                                                                                                                                                                                                           >  2:1055352/1‑146 (MQ=32)
                                                                                                            aGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACt                                                                                                                  >  1:230569/1‑236 (MQ=33)
                                                                                                              gTCCCCATCCCCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCa                                                                                                                                                                                                              <  2:291862/142‑1 (MQ=17)
                                                                                                                  ccATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGt                                                                                                                            <  2:686805/220‑1 (MQ=12)
                                                                                                                     tCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGt                                                                                                                            >  2:612169/1‑217 (MQ=34)
                                                                                                                     tCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGAc                                                                                                          >  1:1004603/1‑235 (MQ=38)
                                                                                                                                    cTTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGc                                                                                                                      <  2:506561/208‑1 (MQ=14)
                                                                                                                                                agaACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCaa                                                                               >  1:108139/1‑235 (MQ=38)
                                                                                                                                                agaACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCaa                                                                               >  1:430622/1‑235 (MQ=38)
                                                                                                                                                  aaCGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCtt                                                                                                       <  1:122291/209‑1 (MQ=38)
                                                                                                                                                         aaCCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTgcg                                                                                                                                                                                       >  1:412929/1‑122 (MQ=32)
                                                                                                                                                                  tCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGAcc                                                                                                         >  2:1060511/1‑191 (MQ=37)
                                                                                                                                                                   cACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCCCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCaa                                                                  >  2:799541/1‑229 (MQ=39)
                                                                                                                                                                      tGTTCACGTCACCCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCGAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATc                                                                                                                                                   <  2:213951/145‑1 (MQ=17)
                                                                                                                                                                            cGTCAACCAGCCGCAGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTa                                                                                                                                                                      <  2:275369/120‑1 (MQ=17)
                                                                                                                                                                             gTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCAt                                                                                                                                        <  2:389842/149‑1 (MQ=34)
                                                                                                                                                                             gTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCtt                                                                                                       >  1:405919/1‑182 (MQ=37)
                                                                                                                                                                             gTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCt                                                                                                        <  2:409577/181‑1 (MQ=37)
                                                                                                                                                                               cAACCAGCCGCTGGCCCACGACCTGGGCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCaa                                                                  <  1:293690/217‑1 (MQ=255)
                                                                                                                                                                                accccgCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCg                                                                                                     <  2:349498/176‑1 (MQ=37)
                                                                                                                                                                                aaCCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACAccc                                                                                                                                                                                               >  1:906142/1‑91 (MQ=34)
                                                                                                                                                                                aaCCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACg                                                                                                                                                                 >  2:1129820/1‑121 (MQ=18)
                                                                                                                                                                                aaCCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGtc                                                                                                                           >  2:404946/1‑158 (MQ=18)
                                                                                                                                                                                aaCCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGGTTACAACGATGCTGGCCAGATCAGCCTggg                                                                                                                                           >  1:528677/1‑143 (MQ=17)
                                                                                                                                                                                    aGCCGCGGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAgtg                                          <  1:437543/236‑1 (MQ=38)
                                                                                                                                                                                     gCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCAc                                                                                                                                                                                             >  1:823504/1‑88 (MQ=35)
                                                                                                                                                                                     gCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCAc                                                                                                                                                                                             <  2:823499/88‑1 (MQ=33)
                                                                                                                                                                                     gCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAgtgg                                         >  1:1137806/1‑236 (MQ=39)
                                                                                                                                                                                     gCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACAccc                                                                                                                                                                                               <  2:546664/86‑1 (MQ=18)
                                                                                                                                                                                        gCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGGCCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGacac                                                       <  1:130159/219‑1 (MQ=39)
                                                                                                                                                                                         cTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGACGCTGGcc                                                                                                                                                        >  1:1186602/1‑121 (MQ=9)
                                                                                                                                                                                           ggCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGcc                                    <  1:201957/235‑1 (MQ=38)
                                                                                                                                                                                           ggCTCACGACCTGGTCGTGACCATGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCa                                   <  1:621914/236‑1 (MQ=38)
                                                                                                                                                                                            gCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCaa                                  <  1:810624/236‑1 (MQ=255)
                                                                                                                                                                                                  gACCTGGTCGTGACCCTGAGCAACACTGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCtt                                                                                                       <  1:393592/161‑1 (MQ=38)
                                                                                                                                                                                                         tCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACAccc                                                                                                                                                                                               >  2:719343/1‑66 (MQ=37)
                                                                                                                                                                                                         tCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGATGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGc                                                                                                                      >  1:1203869/1‑139 (MQ=21)
                                                                                                                                                                                                         tCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGATGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGc                                                                                                                      <  2:1203864/139‑1 (MQ=25)
                                                                                                                                                                                                           gTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGAccc                    >  1:1229114/1‑235 (MQ=255)
                                                                                                                                                                                                           gTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAgtggt                                        <  1:1052704/215‑1 (MQ=12)
                                                                                                                                                                                                               cccTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCCGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGggcgg                                                                                  <  1:587363/169‑1 (MQ=255)
                                                                                                                                                                                                                 cTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGAt                                                                                                                                                                              >  1:442851/1‑75 (MQ=34)
                                                                                                                                                                                                                 cTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCgg              >  1:631102/1‑235 (MQ=40)
                                                                                                                                                                                                                 cTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTCCAACTGGGCGGCac                                                                               >  2:922035/1‑169 (MQ=255)
                                                                                                                                                                                                                    aGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGcc                                    <  1:744860/210‑1 (MQ=40)
                                                                                                                                                                                                                     gCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGggc                                                                                    >  1:381688/1‑161 (MQ=40)
                                                                                                                                                                                                                     gCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGAcc         >  1:764621/1‑236 (MQ=255)
                                                                                                                                                                                                                     gCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCa                                                                                                                                                       <  2:818335/94‑1 (MQ=18)
                                                                                                                                                                                                                      caacaaTGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCg        >  1:816650/1‑236 (MQ=255)
                                                                                                                                                                                                                      caacaaTGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGGGGCCAAGCTGACTGCGACCCCGTCGGTCAcc         >  1:259710/1‑235 (MQ=255)
                                                                                                                                                                                                                       accaaCGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTg                                                                                                                 <  1:1022866/128‑1 (MQ=25)
                                                                                                                                                                                                                       aacaaTGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCtcgagct                                                                                        >  1:1109440/1‑149 (MQ=255)
                                                                                                                                                                                                                            tGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAggcg  <  1:899614/236‑1 (MQ=39)
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           
GCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCG  >  AE015451/210501‑210956

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