Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 353,067 A→G intergenic (+34/+84) hisF CDS (Imidazole glycerol phosphate synthase subunit HisF) → / ← cbcV CDS (choline _ betaine _ carnitine ABC transporter ‑ ATP binding subunit) –/–

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM1353,0670AG99.7% 1459.8 / NA 392intergenic (+34/+84)hisF CDS (Imidazole glycerol phosphate synthase subunit HisF)/cbcV CDS (choline _ betaine _ carnitine ABC transporter ‑ ATP binding subunit)–/–
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (182/209);  total (183/209)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 394 total aligned reads displayed.

TGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGA  >  AE015451_phaM1/352934‑353200
                                                                                                                                     |                                                                                                                                      
tGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGc                                                                                                                                     <  2:7096875/137‑1 (MQ=255)
 gCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGCTCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCa                                                                                                                                    >  2:18412011/1‑137 (MQ=255)
   cAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGTCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACg                                                                                                                                  >  1:3465488/1‑137 (MQ=255)
   cAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCNTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACg                                                                                                                                  >  1:17174340/1‑137 (MQ=255)
   cAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACg                                                                                                                                  >  1:10286342/1‑137 (MQ=255)
   cAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACg                                                                                                                                  >  2:10465853/1‑137 (MQ=255)
   cAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACg                                                                                                                                  >  1:6452475/1‑137 (MQ=255)
   cAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACg                                                                                                                                  <  2:16089033/137‑1 (MQ=255)
    aGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGc                                                                                                                                 <  1:13387117/137‑1 (MQ=255)
    aGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGc                                                                                                                                 <  1:5150854/137‑1 (MQ=255)
    aGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTACCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGc                                                                                                                                 <  2:74193/137‑1 (MQ=255)
         ggTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGGACGCCCGcc                                                                                                                            >  1:6078938/1‑136 (MQ=255)
         ggTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCt                                                                                                                            >  1:15025577/1‑137 (MQ=255)
         ggTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCACt                                                                                                                            >  1:3805455/1‑137 (MQ=255)
          gTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTc                                                                                                                          <  1:7183236/137‑1 (MQ=255)
          gTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTc                                                                                                                          <  1:3102856/137‑1 (MQ=255)
          gTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTc                                                                                                                          <  1:258178/137‑1 (MQ=255)
          gTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTc                                                                                                                          <  2:8782976/137‑1 (MQ=255)
               ggcggcCAGCATCTTCCACTTTGGCGATTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCcaca                                                                                                                     <  2:8056372/137‑1 (MQ=255)
               ggcggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCcaca                                                                                                                     <  1:16133569/137‑1 (MQ=255)
               ggcggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCcaca                                                                                                                     <  2:15015697/137‑1 (MQ=255)
               ggcggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCcaca                                                                                                                     >  2:1296185/1‑137 (MQ=255)
                gcggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAg                                                                                                                    >  2:12622134/1‑137 (MQ=255)
                   gcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGtt                                                                                                                 <  1:9272202/137‑1 (MQ=255)
                    ccAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGttt                                                                                                                <  2:16422437/137‑1 (MQ=255)
                    ccAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGttt                                                                                                                <  1:3186746/137‑1 (MQ=255)
                       gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                             <  2:8657514/137‑1 (MQ=255)
                       gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                             <  1:12385842/137‑1 (MQ=255)
                       gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                             <  1:14318902/137‑1 (MQ=255)
                       gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                             <  1:7520421/137‑1 (MQ=255)
                       gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                             <  2:6586000/137‑1 (MQ=255)
                       gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                             <  2:14349550/137‑1 (MQ=255)
                       gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                             <  2:10192370/137‑1 (MQ=255)
                       gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                             <  2:11709440/137‑1 (MQ=255)
                       gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                             <  2:15785675/137‑1 (MQ=255)
                        cATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGGTg                                                                                                            >  1:7845269/1‑137 (MQ=255)
                        cATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGAACGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCGCTCCCACAGGTTTGCTg                                                                                                            >  2:18553988/1‑137 (MQ=255)
                         aTCTTCCACTTTGGCGAGTACACAGTGCCGTAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGGTCGCTGAGCCATTGACCGGGGCGGCTTTGCAGTCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  2:8919964/137‑1 (MQ=38)
                         aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  2:10716739/137‑1 (MQ=255)
                         aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  1:7368171/137‑1 (MQ=255)
                         aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  1:9869686/137‑1 (MQ=255)
                         aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  1:7052924/137‑1 (MQ=255)
                         aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  2:7769999/137‑1 (MQ=255)
                         aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  2:9322806/137‑1 (MQ=255)
                         aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  1:5040515/137‑1 (MQ=255)
                         aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  1:4074935/137‑1 (MQ=255)
                         aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  1:11049474/137‑1 (MQ=255)
                         aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  1:15617241/137‑1 (MQ=255)
                         aTCTTCCACTGTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTTGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  2:2183469/137‑1 (MQ=255)
                         aTCCTCCACTTTTGCGAGTACACAGTGCCGTAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGCCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACGGGTTTGCTga                                                                                                           <  2:16720757/137‑1 (MQ=38)
                           cTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTgaga                                                                                                         <  1:4421364/137‑1 (MQ=255)
                            ttCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGAt                                                                                                        <  1:9518203/137‑1 (MQ=255)
                             tCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  1:8593777/133‑1 (MQ=255)
                               cACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGAtctc                                                                                                     <  1:1910188/137‑1 (MQ=255)
                               cACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGAtctc                                                                                                     <  1:15552985/137‑1 (MQ=255)
                                aCTTTGGCGCGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACACCCGCTCCCCCCGGGTTTGCGAGAtctcc                                                                                                    >  1:6690239/1‑136 (MQ=38)
                                aCTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACCTGGCTTCGCGCGGGATCGTCGTTCGCTGATCCACTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGTACGCCCACCCCCCCAAAGTTGTTGAAAtctcc                                                                                                    >  1:5326206/1‑136 (MQ=255)
                                     ggCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAAGTTTGCTGAGATCTCTGTggg                                                                                               >  1:8372943/1‑137 (MQ=255)
                                        gAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGaac                                                                                            >  2:3129295/1‑135 (MQ=255)
                                              acaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGGCCGCTCCCACAAGTTTTTTGAGATCTCCGTGGGAGCGGGCGt                                                                                      >  1:1963045/1‑137 (MQ=38)
                                                   gCCGGAGGCCAAGGCCTACGTGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCACTCCCAAAGGGTTGCTGAGATCTCTGGGGGAGAGGGCGTGCCcg                                                                                 >  2:1363209/1‑137 (MQ=38)
                                                   gCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCCCAGGTTTTTTGGGGTCTCTGTGGGAGTGGGCGTGCCcg                                                                                 >  1:17683533/1‑137 (MQ=38)
                                                   gCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCcg                                                                                 >  2:8413735/1‑137 (MQ=255)
                                                   gCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCcg                                                                                 >  1:13739763/1‑137 (MQ=255)
                                                   gCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCACCCCCAAAGGTTTGCTGAGGGCTCTGTGGGAGCGGGCGTGccca                                                                                 >  1:11601714/1‑136 (MQ=38)
                                                      ggAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCACTCCCCCAAGTTTTTTGGGATCTCTTTGGGAGCGGGCGTGCCCGCGa                                                                              >  1:16414120/1‑137 (MQ=38)
                                                       gAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTAGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTGTGTGGGAGCGGGCGTGCCCGCGaa                                                                             >  2:3990783/1‑137 (MQ=255)
                                                        aGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCCCCCACAGGTTTGCTGAGGTTTCCGTGGGAGCGGGCGTTCCCGAGAAg                                                                            >  1:3609219/1‑137 (MQ=38)
                                                        aGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGACCGCTCCCACAAGTTTGCTGAGGTCTCTGTGGGAGCTGGCGTTcccacaaaa                                                                            >  1:15681444/1‑131 (MQ=38)
                                                           ccAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAAGTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGc                                                                         >  2:11469485/1‑137 (MQ=255)
                                                                tccTACATGGCTTCGCGTGGGATTTTTTTTCGTTGACCCATTGCCTGGGCCTTCTTTGCATCCTTTTGTGGGCACCCCCGCCCCCACAGTTTTTCTGAGCTCTCTGTGGGAGTGGACGTGCCCGCGAAGGGCCGCGa                                                                    <  2:14069625/136‑1 (MQ=255)
                                                                 ccTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGACCACTCCCAAAAGGTTGGTGGGAGCTCTGTGGGAGCGGGCGTGCCCACAGAGGGACGCGaa                                                                   >  1:14634247/1‑137 (MQ=25)
                                                                       tGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTTCTTAGATCTCTGTGGGAGCAGGCATGCCCGCAAAAGGCCGCGAAACGGcc                                                             >  1:17432587/1‑137 (MQ=37)
                                                                       tGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGACCGCCCACTCCCACAAGTTTGCTGAGAGCTCTTTGGGGGCGGGCGTGCCCACAAAGGGCCGCGAAGCGAcc                                                             >  1:15971614/1‑137 (MQ=37)
                                                                             cgcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGGACAACCGCCCCCCCAGGTTTGCTGGGATCTCCGTGGGAGCGGGCGTGCCCACGAAAGGCCCCCAAACGGaccccact                                                       >  1:10001700/1‑133 (MQ=255)
                                                                             cgcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGTCCGCTCCCCAAAGTTTTTTGGGGTCTCTTTTGGAGTTGGTGTTCCCGCGAAGGGACGCGAAACGGGCCCaact                                                       >  1:12611020/1‑135 (MQ=19)
                                                                             cgcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGTTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCCAtt                                                       >  2:11045585/1‑137 (MQ=255)
                                                                             cgcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGGGATCTCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCaagt                                                       >  1:11808805/1‑135 (MQ=255)
                                                                             cgcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAAtt                                                       >  2:17986023/1‑137 (MQ=255)
                                                                             cgcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCACGAAAGGCCGCGAAGCAGCCCCaact                                                       >  2:9493176/1‑135 (MQ=38)
                                                                             cgcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCACCCCCCCCAGTTTTTCGAGATCTCTTTTGGAGAGGGCTTGCCCCCGAAGGGCCACCAAACGACCCCaacc                                                       >  1:1277652/1‑135 (MQ=19)
                                                                             cgcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGACCGCTCCCACAAGGTTGGTTAGATCTCTGTGGGGGAGGGAGTGCCCCAAGAGGGCCGCAAAACGGcccccact                                                       >  1:3415800/1‑133 (MQ=21)
                                                                             cgcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACACCCACTCCCCCAAGTTTGCTGAGGTCTCTTTGGGGGTGGGTGTTCCCGCGAAGGGCCACGGAGCGGCCCCAAtt                                                       >  1:9957904/1‑137 (MQ=37)
                                                                             cgcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGAACGACCACCCCCCCAAGTTTTTTGGGATCTCTTTTGGAGCCGGCGTGCCCCCAAAGGACCGCGAAGCAGacccccct                                                       >  1:10660223/1‑133 (MQ=255)
                                                                                gcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACACCCGCTCCCCCAGGTTTTCTGAGATCTCTGTGGGGGCGGGCGTGCCCACAAAAGGCCCCGAAGCGGACCCaaccgtt                                                    >  1:15447885/1‑132 (MQ=37)
                                                                                  gggATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGGGGTCTCTGTGGGAGCGGGCGTGCCCGCAAAGGGCCGCGAAACGGCCCCAATttcttc                                                  >  2:2884654/1‑137 (MQ=37)
                                                                                   ggATCGTTGTTCGCGGCTCCATTGACTGGGTCTGTTTTGCATCCTTTTGCGGGCACTCCCGCCCCAACGGTTTTGTTGAGATCTCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCa                                                 <  2:17683533/137‑1 (MQ=21)
                                                                                   ggATCGTCGTTCGCTGACCCATTAACTGGGGCTGCTTTGCAGCCTTTCGCGGGCACTCCCGCTCCCAAAGGTTTGCTGAGATCTCTGTGGGAGTGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCa                                                 <  2:17273882/137‑1 (MQ=38)
                                                                                               gCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCAAAGGTTTGCTGAGGTCTCTGTGGGAGCGGGCGTGCCCGCAAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCtt                                     >  2:15442020/1‑137 (MQ=255)
                                                                                                 tGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGTCCGTCCGCTCCCACAGGGTTGTTGGGAGCTCTGTGGGGGTGGGCATGCCCACGAAAGGCCGAAGAGCGGCCCCAATTGTTTCAACCGTGGTTTTTgt                                   >  1:11865797/1‑136 (MQ=19)
                                                                                                  gATCCATTGACGGGGTTTTTTTTGCGGCCTTTCGTGGGCACTCCCTCCCCCACGGTTTTGTGGAGCTTTTTGGGGGGGCGGTCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGcc                                  <  1:13063106/137‑1 (MQ=16)
                                                                                                  gATCCATTACCTTGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCACTCCCACAGTTTTGCTGAGGTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGcc                                  <  2:4180191/137‑1 (MQ=25)
                                                                                                              ggggCTGCTTTGCAGCCCTTCGCGGGCACGCCCACTCCCACAGGTTTTCTGAGATCTCCGTGGGAGCGGGCGTGACCGCAAAAGGCCACAAAACGGAACCAACTGCCTCAACCATGGTTCTTGCCGAGCAACGCGTg                      >  1:18037772/1‑137 (MQ=21)
                                                                                                                  ctgctTTGCAGCCCTTCGCGGGCACGCCCACTCCCACAGGTTTGTTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCCATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTaa                  >  2:13211464/1‑137 (MQ=255)
                                                                                                                     ctTTGCAGCCCTTCTCGGGCACGCCCGCTCCCACAGTTTTGCTGAAATCTCTTTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGc               <  2:3363873/137‑1 (MQ=38)
                                                                                                                             cccTTCGCGGGCACGCCCGCTCCCACAAGTTTGCTGAGGTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAACGGCCCCAATTGATTCAACCATGGTTCTTGCCGAGCAGCGCGTGGTAAAACCCGGTATc       >  1:485576/1‑137 (MQ=37)
                                                                                                                             cccTTCGCGGGCACGCCCGCTCCCAAAGGGTTGCTGAGATCTCTGTGGGAGCGGGCGTGTCCGCAAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTTTTGTCAAGCAGCGCGTGGTAAAGCTCGGTATc       >  2:3724391/1‑137 (MQ=38)
                                                                                                                             cccTTCGCGGGCACGCCCACTCCCACAAGTTTGTTGAGATCTCTGTGGGAGAGGGCGTGCCCACGAAAGGACGCGAAAAGGACCCAACTGCTTCCACCACGGTTCTTTCCAAGCAGCGCGTGGTAAAGCCCGGTATc       >  2:6352924/1‑137 (MQ=255)
                                                                                                                                 tCGCGGGCACTCCCTCTCCCACAGGTTTGCTGAGCCCCCTGTGGGAGTGGGTGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCg   <  2:5375154/137‑1 (MQ=38)
                                                                                                                                  cgcgGGCACGCCCGCTCCCACCGGTTTGCTGAGAGCTCTGGGGGAGCAGGCATGCCCGCGAAGGGCCCAGAAGCGGACCCAACTGCTTCAACCGTGGTTCTTGCCGAACAGCGTGTGGTAAAGATCGGTATccccaa  >  2:5791922/1‑135 (MQ=255)
                                                                                                                                     |                                                                                                                                      
TGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGA  >  AE015451_phaM1/352934‑353200

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