Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 196,006 T→C 11.0% N504N (AAT→AAC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451196,0060TC11.0% 255.6 / 6.8 92N504N (AAT→AACPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (33/48);  new base C (3/8);  total (36/56)
Fisher's exact test for biased strand distribution p-value = 5.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACAC  >  AE015451/195779‑196230
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                
ggCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCa                                                                                                                                                                                                                          <  1:663874/236‑1 (MQ=38)
        cGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGACTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCgg                                                                                                                                                                                                                  >  1:836649/1‑236 (MQ=38)
           tttCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTaa                                                                                                                                                                                                               <  1:520176/236‑1 (MQ=38)
            ttCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTaa                                                                                                                                                                                                               <  1:728698/235‑1 (MQ=38)
             tacAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATcc                                                                                                                                                                                                                                <  1:303318/215‑1 (MQ=38)
                       ggCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGc                                                                                                                                                                                                                     <  1:1013457/218‑1 (MQ=38)
                                    tACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAAc                                                                                                                                                                                      <  1:1023956/236‑1 (MQ=40)
                                    tACAACAACGGCTCGACCGTCAGCGGCACCAGCGAGAGCGCCACTGGCGGCAACTGCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                                                                                                    <  1:295210/206‑1 (MQ=255)
                                                tCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGAGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                          <  1:683349/236‑1 (MQ=38)
                                                 cGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGc                                                                                                                                                                                                                     <  2:533629/192‑1 (MQ=35)
                                                  gACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                                                                                                    <  2:737000/192‑1 (MQ=17)
                                                       tCAGCGTCACCATCGAGAGCGCCACCGGCGGCAACTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCa                                                                                                                                                                                                                          <  2:1428755/181‑1 (MQ=9)
                                                          gCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                                                                                                    <  2:579904/184‑1 (MQ=25)
                                                          gCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                                                                                                    >  1:579906/1‑184 (MQ=37)
                                                                 cATCGAGACCGCCACTGGCGGCAATTTCGAGCAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACg                                                                                                                                                                                     <  2:1064229/208‑1 (MQ=255)
                                                                           gcCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACtg                                                                                                                                                                                                            >  1:1193660/1‑175 (MQ=38)
                                                                           gcCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACtg                                                                                                                                                                                                            <  2:1193657/175‑1 (MQ=38)
                                                                           gcCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGt                                                                                                                                                                                                   <  2:1038486/184‑1 (MQ=37)
                                                                           gcCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGt                                                                                                                                                                                                   >  1:1038489/1‑184 (MQ=37)
                                                                              aCTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGcc                                                                                                                                                                                                                            >  1:275330/1‑156 (MQ=33)
                                                                                tGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                     >  2:1332626/1‑209 (MQ=40)
                                                                                tGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                     <  1:1332629/209‑1 (MQ=255)
                                                                                         cTTCGAGCAGCTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCa                                                                                                                                                                                                                          <  2:676225/147‑1 (MQ=21)
                                                                                           tcaaaCAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGGCGGCGAGCCCCTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACgg                                                                                                                                                                                            <  1:999645/172‑1 (MQ=255)
                                                                                           tCGAGCACCTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                     <  2:303056/198‑1 (MQ=21)
                                                                                               aCAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGCAAAACCCAgg                                                                                                                                                         >  1:346982/1‑206 (MQ=25)
                                                                                                 aGCTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACCCCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAAc                                                                                                                                                                                      <  1:454081/175‑1 (MQ=255)
                                                                                                        cccGTACCCGACGCCGGCTCTGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGCGACCCTGAGCAATCCCGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                     <  1:836681/185‑1 (MQ=255)
                                                                                                             aaccGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                                                                                                    <  2:756904/131‑1 (MQ=21)
                                                                                                               gcGACCCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACCACCGCGGCCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                                                                                                    <  1:202140/130‑1 (MQ=12)
                                                                                                               ccGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAAccccg                                                                                                                                                          >  2:55257/1‑187 (MQ=39)
                                                                                                                cGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                     >  2:853483/1‑177 (MQ=25)
                                                                                                                cGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAAc                                                                                                                                                                                      <  1:851162/160‑1 (MQ=255)
                                                                                                                cGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                     <  1:853485/177‑1 (MQ=25)
                                                                                                                cGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCAACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAAc                                                                                                                                                                                      >  2:851160/1‑160 (MQ=255)
                                                                                                                  aCGCCGGCTCAGGCCACGAGCAGCGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTc                                                                                                                                                                                                  <  1:493503/146‑1 (MQ=255)
                                                                                                                    gCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                     <  1:478565/173‑1 (MQ=37)
                                                                                                                      cGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACCATTACCGTTGa                                                                                                                                                                          <  2:71188/166‑1 (MQ=14)
                                                                                                                               ccacagtCACCGACTCGGTCGACACCACCACCGCGACCCTGACCGCGAGCCCGTCCGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                          <  1:1157426/151‑1 (MQ=255)
                                                                                                                                 aCGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                                                                                                    >  1:183736/1‑113 (MQ=35)
                                                                                                                                         cgtcTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGACCACCTACGCCGTGACCCGGAGCAATCCTGCCTAGACGCCGGTAACTGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                     <  1:594506/149‑1 (MQ=11)
                                                                                                                                          gACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCa                                                                                                                                                                                                                          <  2:160519/98‑1 (MQ=37)
                                                                                                                                           aCTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                                                                                                    <  1:1203669/103‑1 (MQ=35)
                                                                                                                                           aCTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGccc                                                                                                                                                                                                                           <  2:755874/96‑1 (MQ=18)
                                                                                                                                           aCTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTaa                                                                                                                                                                                                               >  1:1395756/1‑108 (MQ=17)
                                                                                                                                                gTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGcc                                                                                                                                                                                                                            <  2:552862/90‑1 (MQ=39)
                                                                                                                                                 tCGACACCACCACCGCGACCCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                                                                                                    <  1:185876/97‑1 (MQ=255)
                                                                                                                                                  cGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                                                                                                    <  2:1328804/96‑1 (MQ=35)
                                                                                                                                                  cGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACCATTACCGTTGa                                                                                                                                                                          >  1:1001873/1‑138 (MQ=37)
                                                                                                                                                         gccaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGc                                                                                                                                                    <  2:261728/152‑1 (MQ=38)
                                                                                                                                                             ccGCGACCCTGACGGCGAGCCCGTCAGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCACCGTCGATTTCCAGACCCCGGCCAAGGACGTCTACAACAACGGTTCGa                                                                                                      >  2:795109/1‑195 (MQ=255)
                                                                                                                                                                                 ccGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCaa                                                                 >  1:220046/1‑212 (MQ=12)
                                                                                                                                                                                  cGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                     <  1:864905/111‑1 (MQ=255)
                                                                                                                                                                                  cGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                     >  2:864903/1‑111 (MQ=255)
                                                                                                                                                                                       gtgaCCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACtgtg                                                                                                                                                                                                          <  2:471413/66‑1 (MQ=1)
                                                                                                                                                                                       gTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACtgtg                                                                                                                                                                                                          <  1:140866/69‑1 (MQ=40)
                                                                                                                                                                                       gTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                     <  1:1117487/106‑1 (MQ=255)
                                                                                                                                                                                        tgaCCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACg                                                                              <  2:907747/190‑1 (MQ=255)
                                                                                                                                                                                        tCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTaa                                                                                                                                                                                                               <  2:410194/63‑1 (MQ=33)
                                                                                                                                                                                          tccGAAGGCGGCGTGATCACCGACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCgg                                 <  2:1409266/234‑1 (MQ=255)
                                                                                                                                                                                             gAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACtgtg                                                                                                                                                                                                          <  2:798658/63‑1 (MQ=37)
                                                                                                                                                                                                        gATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAAc                                                                                                                                                                                      <  1:952837/72‑1 (MQ=39)
                                                                                                                                                                                                        gATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAAc                                                                                                                                                                                      >  2:952835/1‑72 (MQ=39)
                                                                                                                                                                                                              cTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCa                                                                                                                              <  2:519773/122‑1 (MQ=18)
                                                                                                                                                                                                                    cGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGAc                                                                                                                         <  2:1476326/121‑1 (MQ=40)
                                                                                                                                                                                                                    cGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGAc                                                                                                                         >  1:1476329/1‑121 (MQ=40)
                                                                                                                                                                                                                    cGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGt        <  1:40758/234‑1 (MQ=255)
                                                                                                                                                                                                                     gTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAAc                                                                                                                                                                                              >  1:717842/1‑51 (MQ=35)
                                                                                                                                                                                                                       gACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAAc                                                                                                                                                                                      <  2:1387219/57‑1 (MQ=255)
                                                                                                                                                                                                                       gACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACggg                                                                       >  1:519328/1‑168 (MQ=255)
                                                                                                                                                                                                                       gACCCTGAGCAACCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAg                                                                                                                                                                        >  2:1117037/1‑71 (MQ=21)
                                                                                                                                                                                                                         cccTGAGCAATCCTGCCCAGACGCCGGTCACTGCGACACTTTTCCACGCCCAAGCGATTCCCGTTGaaccc                                                                                                                                                                      >  1:1003277/1‑68 (MQ=9)
                                                                                                                                                                                                                         cccTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAAc                                                                                                                                                                                              >  1:534309/1‑47 (MQ=38)
                                                                                                                                                                                                                         cccTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACgg                                                                                                                                                                                            >  1:1436742/1‑49 (MQ=38)
                                                                                                                                                                                                                         cccTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGATGTCTACAACGACGGtac                                                                                                        >  1:492754/1‑131 (MQ=255)
                                                                                                                                                                                                                         cccTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCg                                                                                                       >  1:405643/1‑134 (MQ=40)
                                                                                                                                                                                                                         cccTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCCGACCa                        >  1:282123/1‑213 (MQ=255)
                                                                                                                                                                                                                         cccTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGatgc  >  1:1411232/1‑232 (MQ=255)
                                                                                                                                                                                                                         cccTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGacac  >  1:1081394/1‑235 (MQ=255)
                                                                                                                                                                                                                         cccTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGacac  >  1:688937/1‑235 (MQ=255)
                                                                                                                                                                                                                         cccTGAGCAACCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAAc                                                                                                                                                                                              <  2:50818/47‑1 (MQ=255)
                                                                                                                                                                                                                         cccTGAGCAACCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAAc                                                                                                                                                                                              >  1:50818/1‑47 (MQ=255)
                                                                                                                                                                                                                          ccTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGa                                       >  1:103229/1‑197 (MQ=255)
                                                                                                                                                                                                                           cTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGt                                                                                                                                                                                                   >  1:541325/1‑40 (MQ=40)
                                                                                                                                                                                                                           cTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCa                                                                                                                              <  2:1296230/109‑1 (MQ=40)
                                                                                                                                                                                                                           cTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCa                                                                                                                              >  1:1296233/1‑109 (MQ=40)
                                                                                                                                                                                                                           cTGAGCAACCCTGCCCAGACGCCGGTAACTGTGACCCTGt                                                                                                                                                                                                   <  2:237778/40‑1 (MQ=255)
                                                                                                                                                                                                                           cTGAGCAACCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACgg                                                                                                                                                                                            <  2:976051/47‑1 (MQ=255)
                                                                                                                                                                                                                           cTGAGCAACCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                        <  2:1007244/99‑1 (MQ=18)
                                                                                                                                                                                                                           cTGAGCAACCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCa                                                        >  2:1051264/1‑179 (MQ=255)
                                                                                                                                                                                                                           cTGAGCAACCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCa                                                        <  1:1051267/179‑1 (MQ=255)
                                                                                                                                                                                                                               gCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGCAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGAcccc                       >  1:212677/1‑206 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                
GGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACAC  >  AE015451/195779‑196230

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: