Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 212,155 T→C 6.9% G5887G (GGT→GGC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451212,1550TC 6.9% 209.5 / 3.1 70G5887G (GGT→GGCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (30/35);  new base C (2/3);  total (32/38)
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

CAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGATTAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCG  >  AE015451/211924‑212387
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        
cAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGTATc                                                                                                                                                                                                                                       >  2:576321/1‑235 (MQ=17)
                     ggTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGTGGTGGACGCCACTg                                                                                                                                                                                                                 <  1:1221173/236‑1 (MQ=18)
                     ggTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTg                                                                                                                                                                                                                 <  1:926571/236‑1 (MQ=33)
                                                  tAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCt                                                                                                                                                                                        <  2:1136997/232‑1 (MQ=33)
                                                          aaCGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAg                                                                                                                                                                                          >  2:1138201/1‑222 (MQ=32)
                                                           aCGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCt                                                                                                                                                                           <  1:576324/236‑1 (MQ=255)
                                                                     gACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGCCg                                                                                                                                                                                                                       >  1:912203/1‑182 (MQ=18)
                                                                              tGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGac                                                                                                                                                         >  1:1211140/1‑235 (MQ=255)
                                                                                ttCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCt                                                                                                                                                                                        <  2:293027/202‑1 (MQ=33)
                                                                                tgcACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGacac                                                                                                                                                       <  1:824030/233‑1 (MQ=11)
                                                                                   aCGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGAcaccac                                                                                                                                                    >  1:829176/1‑235 (MQ=255)
                                                                                   aCGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGCCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGCGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACaccacc                                                                                                                                                   >  1:1268345/1‑236 (MQ=255)
                                                                                    cGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCtt                                                                                                                                                                                                       <  1:948967/183‑1 (MQ=33)
                                                                                    cGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACcg                                                                                                                                                  <  1:1055175/236‑1 (MQ=255)
                                                                                          ccAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAag                                                                                                                                            <  1:1233573/236‑1 (MQ=255)
                                                                                                       cACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGccg                                                                                                                                                                             >  2:406123/1‑188 (MQ=11)
                                                                                                        aCGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGTGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCt                                                                                                                               <  1:688697/235‑1 (MQ=255)
                                                                                                                      aCCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAagtagt                                                                                                                                        <  1:371980/212‑1 (MQ=255)
                                                                                                                         cTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGCCGTCTACCACGATGCTGACCAGACCAGCCTGGGTAGCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACcgacg                                                                                                                                               >  1:986467/1‑202 (MQ=255)
                                                                                                                              caacaaTGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGcgaaggtg                                                                                                                                                                     >  1:626832/1‑168 (MQ=2)
                                                                                                                              caacaaTGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCg                                                                                                        >  1:405428/1‑236 (MQ=255)
                                                                                                                                    cgCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTcc                                                                                                                                                                    >  2:962935/2‑170 (MQ=2)
                                                                                                                                    cgCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTcc                                                                                                                                                                    <  1:962938/169‑1 (MQ=2)
                                                                                                                                       ccAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGt                                                                                                <  1:709112/235‑1 (MQ=255)
                                                                                                                                       ccAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTg                                                                                               <  1:637972/236‑1 (MQ=255)
                                                                                                                                       ccAGGTCACCAGCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGa                                                                                                       <  1:535764/228‑1 (MQ=255)
                                                                                                                                              aCCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCAc                                                                                        <  1:1026265/236‑1 (MQ=255)
                                                                                                                                                cATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCAc                                                                                        <  1:702265/234‑1 (MQ=255)
                                                                                                                                                   cAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGc                                                                                                                                                                            >  2:418063/1‑147 (MQ=21)
                                                                                                                                                        ccGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGc                                                                                                                                                                                                             >  2:926589/1‑109 (MQ=35)
                                                                                                                                                        ccGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCAc                                                                                        >  1:98439/1‑226 (MQ=255)
                                                                                                                                                         cGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCAtt                                                                              <  1:829510/235‑1 (MQ=255)
                                                                                                                                                            tGACACCAGCGCGCCGTACACCCTCGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACGCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCg                                                                                                        <  2:221164/206‑1 (MQ=255)
                                                                                                                                                             gAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGt                                                                                                                                                                >  2:713607/1‑149 (MQ=255)
                                                                                                                                                              aaaCCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGGGCGACCTTCGAGAACCTCGAGcg                                                                                                                                                                                        >  1:114494/1‑123 (MQ=25)
                                                                                                                                                                 cacgcgcgcCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCgg                                                                                                                                                                             <  2:250579/129‑1 (MQ=14)
                                                                                                                                                                    gcgcgcCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTg                                                                                               <  2:829173/207‑1 (MQ=255)
                                                                                                                                                                          cGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGt                                                                                                <  1:836556/200‑1 (MQ=255)
                                                                                                                                                                                   aCGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGt                                                                                                             <  1:1057673/178‑1 (MQ=255)
                                                                                                                                                                                      cTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGa                                                                       >  1:992026/1‑213 (MQ=255)
                                                                                                                                                                                       tGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGa                                                                                                       >  2:909686/1‑180 (MQ=9)
                                                                                                                                                                                       tGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGAtt                                                >  2:873159/1‑235 (MQ=255)
                                                                                                                                                                                       tGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCgg                                                                                                                                                                             <  1:1267748/110‑1 (MQ=17)
                                                                                                                                                                                         cgcAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGa                                                                                                           >  1:130415/1‑174 (MQ=12)
                                                                                                                                                                                         cgcAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCgg                                                                                                                                                                             <  2:1143860/108‑1 (MQ=40)
                                                                                                                                                                                         cgcAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCgg                                                                                                                                                                             >  1:1143863/1‑108 (MQ=40)
                                                                                                                                                                                         cgcAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGc                                                     <  1:1136644/228‑1 (MQ=255)
                                                                                                                                                                                          gcAGGGTGATGACGTCTACAGCGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTAGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGt                                                                                                <  1:159664/184‑1 (MQ=255)
                                                                                                                                                                                            aGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACAcc                                                                                 >  1:398558/1‑197 (MQ=255)
                                                                                                                                                                                              ggTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACACAGACGGTCTGCCGATTAACAACca                                        >  1:1116240/1‑236 (MQ=255)
                                                                                                                                                                                              ggTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGATTAACAAcc                                         <  1:873162/235‑1 (MQ=255)
                                                                                                                                                                                              ggTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTg                                                                                               <  1:287024/181‑1 (MQ=255)
                                                                                                                                                                                              agTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGt                                                                                                <  1:550835/179‑1 (MQ=255)
                                                                                                                                                                                                    gaCGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGATTAACAACcaca                                      >  1:1115138/1‑232 (MQ=255)
                                                                                                                                                                                                     aCGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGATTAACAACCACAGTGc                                  <  1:264975/235‑1 (MQ=255)
                                                                                                                                                                                                      cGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGAcacc                                                                                                                                                      >  1:231072/1‑118 (MQ=255)
                                                                                                                                                                                                      cGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGAcacc                                                                                                                                                      <  2:231072/118‑1 (MQ=255)
                                                                                                                                                                                                      cGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCAc                                                                                        >  1:560864/1‑180 (MQ=255)
                                                                                                                                                                                                      cGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCAc                                                                                        <  2:560861/180‑1 (MQ=255)
                                                                                                                                                                                                         cTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGATTAACAACCACAGTGCATTGa                             <  1:1053187/236‑1 (MQ=255)
                                                                                                                                                                                                          tACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGATTAACAACCACAGTGCATTGAc                            <  1:471683/236‑1 (MQ=255)
                                                                                                                                                                                                           aCAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCACGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCAt                                                                               >  1:44777/1‑184 (MQ=255)
                                                                                                                                                                                                           aCAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGATTaacaac                                          >  1:447445/1‑221 (MQ=255)
                                                                                                                                                                                                                gacgCTGGGCGGATCAGCCGGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCAtt                                                                              <  2:130415/177‑1 (MQ=255)
                                                                                                                                                                                                                    cTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCg                                                                                                                                                                                   <  1:781216/75‑1 (MQ=34)
                                                                                                                                                                                                                     tGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGc                                                                                                                                                                               >  2:1078724/1‑78 (MQ=34)
                                                                                                                                                                                                                     tGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGGGCGGCTTCGGGCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACGGCGACCCCGTCGGGGACCGAAGGCGGGgaga                                                                                            >  1:380764/1‑161 (MQ=255)
                                                                                                                                                                                                                          aGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGATTAACAACCACAGTGCATTGACCTTCACCCTGAGCGAc            >  1:632927/1‑236 (MQ=255)
                                                                                                                                                                                                                             tCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGATTAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGc         <  1:1210955/236‑1 (MQ=255)
                                                                                                                                                                                                                                    gggTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGATTAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCg  >  1:592878/1‑236 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        
CAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGATTAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCG  >  AE015451/211924‑212387

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