Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 197,551 C→G 11.0% T1019T (ACC→ACG PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451197,5510CG11.0% 214.7 / 14.0 80T1019T (ACC→ACGPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (45/26);  new base G (6/3);  total (51/29)
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.87e-01

AACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGT  >  AE015451/197318‑197778
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                   
aacaacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCa                                                                                                                                                                                                                                    <  1:35214/235‑1 (MQ=35)
           gACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGtt                                                                                                                                                                                                                            <  1:1199564/232‑1 (MQ=37)
           gACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGtt                                                                                                                                                                                                                            >  2:1199561/1‑232 (MQ=35)
              cGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                                                                                     <  1:618063/236‑1 (MQ=37)
                          cATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGACCACCTACACCGTGATCCTGAGCAATCCGGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGCAAAACCCAg                                                                                                                                                                                                          >  1:248547/1‑235 (MQ=255)
                                                                       cccGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCCGACCCCGGCCaa                                                                                                                                                                             >  2:1273834/1‑219 (MQ=21)
                                                                         cGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                                                                                     >  2:1322098/1‑177 (MQ=37)
                                                                         cGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                                                                                     <  1:1194332/177‑1 (MQ=37)
                                                                         cGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                                                                                     <  1:1322101/177‑1 (MQ=38)
                                                                                ggCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGtt                                                                                                                                                                                                                            >  1:1238796/1‑163 (MQ=17)
                                                                                              tCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACACTGTCCAACGGCCAAACCATTACCGTTGCAGCCGGCAAGaac                                                                                                                                                                                                             >  1:212778/1‑162 (MQ=255)
                                                                                                    aCTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGttctc                                                                                                                                                                                                                         >  1:1285698/1‑143 (MQ=38)
                                                                                                    aCTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAAt                                                                                                                                                                            >  2:1367338/1‑191 (MQ=25)
                                                                                                    aCTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGaa                                                                                                                                >  1:1021241/1‑235 (MQ=18)
                                                                                                         gTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCcac                                                                                                                          >  1:317620/1‑236 (MQ=17)
                                                                                                         gTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGATCGACTg                                                                                                                                                   >  1:595452/1‑211 (MQ=25)
                                                                                                           cGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCac                                                                                                                                                                                                                 >  1:1175979/1‑146 (MQ=37)
                                                                                                            gATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCa                                                                                                                                                                                                                                    <  2:796362/127‑1 (MQ=38)
                                                                                                             aTGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                                                                                     <  2:908283/141‑1 (MQ=37)
                                                                                                             aTGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                                                                                     >  1:908285/1‑141 (MQ=17)
                                                                                                                       cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGAcct                                                                                                                                                                       >  1:161433/1‑175 (MQ=25)
                                                                                                                       cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATGTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCg                                                                                                                                             >  2:403228/1‑203 (MQ=17)
                                                                                                                           aCCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                                                                                     >  2:1105417/1‑127 (MQ=17)
                                                                                                                            gccTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCcactgg                                                                                                                       <  1:845572/219‑4 (MQ=255)
                                                                                                                            cccTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCa                                                                                                                                                                                                      <  2:576575/141‑1 (MQ=12)
                                                                                                                                  cGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTc                                                                                                                                                                      >  2:83923/1‑167 (MQ=35)
                                                                                                                                  cGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAcc                                                                                                 >  1:972113/1‑236 (MQ=35)
                                                                                                                                              cGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGaa                                                                                                                                                                                                                         >  2:933012/1‑104 (MQ=40)
                                                                                                                                               ggTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGa                                                                                       <  2:1021238/233‑1 (MQ=14)
                                                                                                                                                 tCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCggc                                                                                                                   >  1:943081/1‑202 (MQ=9)
                                                                                                                                                     cGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCg                                                                                                                                                       <  2:346980/163‑1 (MQ=14)
                                                                                                                                                          gcggcgTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGCCCCCGGCCAATGACGTCTa                                                                                                                                                                    >  2:1064544/1‑145 (MQ=255)
                                                                                                                                                               gTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGa                                                                                                                                                                                                                          >  1:1524823/1‑86 (MQ=39)
                                                                                                                                                                tGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACCATTACCGTTGAAGcgg                                                                                                                                                                                                                     >  2:1479185/1‑88 (MQ=255)
                                                                                                                                                                           acCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGc                                                                                                                                                                                                       >  1:992897/1‑93 (MQ=11)
                                                                                                                                                                            ccGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAACCCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGCTTTCCAGACCCCGGCGAATGACGCCTACAACCACGGCTCGACTGTCAg                                                                                                                                               >  2:19877/1‑148 (MQ=14)
                                                                                                                                                                             cGTGACCCTGAGCAATCCTGCCCAGACGCCGGTACCCGTGACCCGGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGCCCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCAGCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGGCCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCg                                                      <  1:393191/236‑1 (MQ=255)
                                                                                                                                                                             cGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATCACCGTTGAAGCCg                                                                                                                                                                                                                     >  1:701070/1‑77 (MQ=255)
                                                                                                                                                                                gacccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGACTGACGTCTACaacaac                                                                                                                                                             >  2:711264/4‑130 (MQ=0)
                                                                                                                                                                                gacccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCg                                                                                  >  1:5686/4‑205 (MQ=1)
                                                                                                                                                                                gATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATc                                                                                                                                     >  1:478041/1‑154 (MQ=18)
                                                                                                                                                                                  ccttGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGc                                                                                                                                                                                                                   >  2:804588/4‑74 (MQ=255)
                                                                                                                                                                                  cccTGAGCAATCCTGTCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGGCTacaaca                                                                                                                                                               >  2:856662/2‑126 (MQ=11)
                                                                                                                                                                                  cccTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATCACCGTTGAAGCGGGCa                                                                                                                                                                                                                  >  2:354836/2‑75 (MQ=11)
                                                                                                                                                                                   ccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTAc                                                                                                                                                                                                                                <  2:341622/60‑1 (MQ=35)
                                                                                                                                                                                    cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTAc                                                                                                                                                                                                                                >  1:324829/1‑59 (MQ=35)
                                                                                                                                                                                    cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGCTCGACTGTCAGCGTCACCATCGAGAACGCCACGGGCg                                                                                                                     >  1:723913/1‑166 (MQ=25)
                                                                                                                                                                                       aGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATCACCGTTGAAGCCGGCAAAACCCAGGGCAGc                                                                                                                                                                                                    <  2:1415568/84‑1 (MQ=11)
                                                                                                                                                                                       aGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTc                                                                                                             <  2:4149/171‑1 (MQ=14)
                                                                                                                                                                                        gCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCacta                                             >  2:750862/1‑234 (MQ=2)
                                                                                                                                                                                         cAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCactac                                            >  2:1400601/1‑234 (MQ=9)
                                                                                                                                                                                          aaTCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                                                                                     >  1:1291362/1‑64 (MQ=255)
                                                                                                                                                                                             ccTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAAt                                                                                                                                                                            <  2:1346861/102‑1 (MQ=9)
                                                                                                                                                                                              cTGCCCAGACGCCGGTAACCGTGACCCGGTCCAACGGCCAAACCATTACCGTTGAAGCAGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAGCGcc                                                                                                                            >  2:824240/1‑149 (MQ=255)
                                                                                                                                                                                               tGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGt                                                                                                                   >  1:393519/1‑157 (MQ=9)
                                                                                                                                                                                               tGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAAcc                                                                                                                                                                                                                                     <  1:922097/43‑1 (MQ=38)
                                                                                                                                                                                               tGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGAcc                                     >  1:1415099/1‑235 (MQ=255)
                                                                                                                                                                                                 gccAGACGCCGGTAACCGTGACCCTGTCCAACGCCCAAACCATTACCGTTGAAGCCGGCAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAAt                                                                                                                                                                            <  2:771456/96‑1 (MQ=255)
                                                                                                                                                                                                 cccAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGc                                                                                                                                                                                                    <  2:641555/74‑1 (MQ=255)
                                                                                                                                                                                                    aGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTg                                                                                                                                                                                                                           >  2:1132632/1‑48 (MQ=38)
                                                                                                                                                                                                    aGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGa                                                                                                                                                                                                                          >  1:1482161/1‑49 (MQ=38)
                                                                                                                                                                                                    aGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGc                                                                                                                                                                                                    >  1:735053/1‑71 (MQ=39)
                                                                                                                                                                                                       cGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCaaa                                                                                                                                                                                                                >  1:319966/1‑55 (MQ=14)
                                                                                                                                                                                                         ccGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAAt                                                                                                                                                                            <  1:1042519/90‑1 (MQ=38)
                                                                                                                                                                                                         ccGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTacaaca                                                                                                                                                               <  2:927965/103‑1 (MQ=35)
                                                                                                                                                                                                         ccGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTacaaca                                                                                                                                                               >  1:927967/1‑103 (MQ=35)
                                                                                                                                                                                                           ggTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCg                                                                                                                                                                                                                     <  1:1372147/47‑1 (MQ=11)
                                                                                                                                                                                                              aaCCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGa                                                                                                                                                                              <  2:1497097/83‑1 (MQ=35)
                                                                                                                                                                                                              aaCCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGa                                                                                                                                                                              >  1:1497100/1‑83 (MQ=35)
                                                                                                                                                                                                               aCCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCaa                     <  1:1318100/235‑1 (MQ=255)
                                                                                                                                                                                                               aCCGTGACCCTGTCCAACGGCCAAACCATCACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGcc                                                                                                                            <  2:1047285/132‑1 (MQ=25)
                                                                                                                                                                                                                 cGTGACCCTGTCCAACGGCCAAGCCATCACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                                                                        <  2:771791/70‑1 (MQ=255)
                                                                                                                                                                                                                   tGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAGCGCCACAGGCGGCAATTTCGAGCAACTGACCCCGAACCCGACCCCAGcctcgac                                                                          >  1:207677/1‑172 (MQ=11)
                                                                                                                                                                                                                      accTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTg                                          <  1:35911/206‑1 (MQ=255)
                                                                                                                                                                                                                       ccTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAc                                                                                                  >  1:1022722/1‑150 (MQ=12)
                                                                                                                                                                                                                           tCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCg                                                                                        <  1:1008436/156‑1 (MQ=16)
                                                                                                                                                                                                                             cAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCa                                                                                                        >  1:369282/1‑138 (MQ=12)
                                                                                                                                                                                                                              aaCGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAACTTCGAACAACTGACCCCGAACCCGACCCCAGCttcga                                                                           >  1:1175897/1‑162 (MQ=255)
                                                                                                                                                                                                                              aaCGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGCTCGACTGTCAGCGTCACCATCGAGAGCGCCACAGGCggc                                                                                                                   >  2:947156/1‑125 (MQ=12)
                                                                                                                                                                                                                               aCGGCCAAGCCATCACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAggcggc    <  1:1446505/236‑1 (MQ=255)
                                                                                                                                                                                                                                cGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAggcggc    <  2:1415096/235‑1 (MQ=255)
                                                                                                                                                                                                                                    cagatcATCACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGt  <  2:913830/228‑1 (MQ=255)
                                                                                                                                                                                                                                        ccATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCCTCGAGAACGCCACAGGCg                                                                                                                     >  2:72424/1‑114 (MQ=17)
                                                                                                                                                                                                                                        ccATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGCAATTTCGAACAACTGACCCCGAACCCGACGCCAGCCCCGACCACCATCACCGa                                                              >  1:1358842/1‑169 (MQ=11)
                                                                                                                                                                                                                                        ccATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCCCCGACTCGGTCGAcaccacca                                             >  2:35911/1‑184 (MQ=12)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                   
AACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGT  >  AE015451/197318‑197778

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