Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 203,674 T→C 7.6% D3060D (GAT→GAC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451203,6740TC 7.6% 258.3 / 4.6 85D3060D (GAT→GACPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (34/44);  new base C (3/4);  total (37/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACTACCGTCA  >  AE015451/203444‑203906
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        
acCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCt                                                                                                                                                                                                                                     <  1:589830/236‑1 (MQ=38)
  ccACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATGACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACACAGATGGTCTGc                                                                                                                                                                                                                                   <  1:1062234/236‑1 (MQ=255)
  ccACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGc                                                                                                                                                                                                                                   <  1:309480/236‑1 (MQ=37)
          cgcAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATcaa                                                                                                                                                                                                                            <  1:842972/235‑1 (MQ=40)
          cgcAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCAACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATc                                                                                                                                                                                                                              <  2:608345/233‑1 (MQ=17)
              gggTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAAcc                                                                                                                                                                                                                       <  1:1092638/236‑1 (MQ=40)
               ggTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAAcc                                                                                                                                                                                                                       <  1:708474/235‑1 (MQ=38)
                      aCGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGc                                                                                                                                                                                                                <  1:879198/235‑1 (MQ=38)
                       cGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTg                                                                                                                                                                                                                                    >  2:1013779/1‑214 (MQ=38)
                            atcaggacgCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATcaa                                                                                                                                                                                                                            <  1:1096218/212‑1 (MQ=255)
                                         ccAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCAACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGc                                                                                                                                                                                                                <  1:788519/216‑1 (MQ=25)
                                               cAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCg                                                                                                                                                                                                                                 <  1:559990/193‑1 (MQ=39)
                                               cAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGc                                                                                                                                                                                       <  1:1034498/235‑1 (MQ=38)
                                                aGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGa                                                                                                                                                                                                                                >  2:788518/1‑193 (MQ=25)
                                                aGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGc                                                                                                                                                                                       <  1:459707/234‑1 (MQ=40)
                                                         aTCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTg                                                                                                                                                                                                                 >  1:858725/1‑199 (MQ=40)
                                                              cTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGAt                                                                                                                                                                                                                               >  1:283502/1‑180 (MQ=39)
                                                               tCGGCGGTGGACGCCACTGGTGCAACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGt                                                                                                                                                                      >  1:447376/1‑236 (MQ=17)
                                                                           gCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAAc                                                                                                                                                            >  1:1066662/1‑234 (MQ=37)
                                                                            ccACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTg                                                                                                                                                                     <  1:995143/224‑1 (MQ=21)
                                                                                tGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATcaaca                                                                                                                                                                                                                          <  1:925208/167‑1 (MQ=255)
                                                                                tGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATcaaca                                                                                                                                                                                                                          >  2:925207/1‑167 (MQ=40)
                                                                                 ggTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATcaaca                                                                                                                                                                                                                          <  2:963631/166‑1 (MQ=255)
                                                                                   tGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCAcc                                                                                                                                                      <  2:795650/232‑1 (MQ=39)
                                                                                    gcaaCCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGt                                                                                                                                                    <  2:447375/230‑1 (MQ=37)
                                                                                     cGACCTTCGTGAACCTTGAGCTGGGCGGTGCCGCGAAGGTGGATGTCATCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCg                                                                                                                                                                                                                                 <  2:1048370/155‑1 (MQ=255)
                                                                                     cGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCAt                                                                                                                                                                                                              >  1:747637/1‑174 (MQ=39)
                                                                                     cGACCTCCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGc                                                                                                                                                                                                                                   <  2:1090841/153‑1 (MQ=21)
                                                                                     cGACCATCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCCTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGc                                                                                                                                                                                                                                   <  2:627524/153‑1 (MQ=255)
                                                                                      gACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTTCATGTCACCGACACCACCGACGAAGTCGTGGCCAAGCTGACCGCCACCCCGTCGGTCACCGAAGGTGGCGAGATCACTTACACCATCACCCTGACCAACAAAGACGGCCTGCCg                                                                                                                                                                                                                                 <  2:549202/154‑1 (MQ=1)
                                                                                          ttCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGc                                                                                                                                                                                                                                   <  2:958737/148‑1 (MQ=25)
                                                                                          ttCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGCCCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCaa                                                                                                                                                                                     <  2:1045379/194‑1 (MQ=38)
                                                                                                aaCCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCg                                                                                                                                      <  1:574277/235‑1 (MQ=25)
                                                                                                  ccTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTg                                                                                                                                    <  1:932470/235‑1 (MQ=40)
                                                                                                    gggAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGc                                                                                                                                                <  1:221361/220‑1 (MQ=37)
                                                                                                       cgCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTCCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTggg                                                                                                                                                 <  1:1110309/216‑1 (MQ=38)
                                                                                                             gCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCgg                                                                                                                          <  2:860545/234‑1 (MQ=255)
                                                                                                               ggTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGGCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGGGCATTGACCTTCACGCTGAGCGACGGCAAAACCGt                                                                                                                                                                               >  1:127605/1‑179 (MQ=255)
                                                                                                               ggTGCCGCGAAGGTGGATGTCACCGACACCAACGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGACGGACTGcc                                                                                                                                                                                                                                  >  1:938550/1‑128 (MQ=12)
                                                                                                                gTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACa                                                                                                                       >  2:675018/1‑234 (MQ=255)
                                                                                                                        aaGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACg                                                                                                              <  2:560647/235‑1 (MQ=255)
                                                                                                                         aGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTc                                                                                                            >  1:17049/1‑236 (MQ=255)
                                                                                                                          ggTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCa                                                                                                                                                              >  1:874519/1‑185 (MQ=39)
                                                                                                                          ggTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCa                                                                                                                                                              <  2:874518/185‑1 (MQ=38)
                                                                                                                           gTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCgg                                                                                                          <  1:1114439/236‑1 (MQ=255)
                                                                                                                              gATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGt                                                                                                         >  1:739066/1‑234 (MQ=255)
                                                                                                                                    aCCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGa                                                                                                 <  1:284588/236‑1 (MQ=40)
                                                                                                                                      cGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCa                                                                                                                                                              <  1:149147/173‑1 (MQ=39)
                                                                                                                                        acacCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGt                                                                                                                                                                      >  2:952852/1‑163 (MQ=38)
                                                                                                                                        acacCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGACGGTTTGCCGATCAATAACCACAGTGc                                                                                                                                                                                                                >  2:202273/1‑121 (MQ=255)
                                                                                                                                                                 gTTGACCGCGACCCCTTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCTCCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCgg                                                                                                                          >  1:401848/1‑182 (MQ=255)
                                                                                                                                                                          cACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTg                                                                    >  2:509280/1‑227 (MQ=37)
                                                                                                                                                                          cACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCg                                                            <  1:1079971/235‑1 (MQ=37)
                                                                                                                                                                          cACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCg                                                            <  1:152178/235‑1 (MQ=35)
                                                                                                                                                                          cACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCg                                                            <  2:770113/235‑1 (MQ=2)
                                                                                                                                                                              ccGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAAc                                                                                                                                                                                  <  2:1020572/113‑1 (MQ=37)
                                                                                                                                                                                gTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGt                                                       <  2:903228/234‑1 (MQ=37)
                                                                                                                                                                                  cGGTCACCGAAGGCGGCGAGATCTCTTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGCCTGCCCCGGACAACGCCTACGTCGGTACCAACGACCCTGTGATCaaac                                                                                   >  2:140185/1‑203 (MQ=255)
                                                                                                                                                                                  cGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGACCGCAACCg                                                                       >  2:481789/1‑216 (MQ=25)
                                                                                                                                                                                    gTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGATCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTaa                                                 >  1:637339/1‑236 (MQ=38)
                                                                                                                                                                                            aGGCGGCGAGATCACTTACACCATCACCCTGACCATCAAAGATGGTCTGCCTATCAACAAACACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAATGTGCGGACGTCGGTAAGTTCGAGc                                         <  1:627072/236‑1 (MQ=255)
                                                                                                                                                                                              gcggcgAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCg                                                                                                                           >  1:459852/1‑152 (MQ=255)
                                                                                                                                                                                              gcggcgAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCg                                                                                                                           <  2:459851/152‑1 (MQ=40)
                                                                                                                                                                                               cggcgAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAAc                                      <  1:640011/236‑1 (MQ=14)
                                                                                                                                                                                                  cgAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTg                                    >  1:700747/1‑235 (MQ=39)
                                                                                                                                                                                                     gaTCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCAACGTGCCGGCCaa                                                                                                                                                             >  1:317930/1‑111 (MQ=25)
                                                                                                                                                                                                              cacCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGa                                                                   >  2:495643/1‑192 (MQ=35)
                                                                                                                                                                                                               acCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCGCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGAc                      <  1:315496/236‑1 (MQ=18)
                                                                                                                                                                                                                ccATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGAc                      <  1:441839/235‑1 (MQ=37)
                                                                                                                                                                                                                 cttCACCCGGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCg                                                                                                                                      <  1:244966/120‑1 (MQ=25)
                                                                                                                                                                                                                  aTCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGcc                   <  1:906238/236‑1 (MQ=37)
                                                                                                                                                                                                                      cccTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCCTCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCtg                                                                                            >  2:911882/1‑159 (MQ=25)
                                                                                                                                                                                                                      cccTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTCCGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGGTCGAGCaa                                       >  2:280530/1‑212 (MQ=25)
                                                                                                                                                                                                                      cccTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAat                                                                                   >  2:336392/1‑168 (MQ=35)
                                                                                                                                                                                                                      cccTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTc                                             >  1:536305/1‑206 (MQ=37)
                                                                                                                                                                                                                      cccTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACg                                 >  1:969825/1‑218 (MQ=37)
                                                                                                                                                                                                                      cccTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGt                >  1:366283/1‑235 (MQ=35)
                                                                                                                                                                                                                       ccTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTcg              >  2:877176/1‑235 (MQ=37)
                                                                                                                                                                                                                        cTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGcc                   >  2:752975/1‑230 (MQ=35)
                                                                                                                                                                                                                        cTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGcc                   <  1:752976/230‑1 (MQ=35)
                                                                                                                                                                                                                            ccAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTa                                                                                                                >  1:111673/1‑133 (MQ=37)
                                                                                                                                                                                                                            ccAACAAAGATGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTa                                                                                                                <  2:111673/133‑1 (MQ=17)
                                                                                                                                                                                                                              aacaaAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGAACCAACGACCCTGTGATCAAACCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACgat                               >  1:858495/1‑210 (MQ=25)
                                                                                                                                                                                                                                  aaGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGtt                                              >  2:283318/1‑193 (MQ=37)
                                                                                                                                                                                                                                   aGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACTACCGTCa  >  1:951183/1‑236 (MQ=35)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        
ACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACTACCGTCA  >  AE015451/203444‑203906

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