Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 200,113 T→C 7.8% N1873N (AAT→AAC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451200,1130TC 7.8% 202.1 / 3.0 62N1873N (AAT→AACPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (22/35);  new base C (2/3);  total (24/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

CATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGAC  >  AE015451/199879‑200340
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                   
cATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATg                                                                                                                                                                                                                                    <  1:215846/236‑1 (MQ=38)
     aCGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAg                                                                                                                                                                                                                               >  1:465658/1‑236 (MQ=40)
                  cgcAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAg                                                                                                                                                                                                                               <  1:719007/223‑1 (MQ=21)
                  cgcAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAg                                                                                                                                                                                                                               >  2:719005/1‑223 (MQ=11)
                    cAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCt                                                                                                                                                                                                                <  1:162277/236‑1 (MQ=40)
                               aaGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGaaa                                                                                                                                                                                                         >  2:311179/1‑232 (MQ=255)
                               aaGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCa                                                                                                                                                                                                      <  2:481379/235‑1 (MQ=39)
                                aGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCa                                                                                                                                                                                                      <  2:689329/234‑1 (MQ=255)
                                aGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAgc                                                                                                                                                                                                    <  1:648447/236‑1 (MQ=255)
                                aGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAgc                                                                                                                                                                                                    <  1:28447/236‑1 (MQ=38)
                                  tttGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAgcg                                                                                                                                                                                                   <  1:596877/235‑1 (MQ=255)
                                                cctgGACAAGACCCCGGGCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAg                                                                                                                                                                                                                               <  1:573819/193‑1 (MQ=255)
                                                           ccccGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAgcgcgc                                                                                                                                                                                                >  2:324042/1‑213 (MQ=255)
                                                                   accACCGAGGTCCCCGCCGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAg                                                                                                                                                                           <  1:155665/224‑1 (MQ=255)
                                                                    gCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAgcg                                                                                                                                                                                                   <  2:719262/201‑1 (MQ=255)
                                                                       ccGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGCACACCCACGATGCGCAAGGCGATGACGTCTAt                                                                                                                                                             >  1:201756/1‑236 (MQ=255)
                                                                        cGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCCGGCGAAACCAgcg                                                                                                                                                                                                   >  2:537256/1‑197 (MQ=255)
                                                                            gTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAgg                                                                                                                                                         <  1:537354/235‑1 (MQ=255)
                                                                                 cgacgaGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCAAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAgggtca                                                                                                                                                                                                                          <  1:79764/165‑5 (MQ=255)
                                                                                 cgacgaGCCAGGCACGCCAGGCACCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAg                                                                                                                                                                                                                               <  1:191225/160‑1 (MQ=255)
                                                                                   acgaGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCa                                                                                                                                                                                                                                <  1:692/157‑1 (MQ=40)
                                                                                     gaGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGTGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGcc                                                                                                                                                <  1:434578/235‑1 (MQ=255)
                                                                                      aGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGcc                                                                                                                                                >  1:751260/1‑234 (MQ=255)
                                                                                               cTCCAGGCAACGAAGGCGATATCGTCAAAGTCACCATTACTGCTGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATg                                                                                                                                                                                                                                    <  2:315046/141‑1 (MQ=11)
                                                                                                tCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAgg                                                                                                                                                         <  2:376427/215‑1 (MQ=40)
                                                                                                         cGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACt                                                                                                                           <  1:543314/236‑1 (MQ=255)
                                                                                                                  cagcGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGGCAAGCCGACCGTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAGTGCTCAGGTCACCATCAAGGCTGGCGAAACCAgcgcg                                                                                                                                                                                                 <  2:54004/154‑1 (MQ=255)
                                                                                                                  cATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAg                                                                                                                  >  1:347430/1‑236 (MQ=255)
                                                                                                                      gTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGc                                                                                                              <  1:324044/236‑1 (MQ=255)
                                                                                                                      gTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACAcc                                                                                                                                                                                        >  2:776036/1‑162 (MQ=14)
                                                                                                                      gTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACAcc                                                                                                                                                                                        <  1:776038/162‑1 (MQ=14)
                                                                                                                         aaGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCa                                                                                                                      <  1:617829/225‑1 (MQ=255)
                                                                                                                         aaGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAgcgcg                                                                                                                                                                                                 <  1:1822/150‑1 (MQ=17)
                                                                                                                                  aTCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCa                                                                                                                      >  1:717743/1‑216 (MQ=255)
                                                                                                                                  aTCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCa                                                                                                                      <  2:717741/216‑1 (MQ=255)
                                                                                                                                  aTCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGAAGCGCAAGGCGATGAc                                                                                                                                                                   >  2:163522/1‑171 (MQ=25)
                                                                                                                                    cACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAACCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTCGTCGTGCCCCTGAGCAACAATGCTCAGGTCACCATCCAGGCTGGCGACACCAg                                                                                                                                                                                                     >  2:106634/1‑135 (MQ=11)
                                                                                                                                         ccGACCAGACTTAGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGACCAGATCAGCCTGGGCATCAACTCGGCAGTa                                                                                                                   >  2:447734/1‑212 (MQ=255)
                                                                                                                                            aCCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCgg                                                                                                                        <  1:86099/204‑1 (MQ=255)
                                                                                                                                            aCCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGtgc                                                                                                     >  2:559996/1‑223 (MQ=255)
                                                                                                                                              cAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTg                                                                                                         >  1:424098/1‑217 (MQ=255)
                                                                                                                                              cAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTg                                                                                                         <  2:424096/217‑1 (MQ=255)
                                                                                                                                                   gtCGGTGGCCGAGACCGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGa                                                                                                                 <  1:258691/203‑1 (MQ=255)
                                                                                                                                                   gtCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTg                                                                                                                                     <  2:736174/183‑1 (MQ=38)
                                                                                                                                                              agaACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAg                                                                                                                                                                           <  1:517695/135‑1 (MQ=21)
                                                                                                                                                                        aCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGtt                                                             >  1:790556/1‑235 (MQ=255)
                                                                                                                                                                                tCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTc                                                       >  2:678456/1‑233 (MQ=2)
                                                                                                                                                                                     gTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGacac                                               <  1:311181/236‑1 (MQ=255)
                                                                                                                                                                                         aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAg                                                                                                                                                                                                                               <  2:132025/56‑1 (MQ=35)
                                                                                                                                                                                         aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGTTGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCt                                            >  1:752602/1‑235 (MQ=255)
                                                                                                                                                                                           gTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCa                                                                                                                                                                                                                                <  2:783687/53‑1 (MQ=35)
                                                                                                                                                                                            tCAACCAGCCGCTGGCCCACGACCTGGTCGGGACCCTGAGCAACAATGCGCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGa                                                  <  2:23198/226‑1 (MQ=255)
                                                                                                                                                                                              aaCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTa                                                                                                                   >  2:519037/1‑159 (MQ=255)
                                                                                                                                                                                              aaCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTa                                                                                                                   <  1:519039/159‑1 (MQ=255)
                                                                                                                                                                                                        tGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCGCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGa                                                  >  2:447046/1‑214 (MQ=255)
                                                                                                                                                                                                                     ggTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCa                                                                                                                      <  1:534918/133‑1 (MQ=255)
                                                                                                                                                                                                                         gTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCAccc            >  1:293008/1‑235 (MQ=255)
                                                                                                                                                                                                                            aCCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAgc                                                                                                                                                                                                    <  2:216002/48‑1 (MQ=38)
                                                                                                                                                                                                                                   gCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAgtgg                                 >  1:692802/1‑204 (MQ=255)
                                                                                                                                                                                                                                    caacaaTGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCt                                                                                <  1:773675/156‑1 (MQ=255)
                                                                                                                                                                                                                                    caacaaTGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCATTCGAGAACCTTGAGCt                                                                                >  2:773673/1‑156 (MQ=255)
                                                                                                                                                                                                                                     aacaaTGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGAc  >  2:452959/1‑233 (MQ=255)
                                                                                                                                                                                                                                          tGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAg                                                                                                                     >  2:711100/1‑113 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                   
CATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGAC  >  AE015451/199879‑200340

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