Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 5,988,763 C→G 12.1% intergenic (+57/‑489) kefB‑III → / → PP_mr64 glutathione‑regulated potassium/H+ antiporter/C4 ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154515,988,7630CG12.1% 272.2 / 22.5 82intergenic (+57/‑489)kefB‑III/PP_mr64glutathione‑regulated potassium/H+ antiporter/C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base C (18/54);  new base G (6/4);  total (24/58)
Fisher's exact test for biased strand distribution p-value = 5.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CCTGAGCGATGAACAGGCGGCGGACCGCATTGCACGGTTTACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCC  >  AE015451/5988541‑5988978
                                                                                                                                                                                                                              |                                                                                                                                                                                                                       
ccTGAGCGGTGAACAGGCGGCGGACCGCATTGCACGGTTTACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGC                                                                                                                                                                                                                    <  2:365944‑M1/228‑1 (MQ=255)
ccTGAGCGATGAACAGGCGGCGGACCGCATTGCACGGTTTACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTcgcg                                                                                                                                                                                                            <  1:570220/236‑1 (MQ=255)
ccTGAGCGATGAACAGGCGGCGGACCGCATTGCACGGTTTACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTcgcg                                                                                                                                                                                                            >  1:1142419/1‑236 (MQ=255)
      cGATGAACAGGCGGCGGACCGCATTGCACGGTTTACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCGCGGGCGGCGCTcgatctcatgg                                                                                                                                                                                                      >  1:552506‑M1/1‑225 (MQ=255)
       gATGAACAGGCGGCGGACCGCATTGCACGGTTTACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTaa                                                                                                                                                                                                       >  2:537437/1‑234 (MQ=255)
        aTGAACAGGCGGCGGACCGCATTGCACGGTTTACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAAcc                                                                                                                                                                                                    <  1:1128083/236‑1 (MQ=255)
           accagGCGGCGGACCGCATTGCACGGTTTACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTcgatctcatgggcgc                                                                                                                                                                                                  <  1:447528‑M1/233‑16 (MQ=255)
           aaCAGGCGGCGGACCGCATTGCACGGTTTACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGc                                                                                                                                                                                                 <  1:346003/236‑1 (MQ=255)
            aCAGGCGGCGGACCGCATTGCACGGTTTACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCg                                                                                                                                                                                                  >  1:637562/1‑234 (MQ=255)
            aCAGGCGGCGGACCGCATTGCACGGTTTACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCt                                                                                                                                                                                                >  1:466371/1‑236 (MQ=255)
                gcggcgGACCGCATTGCACGGTTTACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGc                                                                                                                                                                                                 <  1:704689/231‑1 (MQ=255)
                gcggcgGACCGCATTGCACGGTTTACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTa                                                                                                                                                                                            <  1:731094/236‑1 (MQ=255)
                          gCATTGCACGGTTTACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTc                                                                                                                                                                                  <  1:537442/236‑1 (MQ=255)
                          gCATTGCACGGTTTACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTc                                                                                                                                                                                  <  1:450127/236‑1 (MQ=255)
                           cATTGCACGGTTTACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTCGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTc                                                                                                                                                                                  <  1:863805/235‑1 (MQ=255)
                                 aCGGTTTACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAg                                                                                                                                                                           <  1:379446/236‑1 (MQ=255)
                                   ggTTTACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGGGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCG                                                                                                                                                                                                                   >  2:10937‑M1/1‑194 (MQ=255)
                                                tGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTcgatctcatgggcgcccaaaaactatcgccagg                                                                                                                                                                                >  2:1272597‑M1/1‑183 (MQ=255)
                                                  aCGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGc                                                                                                                                                                              >  1:964812/1‑216 (MQ=255)
                                                  aCGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTc                                                                                                                                                                     <  2:637557/225‑1 (MQ=255)
                                                     aaGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCa                                                                                                                                                       <  1:147806/236‑1 (MQ=255)
                                                     aaGAGGTGCTGGAGGCCCAGGGGCAAGGGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTc                                                                                                                                                                                  <  2:1059440/209‑1 (MQ=255)
                                                      agagGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCa                                                                                                                                                       <  1:951653/235‑1 (MQ=255)
                                                        agGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCAcgc                                                                                                                                                    <  1:1020167/236‑1 (MQ=255)
                                                          gTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGAc                                                                                                                                                            <  1:598991/226‑1 (MQ=255)
                                                                   gCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTgg                                                                                                                                         <  1:1228973/236‑1 (MQ=255)
                                                                   gCCCAGGGGCAAGGGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTg                                                                                                                                          <  1:172636/235‑1 (MQ=255)
                                                                    cccAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGgc                                                                                                                                        <  1:1089519/236‑1 (MQ=255)
                                                                     ccAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTg                                                                                                                                          <  2:1142414/233‑1 (MQ=255)
                                                                     ccAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGgc                                                                                                                                        >  1:828015/1‑235 (MQ=255)
                                                                         gggCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTTCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACg                                                                                                                                   <  1:705037/236‑1 (MQ=255)
                                                                         gggCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACg                                                                                                                                   <  1:1109726/236‑1 (MQ=255)
                                                                          ggCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGgc                                                                                                                                        <  1:1131720/230‑1 (MQ=255)
                                                                              aGTGCGTGATGACCGGGCCAAGGTGAGGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGgc                                                                                                                                        <  1:118597/226‑1 (MQ=255)
                                                                                       tgaCCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCGCGGGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaatttcctgacgatcacctcgccgccatgacgcgaact                                                                                                                       >  2:911461‑M1/1‑144 (MQ=255)
                                                                                       tgaCCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctg                                                                                                                                                                          >  2:385430‑M1/1‑144 (MQ=255)
                                                                                            gggACAAGGTGATGCAGACGGCCGAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAg                                                                                                                                                                           <  2:290002/177‑1 (MQ=255)
                                                                                                 aaGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTc                                                                                                                                                                                  <  2:362545/165‑1 (MQ=255)
                                                                                                       aTGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTGGTTGGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACgtttgttt                                                                                                                                                                                   <  2:293310/158‑1 (MQ=255)
                                                                                                          cAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTcg                                                                                                  <  1:1125708/236‑1 (MQ=255)
                                                                                                            gACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACg                                                                                                                                   <  2:1121100/201‑1 (MQ=255)
                                                                                                                 ccAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGACTGTGAGATTGGGCGCCGCGCGGGCGGCGCTcgatctcatgggcgccca                                                                                                                                                                                               >  1:710015‑M1/1‑118 (MQ=255)
                                                                                                                   aaGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGcc                                                                                          <  1:694326/235‑1 (MQ=255)
                                                                                                                    aggaggCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCt                                                                                         <  1:651942/235‑1 (MQ=255)
                                                                                                                      gaggCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCGCGGGCGGCGCTcgatctcatgggcgcccaaaaactatc                                                                                                                                                                                      <  2:251683‑M1/140‑28 (MQ=255)
                                                                                                                          cgcgGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTTTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGatc                                                                                                                                                <  2:1278767/174‑3 (MQ=255)
                                                                                                                            cgGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGgc                                                                                                                                        <  1:1244849/180‑1 (MQ=255)
                                                                                                                                 ggagctggagcGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACaggagg                                                                              <  2:828010/233‑1 (MQ=255)
                                                                                                                                  gagctggagcGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACaggag                                                                               <  2:466366/231‑1 (MQ=255)
                                                                                                                                  gagctggagcGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGAtt                                                                          <  1:387225/236‑1 (MQ=255)
                                                                                                                                   agctggagcGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTGGGCGCGTGCCTGGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTcgcgcgc                                                                                             <  2:1063046/216‑1 (MQ=255)
                                                                                                                                            ggTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTcgcgcgc                                                                                             <  2:660226/207‑1 (MQ=255)
                                                                                                                                            ggTTGTTTGATTCGGACGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGGGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCt                                                                                         <  1:691782/211‑1 (MQ=255)
                                                                                                                                                gttTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaatttcctgacgatcacctcgccgccatgacgcgaactcacaaaaaccatcacaaatcgtctaattctcagctgggtcctacttccccg                                                                    <  2:612844‑M1/228‑142 (MQ=255)
                                                                                                                                                gttTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaatttcctgacgatcacctcgccgccatgacgcgaactcacaaaaaccatcacaaatcgtctaattctcagctgggtcctacttccccg                                                                    >  1:612849‑M1/1‑87 (MQ=255)
                                                                                                                                                     agtCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCCCg                                                                                                                              <  2:694980/163‑1 (MQ=255)
                                                                                                                                                         ggAGGCGGACTGAGGCCTTTTGGCGCGTGCCTGGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACg                                                                                                                                   <  2:410031/156‑1 (MQ=255)
                                                                                                                                                              gggACGGAGGCCTTTTGGGGCGTGCCTTGGGGGTTGTAGTGCCTGGGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCAcgcg                                                                                                                                                   <  2:1016844/134‑1 (MQ=255)
                                                                                                                                                              cGGACTTAGGCCTTTTGGGCCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaatttcctgacgatcacctcgccgccatgacgcgaactcacaaaaaccatcacaaatcgtctaattctcagctgggtcctacttccccgcgaaccatttcctacgttcaca                                              <  1:598142‑M1/236‑164 (MQ=255)
                                                                                                                                                                         cTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCGCGGGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaatttcctgacgatcacctcgccgccatgacgcgaactcacaaaaaccatcacaaatcgtctaattctcagctgggtcctacttccccgcgaaccatttcctacgttcacaccaatcccaag                                   <  1:911466‑M1/236‑175 (MQ=255)
                                                                                                                                                                         cTTTTGGCGCGGGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaatttcctgacgatcacctcgccgccatgacgcgaactcacaaaaaccatcacaaatcgtctaattctcagctgggtcctacttccccgcgaaccatttcctacgttcacaccaatcccaag                                   <  1:385435‑M1/236‑175 (MQ=255)
                                                                                                                                                                            ttGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaatttcctgacgatcacctcgccgccatgacgcgaactcacaaaaaccatcacaaatcgtctaattctcagctgggtcctacttccccgcgaaccatttcctacgttcacaccaatcccaaggaa                                <  1:1272602‑M1/236‑178 (MQ=255)
                                                                                                                                                                                       cTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTg                                                                                                                                          <  1:383771/119‑1 (MQ=255)
                                                                                                                                                                                       cTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTg                                                                                                                                          >  2:383766/1‑119 (MQ=255)
                                                                                                                                                                                           aGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATc                                                                                                                         <  2:277055/132‑1 (MQ=255)
                                                                                                                                                                                                   aGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaatttcctgacgatcacctcgccgccatgacgcgaactcacaaaaaccatcacaaatc                                                                                                   >  2:598926‑M1/1‑36 (MQ=255)
                                                                                                                                                                                                    gTGCCTGTGAGATTGGGCGCCGCGCGGGCGGCGCTcgatctcatgggcgccccaaaactatcgccaggcacctggaagccatcctgcaatttcctgacgatcacctcgccgccatgacgcgaactcacaaaaa                                                                                                               >  2:1071067‑M1/1‑35 (MQ=255)
                                                                                                                                                                                                       ccTGGGAGATCGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAACCCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGCCCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTcgcgcgc                                                                                             <  1:283705/148‑1 (MQ=17)
                                                                                                                                                                                                          gtgAGATTGGGCGCCGCGCGGGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaatttcctgacgatcacctcgccgccatgacgcgaactcacaaaaaccatcacaaatcgtctaattctcagctgggtcctacttccccgcgaaccatttcctacgttcacaccaatcccaaggaattatcctacacaccgtgtcgacgcccgtc   >  2:770999‑M1/1‑29 (MQ=255)
                                                                                                                                                                                                          gtgAGATTGGGCGCCGCGCGGGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaatttcctgacgatcacctcgccgccatgacgcgaactcacaaaaaccatcacaaatcgtctaattctcagctgggtcctacttcccc                                                                     >  1:25993‑M1/1‑29 (MQ=255)
                                                                                                                                                                                                          gtgAGATTGGGCGCCGCGCGGGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaatttcctgacgatcacctcgccgccatgacgcgaactcacaaaaaccatcacaaatcgtctaattctcagctgggtcctacttcccc                                                                     <  2:25992‑M1/169‑141 (MQ=255)
                                                                                                                                                                                                          gtgAGATTGGGCGCCGCCCGCGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaatttcctgacgatcacctcgccgccatgacgcgaactcacaaaaaccatcacaaatcgtctaattctcagctgggtcctacttccccgcgaaccatttcctacgttcacaccaatcccaaggaattatcctacacaccgtgtc             >  1:730345‑M1/1‑29 (MQ=255)
                                                                                                                                                                                                          gtgAGATTGGGCGCCGCCCGCGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaatttcctgacgatcacctcgccgccatgacgcgaactcacaaaaaccatcacaaatcgtctaattctcagctgggtcctacttccccgcgaaccatttcctacgttcacaccaatcccaag                                   <  1:200076‑M1/203‑175 (MQ=255)
                                                                                                                                                                                                          gtgAGATTGGGCGCCGCCCGCGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaatttcctgacgatcacctcgccgccatgacgcgaactcacaaaaaccatcacaaatcgtctaattctcagctgggtcctacttccccgcgaaccatttcctacgttcacaccaatcccaag                                   >  2:200073‑M1/1‑29 (MQ=255)
                                                                                                                                                                                                          gtgAGATTGGGCGCCGCCCGCGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaattacctgacgatcacctcgccgccatgacgcgaactcacaaaaaccatcacaaatcgtctaattctcagctgggtcctacttccccgcgaaccatttcctacgttcacaccaatcccaaggaattatcctacacaccgtgtcgacgcccgtc   >  1:416895‑M1/1‑29 (MQ=255)
                                                                                                                                                                                                           tgAGATTGGGCGCCGCCCGCGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaatttcctgacgatcacctcgccgccatgacgcgaactcacaaaaaccatcacaaatcgtctaattctcag                                                                                      >  2:572948‑M1/1‑28 (MQ=255)
                                                                                                                                                                                                           tgAGATTGGGCGCCGCCCGCGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaatttcctgacgatcacctcgccgccatgacgcgaactcacaaaaaccatcacaaatcgtctaattctcag                                                                                      <  1:572953‑M1/151‑124 (MQ=255)
                                                                                                                                                                                                              gaTTGGGCGCCGCCCGCGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaatttcctgacgatcacctcgccgccatgacgcgaactcacaaaaaccatcacaaatcgtctaattctcagctgggtcctacttccccgcgaaccatttcctacgttcaca                                              <  2:1185815‑M1/188‑164 (MQ=255)
                                                                                                                                                                                                              gaTTGGGCGCCGCCCGCGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaatttcctgacgatcacctcgccgccatgacgcgaactcacaaaaaccatcacaaatcgtctaattctcagctgggtcctacttccccgcgaaccatttcctacgttcaca                                              >  1:1185820‑M1/1‑25 (MQ=255)
                                                                                                                                                                                                               aTTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCa                                                                     <  2:964807/164‑1 (MQ=255)
                                                                                                                                                                                                                     cgccgcCCGCGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaatttcctgacgatcacctcgccgccatgacgcgaactcacaaaaaccatcacaaatcgtctaattctcagctgggtcctacttccccgcgaaccatttcctacgttcacaccaatcccaaggaattatcctacacaccgtgtcgacgcccgtct  >  1:851239‑M1/1‑18 (MQ=255)
                                                                                                                                                                                                                     cgccgcCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGaaacaa                                                              <  2:52135/164‑1 (MQ=11)
                                                                                                                                                                                                                            cgcgcgGCGCTcgatctcatgggcgcccaaaaactatcgccagg                                                                                                                                                                                <  2:821487‑M1/44‑34 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                       
CCTGAGCGATGAACAGGCGGCGGACCGCATTGCACGGTTTACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCC  >  AE015451/5988541‑5988978

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