Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 215,377 C→T 6.0% A6961A (GCC→GCT PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451215,3770CT 6.0% 229.2 / 4.0 81A6961A (GCC→GCTPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (33/43);  new base T (3/2);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CAAGTTCGAACAGCTTACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCACTTACACCATC  >  AE015451/215152‑215611
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                          
cAAGTTCGAACAGCTTACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCa                                                                                                                                                                                                                                  <  1:20140/236‑1 (MQ=31)
                          aGACGCCGGTCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCag                                                                                                                                                                                                        >  1:371907/1‑235 (MQ=17)
                               ccGGTCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCa                                                                                                                                                                                                                                  <  1:800795/205‑1 (MQ=31)
                               ccGGTCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGc                                                                                                                                                                                                     >  2:458130/1‑234 (MQ=34)
                                cGGTCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGAt                                                                                                                                                                                                  <  1:970891/236‑1 (MQ=37)
                                   tCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGa                                                                                                                                                                                                                                       >  2:992960/1‑196 (MQ=31)
                                         cGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTa                                                                                                                                                                                          <  2:889272/235‑1 (MQ=39)
                                                    gacgaGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCt                                                                                                                                                                                <  2:371904/234‑1 (MQ=9)
                                                    gacgaGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTgcgc                                                                                                                                                                                                          >  2:846425/1‑208 (MQ=32)
                                                     acgaGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCa                                                                                                                                                                                                                                  <  2:38962/183‑1 (MQ=31)
                                                              gTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGGCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGc                                                                                                                                                                                 <  2:1021675/223‑1 (MQ=255)
                                                                    caggCAACGAAGGCGACTGGGTCAAGCTCCCCATCCCCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAgcgcg                                                                                                                                                                                                                             <  1:190987/171‑1 (MQ=9)
                                                                              aGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTcc                                                                                                                                                     <  1:631948/235‑1 (MQ=39)
                                                                              aGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTcc                                                                                                                                                     <  1:602847/235‑1 (MQ=37)
                                                                                              gTCACCATCGCCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAacg                                                                                                                                                                                     <  2:108192/187‑1 (MQ=37)
                                                                                                 aCCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAACCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCCAGGCCGGCGAAccc                                                                                                                                                                                                                                   <  1:819674/138‑1 (MQ=12)
                                                                                                    aTCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCtgt                                                                                                                                                >  1:353466/1‑218 (MQ=38)
                                                                                                     tCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCa                                                                                                                                                                                                                  >  2:857734/1‑151 (MQ=34)
                                                                                                         cgccgACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCtt                                                                                                                           <  2:409004/234‑1 (MQ=37)
                                                                                                         cgccgACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTc                                                                                                                          <  1:640395/235‑1 (MQ=37)
                                                                                                            cgACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAGGGCGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCCGCTGTAGac                                                                                                                                            >  2:636519/1‑213 (MQ=255)
                                                                                                                cAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACAcc                                                                                                                                                                                                                    >  2:71765/1‑138 (MQ=34)
                                                                                                                  gACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAg                                                                                                                                                                          <  2:33132/178‑1 (MQ=38)
                                                                                                                   aCCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGcag                                                                                                                                                                                                             >  2:406618/1‑140 (MQ=21)
                                                                                                                   aCCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCa                                                                                                               >  1:738363/1‑236 (MQ=37)
                                                                                                                                agaACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGa                                                                                                                                                                                                                                       >  1:1006712/1‑103 (MQ=32)
                                                                                                                                agaACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCtg                                                                                                                                                 >  2:26118/1‑189 (MQ=35)
                                                                                                                                agaACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCACCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCt                                                                                                                                                                  >  2:300608/1‑172 (MQ=25)
                                                                                                                                 gaACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCgg                                                                                                 <  1:787165/236‑1 (MQ=34)
                                                                                                                                  aaCGTCGAGCCGACCTTCACCGTGCACGTCAACCCGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCg                                                                                                       <  1:164565/229‑1 (MQ=21)
                                                                                                                                   gcGTCAAGCCGACCTTCACCGCGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCATCATCAAGGCCGGCGAAACCAGCGCGCCATCCACCCACGCTGCGCAAGGCg                                                                                                                                                                                                    <  2:1006956/134‑1 (MQ=255)
                                                                                                                                     gTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAgcgcg                                                                                                                                                                                                                             >  1:553181/1‑108 (MQ=33)
                                                                                                                                     gTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAgcgcg                                                                                                                                                                                                                             <  2:553178/108‑1 (MQ=33)
                                                                                                                                        aaGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCCAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTg                                                                                                           <  2:685131/219‑1 (MQ=21)
                                                                                                                                            cGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGtt                                                                                         <  2:934142/233‑1 (MQ=17)
                                                                                                                                              aCCGTCACCGTGCACGGCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATc                                                                                                                                                                       <  1:745895/153‑1 (MQ=25)
                                                                                                                                                  tCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGACACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGCCTACAacga                                                                                                                                                                                    >  2:1440/1‑136 (MQ=21)
                                                                                                                                                  tCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCCGCTGTAGac                                                                                                                                            >  2:440395/1‑175 (MQ=21)
                                                                                                                                                    gccGCGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGc                                                                                                                                                                                                     <  2:161294/116‑1 (MQ=25)
                                                                                                                                                    aCCGTGCACGGCAACCGGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCa                                                                                                                                                                                                                                  <  2:338554/88‑1 (MQ=255)
                                                                                                                                                       gTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATac                                                                           <  1:220119/236‑1 (MQ=32)
                                                                                                                                                          cACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGc                                                                                                                                                                    <  1:234453/144‑1 (MQ=34)
                                                                                                                                                           aCGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGAt                                                                                                                                                                                                  >  1:839236/1‑113 (MQ=255)
                                                                                                                                                           aCGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGAt                                                                                                                                                                                                  <  2:839233/113‑1 (MQ=18)
                                                                                                                                                           aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATc                                                                                                                                                           <  2:1091055/152‑1 (MQ=33)
                                                                                                                                                           aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCtt                                                                                                                           >  2:798333/1‑184 (MQ=33)
                                                                                                                                                           aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATaccac                                                                        >  1:925960/1‑235 (MQ=33)
                                                                                                                                                           aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGGCTACAACGACGCTGGc                                                                                                                                                                             >  2:550682/1‑134 (MQ=21)
                                                                                                                                                            cGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATaccacc                                                                       >  2:538566/1‑235 (MQ=33)
                                                                                                                                                             gTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACcg                                                                      <  1:837331/235‑1 (MQ=32)
                                                                                                                                                                aaCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCtt                                                                                                                           >  1:477796/1‑179 (MQ=34)
                                                                                                                                                                aaCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCtt                                                                                                                           <  2:477793/179‑1 (MQ=34)
                                                                                                                                                                aaCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGT‑‑‑GCTACCGTTCAAGTGACCGATACCACCGACGAagt                                                               >  1:567272/1‑236 (MQ=18)
                                                                                                                                                                aaCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGACCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTAc                                                                                             <  2:748295/209‑1 (MQ=17)
                                                                                                                                                                          tGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGcc                                                        <  1:84724/236‑1 (MQ=32)
                                                                                                                                                                            gCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTgg                                                          >  2:939982/1‑232 (MQ=33)
                                                                                                                                                                            gCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCa                                                       >  1:1071141/1‑235 (MQ=33)
                                                                                                                                                                             cccACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCaa                                                      <  1:1084116/235‑1 (MQ=33)
                                                                                                                                                                                    ccTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCg                                                                                                                                                                                                    >  1:84526/1‑86 (MQ=33)
                                                                                                                                                                                     cTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACcg                                              <  1:71767/235‑1 (MQ=32)
                                                                                                                                                                                        gTCGTGACCCTGAGCAACAACGCCCAGGTCACCCTCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGAc                                          <  1:118565/236‑1 (MQ=17)
                                                                                                                                                                                          cGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTCGACGGTCGCACCtt                                                                                                                           >  2:317959/1‑153 (MQ=18)
                                                                                                                                                                                          cGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGAcc                                         <  1:538569/235‑1 (MQ=33)
                                                                                                                                                                                           gtaagccTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGAccc                                        >  1:581210/6‑235 (MQ=18)
                                                                                                                                                                                           gTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGAccc                                        >  1:71652/1‑235 (MQ=32)
                                                                                                                                                                                           gTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCGGTGGATGTAGACGGCCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCCACCGGTCAAGCGACCGAGACCACCGACGAagga                                                              >  2:287389/1‑211 (MQ=255)
                                                                                                                                                                                            tGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCACAGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCt                                      <  1:600978/236‑1 (MQ=17)
                                                                                                                                                                                                ccTGAGCAACACCGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCg                                   <  1:898333/235‑1 (MQ=17)
                                                                                                                                                                                                     gCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTa                                                                                                                                                                                          >  2:493071/1‑79 (MQ=21)
                                                                                                                                                                                                     gCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCt                                                   >  2:380754/1‑214 (MQ=33)
                                                                                                                                                                                                      caacaaCGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAag                                                                >  1:34661/1‑200 (MQ=33)
                                                                                                                                                                                                            cGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTctg                                                                                                                                                    <  2:37471/110‑3 (MQ=2)
                                                                                                                                                                                                                 aGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGc                                                                                                                                                                                                     <  2:562923/56‑1 (MQ=35)
                                                                                                                                                                                                                      aCCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTg                                                                                                           <  1:325327/141‑1 (MQ=32)
                                                                                                                                                                                                                      aCCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTg                                                                                                           >  2:325324/1‑141 (MQ=32)
                                                                                                                                                                                                                       ccATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCAc            <  1:992963/235‑1 (MQ=31)
                                                                                                                                                                                                                         aTCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGaaa                                                                                                                      <  2:451276/127‑1 (MQ=9)
                                                                                                                                                                                                                         aTCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCtt                                                                                                                           >  2:618298/1‑122 (MQ=33)
                                                                                                                                                                                                                         aTCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCtt                                                                                                                           <  1:618301/122‑1 (MQ=32)
                                                                                                                                                                                                                            aaGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACcg                                                                      >  1:597178/1‑172 (MQ=12)
                                                                                                                                                                                                                            aaGGCTGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACcg                                                                      <  2:597175/172‑1 (MQ=14)
                                                                                                                                                                                                                                ccGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGc                                                                                                                                                                                                     <  1:572083/41‑1 (MQ=35)
                                                                                                                                                                                                                                ccGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGTCCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATc              <  1:973506/224‑1 (MQ=14)
                                                                                                                                                                                                                                ccGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGTCCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATc              >  2:973503/1‑224 (MQ=14)
                                                                                                                                                                                                                                ccGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCACTTACAcc     >  2:820569/1‑233 (MQ=31)
                                                                                                                                                                                                                                 cGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCACTTACACCATc  <  1:594673/235‑1 (MQ=14)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                          
CAAGTTCGAACAGCTTACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCACTTACACCATC  >  AE015451/215152‑215611

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