Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 6,332 C→A 6.3% R17R (CGG→CGT trmE ← GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154516,3320CA 6.3% 282.4 / 9.0 79R17R (CGG→CGTtrmEGTPase
Reads supporting (aligned to +/- strand):  ref base C (35/39);  new base A (3/2);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

TCCAAGACCACAGGGCCACCGTGTCCCTGCAGTTCAAGCACGTCCTCGCCAGTGAACGAATTAGGCCCTGGAAAGAACAGAGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAA  >  AE015451/6108‑6553
                                                                                                                                                                                                                                |                                                                                                                                                                                                                             
tCCAAGACCACAGGGCCACCGTGTCCCTGCAGTTCAAGCACGTCCTCGCCAGTGAACGAATTAGGCCCTGGAAAGAACAGAGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTgg                                                                                                                                                                                                                      <  1:547769/234‑1 (MQ=255)
   aaGACCACAGGGCCACCGTGTCCCTGCAGTTCAAGCACGTCCTCGCCAGTGAACGAATTAGGCCCTGGAAAGAACAGAGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGt                                                                                                                                                                                                                 >  1:979657/1‑236 (MQ=255)
             ggCCACCGTGTCCCTGCAGTTCAAGCACGTCCTCGCCAGTGAACGAATTAGGCCCTGGAAAGAACAGAGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGc                                                                                                                                                                                                                           >  1:905306/1‑216 (MQ=255)
             ggCCACCGTGTCCCTGCAGTTCAAGCACGTCCTCGCCAGTGAACGAATTAGGCCCTGGAAAGAACAGAGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATgg                                                                                                                                                                                                         >  2:402892/1‑234 (MQ=255)
                     tgtCCCTGCAGTTCAAGCACGTCCTCGCCAGTGAACGAATTAGGCCCTGGAAAGAACAGAGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCg                                                                                                                                                                                                             >  2:458838/1‑222 (MQ=255)
                          cTGCAGTTCAAGCACGTCCTCGCCAGTGAACGAATTAGGCCCTGGAAAGAACAGAGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCAc                                                                                                                                                                                          <  1:991208/236‑1 (MQ=255)
                          cTGCAGTTCAAGCACGTCCTCGCCAGTGAACGAATTAGGCCCTGGAAAGAACAGAGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCAc                                                                                                                                                                                          <  1:586244/236‑1 (MQ=255)
                                            gtCGCCAGTGAACGAATTAGGCCCTGGAAAGAACAGAGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCg                                                                                                                                                                        <  2:732078/235‑1 (MQ=255)
                                            cTCGCCAGTGAACGAATTAGGCCCTGGAAAGAACAGAGCAATCCCTTCGTCCAGCACCAGCCCTTCGTCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCg                                                                                                                                                                        <  1:684390/236‑1 (MQ=255)
                                                       aCGAATTAGGCCCTGGAAAGAACAGAGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGtt                                                                                                                                                                                >  2:885153/1‑217 (MQ=255)
                                                       aCGAATTAGGCCCTGGAAAGAACAGAGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGtt                                                                                                                                                                                <  1:885156/217‑1 (MQ=255)
                                                        cGAATTAGGCCCTGGAAAGAACAGAGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCAc                                                                                                                                                                                          <  2:982292/206‑1 (MQ=255)
                                                                gCCCTGGAAAGAACAGAGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGgcggcg                                                                                                                                                                                                    <  2:1010856/188‑1 (MQ=255)
                                                                   cTGGAAAGAACAGAGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGcaca                                                                                                                                                                                     >  2:457818/1‑200 (MQ=255)
                                                                      gAAAGAACAGAGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGgcggcg                                                                                                                                                                                                    >  1:682188/1‑182 (MQ=255)
                                                                      gAAAGAACAGAGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGgcggcg                                                                                                                                                                                                    <  2:682185/182‑1 (MQ=255)
                                                                       aaaGAACAGAGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGa                                                                                                                                              >  1:361629/1‑235 (MQ=255)
                                                                       aaaGAACAGAGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGAt                                                                                                                                             >  1:1027999/1‑236 (MQ=255)
                                                                           aaCAGAGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACg                                                                                                                                                   <  1:458841/226‑1 (MQ=255)
                                                                                  cAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACg                                                                                                                                   <  1:402895/235‑1 (MQ=255)
                                                                                   aaTCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATg                                                                                                                                                                                                          >  2:498098/1‑163 (MQ=255)
                                                                                   aaTCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAAc                                                                                                                                    <  2:1027996/233‑1 (MQ=255)
                                                                                            gTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAggg                                                                                                                        <  1:1033678/236‑1 (MQ=255)
                                                                                             tCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAgggg                                                                                                                       >  1:688838/1‑236 (MQ=255)
                                                                                             tCCAGCACCAGCCCTTCAGCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACg                                                                                                                                                   <  1:239598/208‑1 (MQ=255)
                                                                                                    ccAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGCGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGAt                                                                                                                                                                                                  >  2:280091/1‑154 (MQ=255)
                                                                                                        cccTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGttttt                                                                                                               <  2:920234/233‑1 (MQ=255)
                                                                                                          cTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCAc                                                                                                                             >  1:898077/1‑217 (MQ=255)
                                                                                                          cTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTAtt                                                                                                           >  1:302301/1‑235 (MQ=255)
                                                                                                          cTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTAtt                                                                                                           >  1:536223/1‑235 (MQ=255)
                                                                                                             cATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTc                                                                                                          <  2:670424/233‑1 (MQ=255)
                                                                                                              aTCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCgg                                                                                                                                                                                                                            <  2:949184/118‑1 (MQ=255)
                                                                                                               tcgtcgCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAgg                                                                                                     >  1:269306/1‑236 (MQ=255)
                                                                                                                cggcgcgGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGc                                                                                                                     <  1:282672/216‑1 (MQ=255)
                                                                                                                  tcgCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCt                                                                                                   >  1:324413/1‑235 (MQ=255)
                                                                                                                   cgcgGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTc                                                                                                                                                                                            <  1:664135/145‑1 (MQ=255)
                                                                                                                        aaTGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCCGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGGTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAGAACGCCCCACTAGGGGGCggt                                                                                                                  >  1:225149/1‑212 (MQ=255)
                                                                                                                            gACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCTACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTAt                                                                                                            <  1:1024262/216‑1 (MQ=255)
                                                                                                                                            ccGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGcccc                                                                                                                               >  2:214701/1‑181 (MQ=255)
                                                                                                                                            ccGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAgggg                                                                                                                       >  2:619502/1‑189 (MQ=255)
                                                                                                                                            ccGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGGTTTTTAt                                                                                                            >  1:314395/1‑200 (MQ=255)
                                                                                                                                              gcgGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAACACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGccc                                                                                      >  2:738547/1‑220 (MQ=255)
                                                                                                                                              gcgGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCttt                                                                      >  1:872017/1‑236 (MQ=255)
                                                                                                                                                   gTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCCCGCACAGTGGTCATGCTCGAAAGCct                                                                                                                                                                 >  2:153901/1‑140 (MQ=255)
                                                                                                                                                      aGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTgg                                                                 <  2:361626/233‑1 (MQ=255)
                                                                                                                                                          tGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGtt                                                                                                                  >  1:736607/1‑180 (MQ=255)
                                                                                                                                                               ccGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGCTGGCGGCGATGGTTTCACGcac                                                                                                                                                                                      >  2:353145/1‑107 (MQ=255)
                                                                                                                                                                 ggTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCg                                                      >  2:865109/1‑233 (MQ=255)
                                                                                                                                                                 ggTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCACGGCCCAGGGCCGGGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCttttt                                                                    >  1:881115/1‑219 (MQ=255)
                                                                                                                                                                   tGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACg                                                 <  1:281195/236‑1 (MQ=255)
                                                                                                                                                                     aTCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGc                                               <  1:920886/236‑1 (MQ=255)
                                                                                                                                                                     aTCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGc                                               <  1:865112/236‑1 (MQ=255)
                                                                                                                                                                        aaCAGGCCGGCCTGGGCAGCCAGAGGGCCGGGCAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACg                                                 <  2:846773/231‑1 (MQ=255)
                                                                                                                                                                          cgggccggccTGGGCCGCCAGAGGGCCGGACAACCTGACAATTCCCCCACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCaaaa                                                                                                                                      <  2:149783/142‑1 (MQ=255)
                                                                                                                                                                             gccggccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCg                                                      <  2:824841/221‑1 (MQ=255)
                                                                                                                                                                             gccggccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGc                                               <  1:498101/228‑1 (MQ=255)
                                                                                                                                                                                  ccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCCCGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGAAAGCAAAAACGCCCAACTAGGGGGCGTTTTTTAt                                                                                                            >  1:25354/1‑162 (MQ=255)
                                                                                                                                                                                  ccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCACGGCCCAGGGCCGGGGCGATGGCGGCGATGGTTTc                                                                                                                                                                                            >  1:479539/1‑82 (MQ=255)
                                                                                                                                                                                  ccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCACGGCCCAGGGCCGGGGCGATGGCGGCGATGGTTTc                                                                                                                                                                                            <  2:479536/82‑1 (MQ=39)
                                                                                                                                                                                    ttGGCAGCCGGAGGGCCGGACAACCTGACAATTCCCACCCCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTAt                                                                                                            <  2:321561/160‑1 (MQ=255)
                                                                                                                                                                                     tGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCActg                               >  2:740002/1‑236 (MQ=255)
                                                                                                                                                                                          gCCAGAGGGCCGGCCAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTAt                                                                                                            <  2:269303/154‑1 (MQ=255)
                                                                                                                                                                                                     ggACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAgg               <  1:1010843/236‑1 (MQ=255)
                                                                                                                                                                                                       aCAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGa                         >  1:784576/1‑224 (MQ=255)
                                                                                                                                                                                                         aaCCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATAc                                                  <  2:61326/197‑1 (MQ=255)
                                                                                                                                                                                                         aaCCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATAc                                                  >  1:61328/1‑197 (MQ=255)
                                                                                                                                                                                                            cTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAgttgt        <  1:949117/236‑1 (MQ=255)
                                                                                                                                                                                                            cTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAgttgt        <  1:664863/236‑1 (MQ=255)
                                                                                                                                                                                                              gACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAgttgtt       >  1:852308/1‑235 (MQ=255)
                                                                                                                                                                                                              gACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTc      >  1:683422/1‑236 (MQ=255)
                                                                                                                                                                                                                  aTTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTcaca   <  1:1073931/235‑1 (MQ=255)
                                                                                                                                                                                                                   ttCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACaa  <  1:740005/235‑1 (MQ=255)
                                                                                                                                                                                                                    tCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTa                                                                                             <  2:891280/143‑1 (MQ=255)
                                                                                                                                                                                                                      ccACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTa                                                               >  2:107648/1‑171 (MQ=255)
                                                                                                                                                                                                                           ccGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAgttgt        <  2:979654/221‑1 (MQ=255)
                                                                                                                                                                                                                            cGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAgttgt        <  1:214703/220‑1 (MQ=255)
                                                                                                                                                                                                                            cGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTca     >  2:657785/1‑223 (MQ=255)
                                                                                                                                                                                                                             gCCACGGCCCAGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCaaaa                                                                                                                                      <  2:537703/93‑1 (MQ=255)
                                                                                                                                                                                                                             gCCACGGCCCAGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCaaaa                                                                                                                                      >  1:537706/1‑93 (MQ=255)
                                                                                                                                                                                                                               cccGGCCCTTGGCCGTGGCGATGGCGGCGACGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAgttgt        <  2:235559/217‑1 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                             
TCCAAGACCACAGGGCCACCGTGTCCCTGCAGTTCAAGCACGTCCTCGCCAGTGAACGAATTAGGCCCTGGAAAGAACAGAGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAA  >  AE015451/6108‑6553

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