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.8% 1804.8 / NA 481intergenic (+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 (221/258);  total (222/259)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 482 total aligned reads displayed.

TGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCG  >  AE015451_phaM1/352934‑353199
                                                                                                                                     |                                                                                                                                     
 gCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCa                                                                                                                                   >  1:255751/1‑137 (MQ=255)
   cAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCAca                                                                                                                                 >  1:3346126/1‑136 (MQ=255)
   cAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACg                                                                                                                                 <  1:16293897/137‑1 (MQ=255)
   cAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACg                                                                                                                                 <  1:9190042/137‑1 (MQ=255)
   cAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACg                                                                                                                                 <  1:12208147/137‑1 (MQ=255)
    aGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGc                                                                                                                                <  2:2683863/137‑1 (MQ=255)
         ggTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGc                                                                                                                                    >  1:7693950/1‑128 (MQ=255)
          gTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTc                                                                                                                          <  2:2300315/137‑1 (MQ=255)
          gTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTc                                                                                                                          <  2:5652897/137‑1 (MQ=255)
          gTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTc                                                                                                                          >  2:13369792/1‑137 (MQ=255)
          gTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGACACGCCCGCTc                                                                                                                          <  2:4317846/137‑1 (MQ=255)
              tGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGTGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTACAGCCCTTCGCGGGCAAGCCCGATCCcac                                                                                                                      >  2:9288935/1‑137 (MQ=255)
               ggcggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCcaca                                                                                                                     <  1:6974618/137‑1 (MQ=255)
               ggcggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCcaca                                                                                                                     <  2:8922935/137‑1 (MQ=255)
               ggcggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCcaca                                                                                                                     <  2:17966904/137‑1 (MQ=255)
               ggcggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCcaca                                                                                                                     <  1:2143351/137‑1 (MQ=255)
                 cggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAga                                                                                                                   <  1:18074106/137‑2 (MQ=255)
                   gcCAGCATCTTCCACTTTGGCGAGTCCACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGtt                                                                                                                 >  2:4608055/1‑137 (MQ=255)
                   gcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGtt                                                                                                                 <  1:8379639/137‑1 (MQ=255)
                     cAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTg                                                                                                               <  2:7430316/137‑1 (MQ=255)
                     cAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCGGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTg                                                                                                               <  2:2979599/137‑1 (MQ=255)
                       gCATCTTCCACTTTGGCGAGTACACAGTGCCTGAGGCCAATGCCTACATGACTTCGCGCGGGATCGTCGTTCGCTGATCCATTGATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                             <  2:18145592/137‑1 (MQ=38)
                       gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                             <  1:416446/137‑1 (MQ=255)
                       gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                             <  1:17955221/137‑1 (MQ=255)
                       gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                             <  1:16070353/137‑1 (MQ=255)
                       gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                             <  2:5265443/137‑1 (MQ=255)
                       gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                             <  1:13580855/137‑1 (MQ=255)
                       gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                             <  2:11700820/137‑1 (MQ=255)
                       gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                             <  2:14792070/137‑1 (MQ=255)
                       gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                             <  1:11792175/137‑1 (MQ=255)
                       gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                             <  2:13291700/137‑1 (MQ=255)
                       gCATCTTCCACTTCGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTGTGCt                                                                                                             <  2:9917819/137‑1 (MQ=255)
                       gCATCTGCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGTGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                             <  2:953643/137‑1 (MQ=255)
                        cATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGACCGCTCCCAAAAGTTTGCTg                                                                                                            >  2:12609468/1‑137 (MQ=255)
                         aTCTTCCACTTTGGCGAGTACATAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTTGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  2:2121971/137‑1 (MQ=255)
                         aTCTTCCACTTTGGCGAGTACACAGTGCCTGGGGTCAAGGCCTACGTGGCTTCGCGCGGGGTCGGCGTTCGCTGATCCATTGACCGGGGCGGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTGGCTga                                                                                                           <  2:16772317/137‑1 (MQ=37)
                         aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTTGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGTTTTGCAGCTCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  1:12830327/137‑1 (MQ=255)
                         aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  1:3740855/137‑1 (MQ=255)
                         aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  2:11257871/137‑1 (MQ=255)
                         aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  2:13006243/137‑1 (MQ=255)
                         aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  2:851518/137‑1 (MQ=255)
                         aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  1:14430457/137‑1 (MQ=255)
                         aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  1:8080513/137‑1 (MQ=255)
                         aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  1:17458894/137‑1 (MQ=255)
                         aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  2:14080455/137‑1 (MQ=255)
                         aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  1:5506720/137‑1 (MQ=255)
                         aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCCGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  2:15040433/137‑1 (MQ=255)
                         aTCTTCCACGTTTGCGAGTACACATTGCCGGAGGCCAAGGCCTCCATGGCTTCGCGCGTGATCGCCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCACTCCCACAGGTTTGCTga                                                                                                           <  2:15401134/137‑1 (MQ=37)
                           cTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTgaga                                                                                                         <  1:7869407/137‑1 (MQ=255)
                             tCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGAtc                                                                                                       <  2:2512536/137‑1 (MQ=255)
                               cACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGAtctc                                                                                                     <  1:4474975/137‑1 (MQ=255)
                               cACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGAtctc                                                                                                     <  2:7203687/137‑1 (MQ=255)
                               cACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGAtctc                                                                                                     <  2:3592045/137‑1 (MQ=255)
                                aCTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCCCAGGTTTGGTGAGGTctct                                                                                                    >  1:9657755/1‑137 (MQ=255)
                                aCTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCAAAGGTTTGCTGGGGTctct                                                                                                    >  2:1097988/1‑137 (MQ=255)
                                     ggCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCCCCCACAGGTTTGCTGAGATCTCTGTggg                                                                                              >  1:13011620/1‑137 (MQ=255)
                                       cGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           >  1:10290940/1‑123 (MQ=255)
                                       cGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCAAAGGTTTGGTGAAATCTCTGTGGGAg                                                                                            >  1:14057348/1‑137 (MQ=255)
                                       cGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACACCCCCCCCCCCAAAGTTGGTGAGGCCTCTGTTGGAg                                                                                            >  1:367709/1‑137 (MQ=255)
                                        gAGTACACAGTGCCGGAGGCCAAGGCATACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACGGGTTTGCTGAGATCTCTGTGGGAGc                                                                                           >  2:4052161/1‑137 (MQ=255)
                                                   gCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGACCACCCCCACAAGTTTGCTGAGATCTCCGTGGGAGAGGGCATGCCcg                                                                                >  1:18413252/1‑137 (MQ=38)
                                                   gCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGACCGACCACTCCCCCAGGTTTGCTGGAAGCTCTGTGGGAGCGGGCATGCCcg                                                                                >  1:10878258/1‑137 (MQ=255)
                                                     cGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGTACGCCCGCTCACACAGGTTTGCTGAGATATCTGTGGGAGCGGGCGTGCCcgca                                                                              >  2:6885747/1‑136 (MQ=39)
                                                       gAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTTAGATCTCTGTGGGAGCGGGCGTGCCCGCGaa                                                                            >  2:18370875/1‑137 (MQ=255)
                                                       gAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGGTTTCTGAGGTCTCTGTTGGAGCGGGTGTGCCCGCGaa                                                                            >  2:5079310/1‑137 (MQ=38)
                                                       gAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACACCCGCTCCCACAAGTTTTATGAGATCTCTGTGGGAGAGGGTGTGCCCGCGaa                                                                            >  1:17071538/1‑137 (MQ=38)
                                                       gAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACACCCACTCCCACAAGGTTGGTGAGGCCTCTGTGGGGGCGGGAGTGCCCACGaa                                                                            >  1:6292485/1‑137 (MQ=37)
                                                        aGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGGTCTCTGTGGGAGCGGGCGTGCCCGCGAAg                                                                           >  2:3408224/1‑137 (MQ=255)
                                                        aGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCCGTGGGAGCGGGCGTGCCCGCGAAg                                                                           >  1:13845514/1‑137 (MQ=255)
                                                          gCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                             >  1:11487334/1‑102 (MQ=255)
                                                               ggCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCcgcg                                                                    >  2:18965762/1‑137 (MQ=255)
                                                                       tGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGTACGCCCGCTCCCACACGTTTGTTTAGGGTTCTTTGGGGGCCGGCGTGCCCGCGAAGGGACGCAAAgaaggc                                                            >  1:14961228/1‑132 (MQ=25)
                                                                        ggCTTCGCGCGGGATCGTCGTTCGCTGATCCCTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACCGGTTTGCTGAGATTTCTGTGGGAGCGGGCGTGTCCAAGAAGGGCCGCGAAGTGGacc                                                           >  2:10160413/1‑137 (MQ=37)
                                                                             cgcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTTCTGAGGTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAAtt                                                      >  2:12112765/1‑137 (MQ=255)
                                                                             cgcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGGTCTCTTTGGGAGCGGGCGTGCCCGCGAAGGGTCGCGAAGCGGCCCCCAtt                                                      >  1:13331163/1‑137 (MQ=38)
                                                                             cgcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGGTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAAtt                                                      >  2:8608238/1‑137 (MQ=255)
                                                                             cgcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGGTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAAtt                                                      >  2:14299065/1‑137 (MQ=255)
                                                                             cgcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCATGCCCACAAAGGGCCACGAAGCGGCCCCAAtt                                                      >  2:6400442/1‑137 (MQ=38)
                                                                                gcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGc                                                                                                                                >  1:15391227/1‑61 (MQ=255)
                                                                                gcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCAAAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCt                                                   >  2:1642847/1‑137 (MQ=255)
                                                                                  gggATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGGAAGACCCCCCCCCCAAGGTTTTTGGGGGCTCTGTGGGAGAGGGGGTGCCCACAAAAGGCCACAAAACGGACCCAActtcttt                                                 >  1:110417/1‑136 (MQ=255)
                                                                                  gggATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  1:7282140/80‑1 (MQ=255)
                                                                                  gggATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCCCTGTGGGAGCAGGCGTGCCCGCGAAGGACCGCGAAGCGGCCCAAATTGCTTc                                                 >  2:7602888/1‑137 (MQ=38)
                                                                                   ggATCGTCGTTCGCGGACCGGTGCACGGGGGCGGTTTTGCAGCCCTGGGAGGACACTCCCGCCCACAAGGGTTTGTGGAGATCTTGTGGGGAGGGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCa                                                <  1:10013827/137‑1 (MQ=12)
                                                                                                      cATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGTGGGCGTGCCCGCGAAGGGCCTCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagc                             >  2:16222436/1‑137 (MQ=255)
                                                                                                       gttGACTGGGGCTGTTTTACATCCTTTCGCGGGCACGCCCGCTCCCACGGGTTTGTGGAGATCTTTGGGGGAGTGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagca                            <  2:11846675/136‑1 (MQ=37)
                                                                                                       atggACGGGGGTTGTTTGGCAGCCTTTGGGGGGCACGCCCTCTCCCACGGGTTTGTGGAGCCCTTTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagca                            <  1:16663346/134‑1 (MQ=21)
                                                                                                         tGACTGGGGTTGTTTTGAAGCCTTTCGTGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCGGGGGGAGTGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagcagc                          <  2:18732576/137‑1 (MQ=37)
                                                                                                         tGACTGGGGCTGCTTTGCATCCCTTCGCGGGCACGCCCACTCCCACGGTTTTGCTGAGATCTTTGGGGGGGGGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagcagc                          <  2:6144167/137‑1 (MQ=37)
                                                                                                                 gctgctTTGCAGCCCTTCGCGGACACGCCCGCTCCCACAGGTTTCCTGAAACCTCTTTGGAGGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTa                  <  2:8150236/137‑1 (MQ=38)
                                                                                                                   tgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCCACCGTGGTTCTTGCCAAGCAACGCGTGGTaaa                >  2:10582835/1‑137 (MQ=255)
                                                                                                                   tgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGGTTGTTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGGAGGGCCACGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTaaa                >  2:17043834/1‑137 (MQ=38)
                                                                                                                     ctTTGCGTCCCTTCGTGGGCACGCCCNCTCCCACAGGTTTGCGGAGCTCTTTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGc              <  1:18741210/137‑1 (MQ=38)
                                                                                                                     ctTTGCAGCCCTTTGCGGGCACACCCGCTCCCACAGGTTTGCTGAGATCTCTTTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGc              <  2:15289029/137‑1 (MQ=255)
                                                                                                                       ttGCATCCCTTCGTGGGCACGCCTGCCCCCACAGGTTTGCGGAGGCCTTTGTGGGAGTGGGCGGGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTc            <  2:17554537/137‑1 (MQ=37)
                                                                                                                        tGCAGCCCTTCGCGGGCACACCCGCTCCCAAAGGTTTGCTGAGATCTCTGTGGGAACGGGCGTGCCCACGAAGGGCCACGAAGCGGCCCCAATTGCTTCCACCGTGGTTCTTGCCGAACAGCGCGTGGTAAAGTTCg           >  2:17470670/1‑137 (MQ=255)
                                                                                                                            cctcTTCGGGGGAACGCCCGCTCCCAAGGGTTTGTGGAGATTTTTGTGGGGGTGGTTGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTAt       <  1:3481337/134‑1 (MQ=255)
                                                                                                                             cccTTCGCGGGCACACCCGCTCCCACAAGGTTTCTGAGAGCTCTGTGGGAGCGGGAGTGCCCGCGAAGGGACACGAAACGGCCCCAACTGCTTCAAACATGGTTCTTGCCGAGCAGCGCGGGGTAAAGCTCGGTATc      >  2:14621454/1‑137 (MQ=255)
                                                                                                                              ccTTCGCGGACACTCCAGCTCCCACGGGTTGGTGGAGCTTTTTTGGGGAGGGGGTGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATcc     <  1:8673056/137‑1 (MQ=25)
                                                                                                                                 tCGCGGGCACGCCTATTCCCACAGGTTTACGGAAACCTCTGTGGGAGCGGTCTTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCg  <  1:2978642/137‑1 (MQ=37)
                                                                                                                                     |                                                                                                                                     
TGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCG  >  AE015451_phaM1/352934‑353199

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