Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 213,739 T→C 6.7% G6415G (GGT→GGC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451213,7390TC 6.7% 280.4 / 4.6 85G6415G (GGT→GGCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (46/33);  new base C (2/4);  total (48/37)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

CGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAAC  >  AE015451/213504‑213971
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                        
cGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGt                                                                                                                                                                                                                                          >  1:1006622/1‑236 (MQ=255)
cGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGt                                                                                                                                                                                                                                          <  1:1155627/236‑1 (MQ=40)
 gACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGt                                                                                                                                                                                                                                          <  1:1303737/235‑1 (MQ=40)
        gagaACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCAc                                                                                                                                                                                                                                  >  1:1277388/1‑236 (MQ=255)
                           gcAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACaa                                                                                                                                                                                                                >  1:1032518/1‑235 (MQ=255)
                                gCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGCCCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTa                                                                                                                                                                                                          <  1:273100/236‑1 (MQ=255)
                                    cGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGt                                                                                                                                                                                                       >  1:1033897/1‑235 (MQ=255)
                                    cGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTc                                                                                                                                                                                                      >  1:1366260/1‑236 (MQ=255)
                                           cgtcaCCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGt                                                                                                                                                                                                       <  2:762043/224‑1 (MQ=17)
                                               aCCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTAGGTACTGCCACCGTGACTGCTCCGGACAACg                                                                                                                                                                                                              >  2:1065590/1‑217 (MQ=255)
                                                   ccaccacCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCACCAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCg                                                                                                                                                                                                                     <  1:289212/205‑1 (MQ=255)
                                                   acacCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGGCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCt                                                                                                                                                                                                                        <  1:289560/203‑1 (MQ=38)
                                                       caccGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCCCCGGACaa                                                                                                                                                                                                                >  2:1400587/1‑207 (MQ=37)
                                                         ccGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGt                                                                                                                                                                                                             >  1:551406/1‑208 (MQ=255)
                                                                agtagtGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGa                                                                                                                                                                          <  1:686368/236‑1 (MQ=40)
                                                                 gtagtGGCCATGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCTTGACCTTCTCGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGa                                                                                                                                                                          <  2:806402/234‑1 (MQ=255)
                                                                      ggCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTg                                                                                                                                                                                                                          >  2:627791/1‑182 (MQ=37)
                                                                                    cgACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGaa                                                                                                                                                            >  2:1422154/1‑230 (MQ=40)
                                                                                     caCCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCg                                                                                                                                                      >  1:11100/2‑235 (MQ=37)
                                                                                       ccccGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTAc                                                                                                                                                                   <  2:523097/220‑1 (MQ=255)
                                                                                       ccccGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCg                                                                                                                                                      <  2:552616/233‑1 (MQ=255)
                                                                                          cGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGTCACCGTGACTGCTcc                                                                                                                                                                                                                      <  2:1419982/166‑1 (MQ=39)
                                                                                          cGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTcc                                                                                                                                                                                                                      >  1:1419985/1‑166 (MQ=39)
                                                                                              ggTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGAGGTGCGGACGTCggta                                                                                                                                            <  1:317378/236‑3 (MQ=37)
                                                                                                     gAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCTGACAACGTCTACGACGGTATCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTc                                                                                                                                       <  2:802384/234‑1 (MQ=255)
                                                                                                                  ttACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACATCCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGCCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACg                                                                                                                                                                                                        <  1:551128/156‑1 (MQ=255)
                                                                                                                          acCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAg                                                                                                                  <  2:895309/234‑1 (MQ=255)
                                                                                                                                 cGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTCCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGt                                                                                                                                                                >  1:518177/1‑181 (MQ=255)
                                                                                                                                  gCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCATGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGGACCAACGACCCTGTCGt                                                                                                                                                                                  >  2:341222/1‑162 (MQ=255)
                                                                                                                                  gCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCg                                                                                                                                                                                                                     <  2:322514/127‑1 (MQ=38)
                                                                                                                                  gCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCaa                                                                                                                                                                                              <  2:889149/150‑1 (MQ=40)
                                                                                                                                  gCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGCCCCTGTCGt                                                                                                                                                                                  >  2:397799/1‑162 (MQ=38)
                                                                                                                                  gCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGt                                                                                                          >  2:1345661/1‑234 (MQ=255)
                                                                                                                                  gCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCa                                                                                                        >  1:227210/1‑236 (MQ=255)
                                                                                                                                  gCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCGGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTc                                                                                                                                                                            <  2:76034/168‑1 (MQ=38)
                                                                                                                                                  ggCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCg                                                                                                                                                      <  2:162644/174‑1 (MQ=255)
                                                                                                                                                   gCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGCCAACGTCTACGTCga                                                                                                                                                                                                    >  2:119949/1‑126 (MQ=37)
                                                                                                                                                   gCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACAGTGACTGCTCCGGACAACGTCTACGTCGGTAc                                                                                                                                                                                                 >  2:1518988/1‑130 (MQ=255)
                                                                                                                                                   gCCTGCCGATCAACAACCACAGGGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCg                                                                                                                                                      <  2:390259/173‑1 (MQ=255)
                                                                                                                                                     cTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGAcccc                                                                                                             >  1:975067/1‑212 (MQ=255)
                                                                                                                                                       gCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGc                                                                                   >  1:438104/1‑236 (MQ=255)
                                                                                                                                                        ccGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTg                                                                                                                                                        >  1:48103/1‑166 (MQ=255)
                                                                                                                                                        ccGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGtt                                                                                                                                        >  1:886697/1‑182 (MQ=40)
                                                                                                                                                         cGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGc                                                                                   <  2:1207505/234‑1 (MQ=255)
                                                                                                                                                                          gCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAg                                                                                                                                          >  2:1112144/1‑162 (MQ=255)
                                                                                                                                                                           cATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACCCGACAAGACCCCGGTCAGCACTTCGGTCGCcga                                                                                        >  2:290408/1‑211 (MQ=255)
                                                                                                                                                                             ttGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGGCAGc                                                                                                      >  2:298140/1‑195 (MQ=255)
                                                                                                                                                                              tGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGt                                                                                                                                                                                                                                          <  2:1395528/62‑1 (MQ=31)
                                                                                                                                                                              tGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTATCAACGACCCTATCGTGATGTCGATCGCTACCGTTGAAGGTGCGGTCGTTGGCAAGTTCGAGCGACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGc                                                            <  1:1436582/236‑1 (MQ=255)
                                                                                                                                                                                aCCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCg                                                           >  1:1452410/1‑235 (MQ=255)
                                                                                                                                                                                aCCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACCCGCCAAGACCCCGGTCAGCACTTCGGTCACCGACGCGcc                                                                                  >  2:292952/1‑212 (MQ=255)
                                                                                                                                                                                aCCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGTGCGGACAGCGGCAAGTCCGAGCAACTGACACTCGACAAGACCCCGGt                                                                                                          >  2:835722/1‑188 (MQ=255)
                                                                                                                                                                                 ccTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCAGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACCTCGGTCACCGACGAGCCGGGTACCCCAGCCAACGAAGGCAAc                                                         >  1:1101117/1‑236 (MQ=255)
                                                                                                                                                                                   ttCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCa                                                                                                                                            <  1:330126/151‑1 (MQ=255)
                                                                                                                                                                                   ttCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCa                                                                                                                                            >  2:330124/1‑151 (MQ=255)
                                                                                                                                                                                   ttCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGc                                                                                                      >  1:652819/1‑189 (MQ=255)
                                                                                                                                                                                      aCGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTg                                                                                                                                                                                      >  1:153493/1‑106 (MQ=40)
                                                                                                                                                                                      aCGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGaa                                                               >  2:547229/1‑225 (MQ=2)
                                                                                                                                                                                       cGCTTAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTc                                                   <  1:869287/236‑1 (MQ=255)
                                                                                                                                                                                        gCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCg                                                                                                                                                                                                                     <  2:319184/73‑1 (MQ=255)
                                                                                                                                                                                        gCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACa                                                                                                                    >  2:829307/1‑170 (MQ=255)
                                                                                                                                                                                        gCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAgg                                                                      <  2:797391/216‑1 (MQ=255)
                                                                                                                                                                                        gCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAgg                                                                      >  1:797393/1‑216 (MQ=255)
                                                                                                                                                                                        gCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCACCGACgg                                                             >  2:448481/1‑225 (MQ=255)
                                                                                                                                                                                        gCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCCGTACCAACGACCCTGGCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCCAGTTc                                                                                                                                       >  2:365941/1‑151 (MQ=255)
                                                                                                                                                                                             gCGACGGCAAAACCGTCATCACCGTGCCGGCCAACTGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGCGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTc                                                                                                                                       >  2:202137/1‑146 (MQ=255)
                                                                                                                                                                                             gCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTGCCAACGACCCTGTCGTGATGTCGATCGCAAcc                                                                                                                                                                  >  1:1151556/1‑119 (MQ=255)
                                                                                                                                                                                             gCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAgg                                                                      >  2:604097/1‑211 (MQ=255)
                                                                                                                                                                                               gACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGt                                                    >  2:1340663/1‑227 (MQ=255)
                                                                                                                                                                                               gACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTc                                             >  2:1068816/1‑234 (MQ=255)
                                                                                                                                                                                                aCGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCAc                                           <  1:1070510/235‑1 (MQ=255)
                                                                                                                                                                                                   gCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCAcc                                                                                          <  1:667629/185‑1 (MQ=255)
                                                                                                                                                                                                     aaaaCCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGt                                                                                                                                                 >  2:855042/1‑128 (MQ=255)
                                                                                                                                                                                                     aaaaCCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCa                                      <  1:756615/235‑1 (MQ=255)
                                                                                                                                                                                                        aCCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCAcc                                                                                          <  2:1276153/180‑1 (MQ=255)
                                                                                                                                                                                                         ccGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCg                                                                                                                                                                                                                     <  1:1350984/56‑1 (MQ=255)
                                                                                                                                                                                                         ccGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACtt                                                                                                  >  1:250380/1‑171 (MQ=255)
                                                                                                                                                                                                             catcaCCGTGCCGTCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCg                                                           <  2:257333/206‑1 (MQ=255)
                                                                                                                                                                                                                   cGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACGCTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCt                       <  1:769432/236‑1 (MQ=255)
                                                                                                                                                                                                                        cGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTg                  >  1:1401857/1‑236 (MQ=255)
                                                                                                                                                                                                                         ggCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACGTCGGTCACCGCCGAGCCGGGTACCCCAGGCAACg                                                                 >  2:596432/1‑188 (MQ=255)
                                                                                                                                                                                                                         ggCCAACGGCACCGTGGGCACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTg                                                                                                                                                              >  1:176855/1‑95 (MQ=255)
                                                                                                                                                                                                                          gCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGc                <  1:1340666/236‑1 (MQ=255)
                                                                                                                                                                                                                                     cGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTc      <  1:1396760/235‑1 (MQ=255)
                                                                                                                                                                                                                                     cGTGGGCACTGCCACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACGCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTc      <  1:1165984/235‑1 (MQ=255)
                                                                                                                                                                                                                                     cGTGGGCACTGCCACCGTCACTGCCCCCGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTc      <  1:1504077/235‑1 (MQ=255)
                                                                                                                                                                                                                                          gTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCggg                                                                               >  2:118117/1‑157 (MQ=255)
                                                                                                                                                                                                                                          gTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAAc  <  2:1032515/234‑1 (MQ=255)
                                                                                                                                                                                                                                          gTACTGCCACCGTCACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTg                                                                                                                                                                                 <  1:1353662/59‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                        
CGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAAC  >  AE015451/213504‑213971

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: