Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,907,411 Δ1 bp 6.7% coding (378/939 nt) PP_4316 ← D‑isomer specific 2‑hydroxyacid dehydrogenase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,907,4110C. 6.7% 205.3 / 8.4 60coding (378/939 nt)PP_4316D‑isomer specific 2‑hydroxyacid dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base C (30/26);  new base . (2/2);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

ACGGCTTCAGCCAAGTGTCGGGGGCTGTCGTACCAGGCATAGGGCACATCCTGGCGTGCTGTGCGGCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGACCTGCTCGTAGCCGGCGCCAATCACGCAGATGATCTGCAATCTGGGTAATGCGTCAATTTCGGTGGCGGATAGCCCCAGCGGCCCACG  >  AE015451/4907208‑4907635
                                                                                                                                                                                                           |                                                                                                                                                                                                                                 
aCGGCTTCAGCCAAGTGTCGGGGGCTGTCGTACCAGGCATAGGGCACATCCTGGCGTGCTGTGCGGCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTTCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCTCGTGCTGGCATCCgc                                                                                                                                                                                                                    >  1:314614/1‑218 (MQ=255)
aCGGCTTCAGCCAAGTGTCGGGGGCTGTCGTACCAGGCATAGGGCACATCCTGGCGTGCTGTGCGGCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCgca                                                                                                                                                                                                   >  1:692270/1‑235 (MQ=255)
            aaGTGTCGGGGGCTGTCCTACCAGGCATAGGGCACATCCTGGCGTGCTCTGCGGGTGTGGTAGCTGACCTGCACGTCGAAGCCCAAGTGGGCCCGCTTGGCTATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACTAGGTTCCATTCACCGCGCCGCGTGCTGGCa                                                                                                                                                                                                                          <  2:599616/201‑1 (MQ=255)
               tgtCGGGGGCTGTCGTACCAGGCATAGGGCACATCCTGGCGTGCTGTGCGGCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCcagcagcag                                                                                                                                                                                   <  1:739223/236‑1 (MQ=255)
                 tCGGGGGCTGTCGTACCAGGCATAGGGCACATCCTGGCGTGCTGTGCGGCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCTGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGgc                                                                                                                                                                                                                    <  1:611924/201‑1 (MQ=255)
                 tCGGGGGCTGTCGTACCAGGCATAGGGCACATCCTGGCGTGCTGTGCGGCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTg                                                                                                                                                                                                                             <  2:1065505/193‑1 (MQ=255)
                 tCGGGGGCTGTCGTACCAGGCATAGGGCACATCCTGGCGTGCTGTGCGGCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGgc                                                                                                                                                                                                                    >  2:611923/1‑201 (MQ=255)
                              tACCAGGCATAGGGCACATCCTGGCGTGCTGTGCGGCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCCTTCACCGCGCCGCGTGCTGGCAt                                                                                                                                                                                                                        >  1:305843/1‑184 (MQ=255)
                              tACCAGGCATAGGGCACATCCTGGCGTGCTGTGCGGCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTgg                                                                                                                                                                                                                            <  2:492991/181‑1 (MQ=255)
                              tACCAGGCATAGGGCACATCCTGGCGTGCTGTGCGGCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGGGTGGTc                                                                                                                                                                    >  1:364357/1‑236 (MQ=255)
                                ccAGGCATAGGGCACATCCTGGCGTGCTGTGCGGCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGt                                                                                                                                                                                                                                 <  2:644495/174‑1 (MQ=255)
                                   ggCATAGGGCACATCCTGGCGTGCTGTGCGGCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCAccgcgccgc                                                                                                                                                                                                                                   >  2:752743/1‑169 (MQ=255)
                                      aTAGGGCACATCCTGGCGTGCTGTGCGGCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTTTGGTCGGCGACgg                                                                                                                                                            >  1:1059431/1‑236 (MQ=255)
                                                 ccTGGCGTGCTGTGCGGCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGCTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCgggggct                                                                                                                                                                                                           >  2:115634/1‑176 (MQ=255)
                                                       gTGCTGTGCGGCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCg                                                                                                                                           >  1:456989/1‑236 (MQ=255)
                                                          cggtgCGGCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGACCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGccc                                                                                                                                             <  1:873447/229‑1 (MQ=255)
                                                           tgtgCGGCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGgcg                                                                                                                                        <  1:1050531/235‑1 (MQ=255)
                                                            gtgCGGCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCAtt                                                                                                                                                  >  1:303899/1‑224 (MQ=255)
                                                             tgCGGCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTgg                                                                                                                                                                                                                            <  2:560571/150‑1 (MQ=255)
                                                             tgCGGCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAgg                                                                                                                                                                            >  1:725568/1‑197 (MQ=255)
                                                               cGGCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTTACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGgcg                                                                                                                                        <  1:659585/231‑1 (MQ=255)
                                                                ggCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGc                                                                                                                                                                                <  1:85530/190‑1 (MQ=255)
                                                                     tgGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACCCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCcagcag                                                                                                                                                                                      >  2:11486/1‑179 (MQ=255)
                                                                       gTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTTTGGTCGGCGACGGCCGCGGCATTGGCCCCg                                                                                                                                           <  2:1059430/220‑1 (MQ=255)
                                                                          gCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGAt                                                                                                                                                                                                           >  2:273211/1‑153 (MQ=255)
                                                                          gCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTa                                                                                                                        >  2:672908/1‑236 (MQ=255)
                                                                          gCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGATGCCAGCACGCGGCGCGGGGGATGTCGCGCAGCAGCGCcagcagcag                                                                                                                                                                                   <  2:500017/177‑1 (MQ=255)
                                                                          gCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGATGCCAGCACGCGGCGCGGGGGATGTCGCGCAGCAGCGCcagcagcag                                                                                                                                                                                   >  1:500018/1‑177 (MQ=255)
                                                                                     tGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCgcagc                                                                                                                                                                                                 >  1:260209/1‑152 (MQ=255)
                                                                                              ccAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAgcggcggcg                                                                                                    >  1:1071841/1‑236 (MQ=255)
                                                                                              ccAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAgcggcggcg                                                                                                    >  1:258451/1‑236 (MQ=255)
                                                                                               caggTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCcgcg                                                                                                                 <  2:401414/219‑1 (MQ=255)
                                                                                                   tGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGTGCAGCAGCGCcagcagcag                                                                                                                                                                                   <  2:847804/152‑1 (MQ=255)
                                                                                                          gCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGAc                                                                                         >  1:1091688/1‑235 (MQ=255)
                                                                                                                   tCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCg                                                                                                                                           <  2:811312/176‑1 (MQ=255)
                                                                                                                                     cgaggccgagAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGACCTGCTCGTAGCCGGCGCCAATCACGCa                                                              <  1:1025955/235‑1 (MQ=255)
                                                                                                                                       aggccgagAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCa                                                                                                             >  2:976854/1‑186 (MQ=255)
                                                                                                                                               aaTGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGATGCCAGCACGCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGACCTGCTCGTAGCCGGCGCCAATCACGCAgatgactgtct                                                   >  1:854496/1‑230 (MQ=255)
                                                                                                                                                aTGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTc                                                                                            <  2:494364/194‑1 (MQ=255)
                                                                                                                                                           cTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGACCTGCTCGTAGCCGGCGCCAATCACGCAGATGATCTGCAATCTGGGTaa                                         >  2:1063757/1‑234 (MQ=255)
                                                                                                                                                               ccACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGATGCCAGCACGCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGACCTGCTCGTAGCCGGCGcc                                                                       <  2:854495/200‑1 (MQ=255)
                                                                                                                                                                cACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGGCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACgg                                                                                                                           >  2:353017/1‑147 (MQ=255)
                                                                                                                                                                     aCCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACg                                                                                                                                                             >  2:276546/1‑108 (MQ=255)
                                                                                                                                                                        gAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGACCTGCTCGTAGCCGGCGCCAATCACGCAGATGATCTg                                                     <  1:526756/209‑1 (MQ=255)
                                                                                                                                                                             gCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCCGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAgcg                                                                                                          >  2:60132/1‑151 (MQ=255)
                                                                                                                                                                                gATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCa                                                                                                                                                    >  2:49778/1‑106 (MQ=255)
                                                                                                                                                                                gATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGACCTGCTCGTAGCCGGCGCCAATCACGCAGATGATCTGCAATCTGGGTAATGCGTCAATTTCGGTGGCGGATa                  >  1:1086861/1‑236 (MQ=255)
                                                                                                                                                                                  tAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGa                                                                                          >  1:1093933/1‑162 (MQ=255)
                                                                                                                                                                                  tAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGACCTGCTCGTAGCCGGCGCCAATCACGCAGATGATCTGCAATCTGGGTAATGCGTCAATTTCGGTGGCGGATAGc                >  2:684823/1‑236 (MQ=255)
                                                                                                                                                                                    aCCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGACCTGCTCGTAGCCGGCGCCAATCACGCAGATGATCTGCAATCTGGGTAATGCGTCAATTTCGGTGGCGGATAGccc              <  1:202680/236‑1 (MQ=255)
                                                                                                                                                                                       cGGTTCCATTCACCGCGCCGCGTGCTgg                                                                                                                                                                                                                            <  2:239613/28‑1 (MQ=255)
                                                                                                                                                                                           tCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGACCTGCTCGTAGCCGGCGCCAATCACGCAGATGATCTGCAATCTGGGTAATGCGTCAATTTCGGTGGCGGATAGCCCCAGCGGc       >  1:559219/1‑236 (MQ=255)
                                                                                                                                                                                            ccATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGACCTGCTCGTAGCCGGCGCCAATCACGCAGATGATCTGCAATCTGGGTAATGCGTCAATTTCGGTGGCGGATAGCCCCAGCGGc       <  1:273212/235‑1 (MQ=255)
                                                                                                                                                                                              aTTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGACCTGCTCGTAGCCGGCGCCAATCACGCAGATGATCTGCAATCTGGGTAATGCGTCAATTTCg                            <  1:619127/212‑1 (MQ=255)
                                                                                                                                                                                              aTTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGACCTGCTCGTAGCCGGCGCCAATCACGCAGATGATCTGCAATCTGGGTAATGCGTCAATTTCg                            >  2:619126/1‑212 (MQ=255)
                                                                                                                                                                                              aTTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGACCTGCTCGTAGCCGGCGCCAATCACGCAGATGATCTGCAATCTGGGTAATGCGTCAATTTCGGTGGCGGACAGCCCCAGCGGCCCa    >  1:390477/1‑236 (MQ=255)
                                                                                                                                                                                               ttCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGACCTGCTCGTAGCCGGCGCCAATCACGCAgatga                                                         <  2:451978/182‑1 (MQ=255)
                                                                                                                                                                                                tCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGACCTGCTCGTAGCCGGCGCCAATCACGCAGATGATCTGCAATCTGGGTAATGCGTCAATTTCGGTGGCGGATAGCCCCAGCGGCCCACg  <  1:459539/236‑1 (MQ=255)
                                                                                                                                                                                                tCACCGCGCCGCGTGCTGCCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGACCTGCTCGTAGCCTGCGCCAATCACGCAGATGATCTGCAATCTGGGTAATGCGTCAATTTCGGTGGCGGATAGCCCCAGCGGCCCACg  <  1:605430/236‑1 (MQ=255)
                                                                                                                                                                                                tCACCGCGCCGCGTGCGGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGACCTGCTCGt                                                                                 <  2:806256/157‑1 (MQ=255)
                                                                                                                                                                                                  aCCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGACCTGCTCGTAGCCGGCGCCAATCACGCAGATGATCTGCAATCTGGGTAATGCGTCAATTTCGGTGGCGGATAGc                >  1:1039495/1‑220 (MQ=255)
                                                                                                                                                                                                    cgcgcggcgTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGACCTGc                                                                                     <  2:381596/149‑1 (MQ=255)
                                                                                                                                                                                                           |                                                                                                                                                                                                                                 
ACGGCTTCAGCCAAGTGTCGGGGGCTGTCGTACCAGGCATAGGGCACATCCTGGCGTGCTGTGCGGCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGACCTGCTCGTAGCCGGCGCCAATCACGCAGATGATCTGCAATCTGGGTAATGCGTCAATTTCGGTGGCGGATAGCCCCAGCGGCCCACG  >  AE015451/4907208‑4907635

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