Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 200,545 G→C 6.2% A2017A (GCG→GCC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451200,5450GC 6.2% 203.2 / 2.9 62A2017A (GCG→GCCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (35/23);  new base C (2/2);  total (37/25)
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.98e-01

GACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGG  >  AE015451/200338‑200777
                                                                                                                                                                                                               |                                                                                                                                                                                                                                        
gACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTc                                                                                                                                                                                                               >  1:609142/1‑235 (MQ=34)
gACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACAAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCg                                                                                                                                                                                                              >  1:384738/1‑236 (MQ=25)
    gAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGCCGGCACCAACGCGCCGGCGGTCAATGCCATCGACGCGGGCAGCGGCGCa                                                                                                                                                                                                                              >  2:88406/1‑216 (MQ=21)
        gCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAaacctg                                                                                                                                                                                                      <  1:31142/236‑1 (MQ=14)
         cGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAaacctga                                                                                                                                                                                                     <  1:674441/236‑1 (MQ=34)
                aTCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCAACGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTgg                                                                                                                                                                                               <  1:271861/235‑1 (MQ=25)
                              cACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCACCGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGt                                                                                                                                                                                                                       >  1:215823/1‑197 (MQ=21)
                                cccTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACaa                                                                                                                                                                                           >  2:208083/1‑223 (MQ=35)
                                 ccTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCATCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGc                                                                                                                                                                              <  1:278722/235‑1 (MQ=37)
                                 ccTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCa                                                                                                                                                                             >  1:475123/1‑236 (MQ=17)
                                   tGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAaacctg                                                                                                                                                                                                      <  1:660299/209‑1 (MQ=35)
                                     aCCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGa                                                                                                                                                                                              >  1:324600/1‑215 (MQ=21)
                                        aacaaAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACTGCACCACCGTCGTCGTGGCAGTCAACAGCACCACTGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGAt                                                                                                                                                                                                                           >  2:404864/1‑183 (MQ=255)
                                        aacaaAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAg                                                                                                                                                                        >  2:264796/1‑234 (MQ=35)
                                             aGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGAcccc                                                                                                                                                                                     <  2:624006/216‑1 (MQ=18)
                                             aGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGAcccc                                                                                                                                                                                     >  1:624007/1‑216 (MQ=18)
                                                cGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACcgac                                                                                                                                                               <  1:264797/235‑1 (MQ=38)
                                                     tGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGcc                                                                                                                                                          <  1:539192/235‑1 (MQ=255)
                                                      gCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACcga                                                                                                                                                                >  2:627959/1‑228 (MQ=255)
                                                         gATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGca                                                                                                                                                     <  1:617547/236‑1 (MQ=40)
                                                          aTCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGca                                                                                                                                                     <  1:208083/235‑1 (MQ=38)
                                                             aacaacCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACAAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGcaca                                                                                                                                                   <  2:384737/234‑1 (MQ=255)
                                                                      tCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGa                                                                                                                                                                    <  2:324599/208‑1 (MQ=38)
                                                                         gAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAgg                                                                                                                                     <  2:376431/236‑1 (MQ=255)
                                                                            cTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAAc                                                                                                                                          >  2:497003/1‑228 (MQ=255)
                                                                             tGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGat                                                                                                                                 >  1:635145/1‑236 (MQ=255)
                                                                                  ttCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCg                                                                                                                             <  1:445199/235‑1 (MQ=255)
                                                                                       gCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCa                                                                                                                                                                             >  1:13760/1‑182 (MQ=37)
                                                                                          gACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTaa                                                                                                                                           >  1:629848/1‑213 (MQ=255)
                                                                                              gACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCActg                                                                                                                                                                          >  1:57559/1‑176 (MQ=38)
                                                                                              gACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGCCACCAt                                                                                                                >  1:203619/1‑236 (MQ=255)
                                                                                              gACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGCACCTGGACAAgc                                                                                                                                                                                         >  2:213168/1‑162 (MQ=255)
                                                                                               aCGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATc                                                                                                               >  1:559593/1‑236 (MQ=255)
                                                                                                             gtcgtGGCAGCCAACAGCACCACGGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAaacctg                                                                                                                                                                                                      <  1:101668/135‑1 (MQ=37)
                                                                                                                gtGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCa                                                                                                                                                                             >  2:49094/1‑157 (MQ=38)
                                                                                                                      gCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGa                                                                                                       >  1:531404/1‑221 (MQ=255)
                                                                                                                                 caccGGTTCGGCCACTGCAACCGCCCCAGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTc                                                                              <  1:592752/235‑1 (MQ=255)
                                                                                                                                 caccGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATTTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTc                                                                              <  1:60406/235‑1 (MQ=255)
                                                                                                                                  accGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCaa                                                                            <  1:497004/236‑1 (MQ=255)
                                                                                                                                      gTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCa                                                                                                                                                                             >  1:631452/1‑135 (MQ=255)
                                                                                                                                      gTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAAcc                                                                         >  1:102066/1‑235 (MQ=255)
                                                                                                                                        tCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTc                                                                              <  2:354950/228‑1 (MQ=255)
                                                                                                                                           gCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGCAACGAAGGCGac                                                                                                                                 >  1:496218/1‑173 (MQ=37)
                                                                                                                                                      cGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACg                                                        <  1:88406/236‑1 (MQ=255)
                                                                                                                                                             gACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAAcc                                                  >  1:24417/1‑235 (MQ=255)
                                                                                                                                                                    tCTATGGCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCt                                                                    <  1:233777/210‑1 (MQ=255)
                                                                                                                                                                     cTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCgc                                                                                                                                                                         >  2:71139/1‑107 (MQ=14)
                                                                                                                                                                     cTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGGCATCGGTAGCCGAGAACGTCAAACCGACCTTCAct                                                               >  1:65808/1‑213 (MQ=255)
                                                                                                                                                                      tATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCgg                                                                                            >  2:360964/1‑184 (MQ=18)
                                                                                                                                                                          tCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCAc                                    <  1:496290/236‑1 (MQ=255)
                                                                                                                                                                           cGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGccc                                      >  2:496289/1‑233 (MQ=255)
                                                                                                                                                                               aCCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACACCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACCGCCGACCAgg                                                                                                  >  2:223853/1‑169 (MQ=12)
                                                                                                                                                                                ccAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAg                                                                                                                        >  2:375512/1‑146 (MQ=255)
                                                                                                                                                                                ccAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAg                                                                                                                        <  1:375513/146‑1 (MQ=255)
                                                                                                                                                                                ccAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGta                                                            >  1:305963/1‑205 (MQ=255)
                                                                                                                                                                                ccAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTg                              >  1:304414/1‑236 (MQ=255)
                                                                                                                                                                                    cgcgCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCg                                             >  2:711120/1‑217 (MQ=255)
                                                                                                                                                                                      cgcCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTc                                                                              >  1:576077/1‑182 (MQ=255)
                                                                                                                                                                                      cgcCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTc                                                                              <  2:576076/182‑1 (MQ=255)
                                                                                                                                                                                       gcCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGCCCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTc                                                                              <  1:229280/181‑1 (MQ=255)
                                                                                                                                                                                                      cATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGGACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCa                                     >  2:21617/1‑207 (MQ=255)
                                                                                                                                                                                                      cATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACg        >  1:591977/1‑236 (MQ=255)
                                                                                                                                                                                                            cgcgGGCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAg                <  2:102066/222‑1 (MQ=255)
                                                                                                                                                                                                             gcgGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCa    <  2:635144/233‑1 (MQ=255)
                                                                                                                                                                                                             gcgGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAgg  <  1:578180/235‑1 (MQ=255)
                                                                                                                                                                                                               |                                                                                                                                                                                                                                        
GACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGG  >  AE015451/200338‑200777

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