Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 1,499,480 2 bp→AC intergenic (+57/‑106) trpS CDS (Tryptophan‑‑tRNA ligase) → / → zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain) –/–

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM11,499,4800CA92.2% 408.6 / NA 194intergenic (+57/‑107)trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)–/–
Reads supporting (aligned to +/- strand):  ref base C (14/1);  new base A (106/71);  total (120/74)
Fisher's exact test for biased strand distribution p-value = 1.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*AE015451_phaM11,499,4810AC97.4% 463.3 / ‑6.3 195intergenic (+58/‑106)trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)–/–
Reads supporting (aligned to +/- strand):  ref base A (4/1);  new base C (116/71);  total (120/75)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 195 total aligned reads displayed.

GTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAG  >  AE015451_phaM1/1499351‑1499611
                                                                                                                                 ||                                                                                                                                        
  tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTg                                                                                                                                  <  1:3964874/137‑1 (MQ=255)
  tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTg                                                                                                                                  <  2:7854379/137‑1 (MQ=255)
  tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTg                                                                                                                                  <  2:16323413/137‑1 (MQ=255)
   cGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGc                                                                                                                                 <  1:15831770/137‑1 (MQ=255)
     acgaCGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCTTTGGCACCACAGCATGACTGGGGCTGCCAAACGTTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg                                                                                                                               <  2:14902189/134‑1 (MQ=37)
      cAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGTCCGCGAAACAGGCGCCGCGGTGTCTGGCACCAGATTTGCCgg                                                                                                                              >  2:427895/1‑137 (MQ=255)
      cAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAACGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGcctg                                                                                                                              >  2:9720241/1‑135 (MQ=255)
         cGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgtgt                                                                                                                          <  2:668486/137‑1 (MQ=255)
         cGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgtgt                                                                                                                          <  1:9424266/137‑1 (MQ=255)
         cGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCAGCTTTGCCGgtgt                                                                                                                          >  2:4809709/1‑137 (MQ=255)
            tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCCCGAAACAAGCGCCACAGTTCATTGCCCCAGCTTTTCCCGTGtttt                                                                                                                     >  1:4971436/1‑136 (MQ=255)
            tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCACGCAAAAGGAGCCTCAGGGCATGGCCCCAGATTTTCCGGTGttct                                                                                                                      >  1:13026964/1‑135 (MQ=255)
            tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGTCCACGAAACAAGCGCCACGGTGCCTGGCACCGGCTTTGCCTGTGttct                                                                                                                      >  1:11701029/1‑135 (MQ=255)
            tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTGCACGGCACCGGCTTTGTCGGTGttcg                                                                                                                      >  2:6301126/1‑135 (MQ=37)
            tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCATGCCCGCAAAACAGGCGCCGAGGTGCATGGCACCGGCTTTGCCGGTGttcg                                                                                                                      >  2:150466/1‑135 (MQ=38)
            tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAAGCGCCACGGTGCACGGCACCGGATTTTTCGGTGtttt                                                                                                                     >  1:2853386/1‑136 (MQ=255)
              aaGCCCGTGATGTGGCCATTCAGAGCCTGAACGAAGCATTGGCACCACCGCAGGACTGAGGCTGCCAAACGGTAGAAGCGGGCGTGTCCACGGAACAGGCGCCACCGTGCACGGCACCAACTTTGCCGGTGttcgcg                                                                                                                    >  2:6108530/1‑133 (MQ=255)
                      gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGTGTGCCCGAGAAAAAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAagcc                                                                                                         >  2:2441653/1‑137 (MQ=25)
                      gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCACCGGATTTGCCGGTGTTCGCGGGCAcgcc                                                                                                         >  2:16218334/1‑137 (MQ=25)
                      gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCCCTGCTTTGTCGGTGTTCGCGGGCAcgcc                                                                                                         >  2:10817124/1‑137 (MQ=37)
                      gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAGGCGCCACGGGGTATGGCACCCGCTTTTTCCGTGTTCTCGGGCAcaac                                                                                                           >  1:5258733/1‑134 (MQ=255)
                      gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAAGCGCCACAGTGTATGGCACCAGCTTTGCCGGTGTTCTCGGGCAcacc                                                                                                         >  2:14170624/1‑137 (MQ=25)
                      gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAAAAAGAGCCGCGGTGCATGGGACCAGATTTTCCGGTGTTCTTGGGCAcacc                                                                                                         >  1:10175655/1‑137 (MQ=255)
                       aTGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTTCCCGCGAAACAGGCGCCGCGGGGCATGGCACCGGCTTTGCCAGTGTTCGCGGGCAcgccc                                                                                                        >  2:9600394/1‑136 (MQ=37)
                       aTGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCACCGGCTTTGTCGGTGTTCGCGGGCAcgccc                                                                                                        >  2:3984472/1‑136 (MQ=25)
                        tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGGCCGCGAAACAAGAGCCACGGTGCCTGGCCCCGGGTTTTCCGGTGTTTTCGGGCAcgaccg                                                                                                        >  1:4777954/1‑133 (MQ=255)
                             ccATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTccc                                                                                                   >  2:6751384/1‑137 (MQ=21)
                             ccATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGTCGGTGTTCGCGGGCACGCCCGCTccc                                                                                                   >  2:3017836/1‑137 (MQ=21)
                             ccATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTccc                                                                                                   >  2:10071397/1‑137 (MQ=21)
                             ccATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTGCATGGCACCGGGTTTGCCGGTGTTCGCGGGCACGCCCGCTccc                                                                                                   >  2:15297598/1‑137 (MQ=18)
                             ccATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCAAAACAAGCGCCATGGTGCACGGCACCAGCTTTTTCGGTGTTTTCGGGCACGCCCGCTccc                                                                                                   >  1:13872817/1‑137 (MQ=255)
                                    gatccTGAACGAAGCTTTGGCACCACAGCAGGACGGAGGCTGCCAAGCGGTGGGAGCGGGCGTGCCCGTGAAACATGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGACACGCCCGCTCCCACAGACt                                                                                            <  2:15394393/134‑1 (MQ=17)
                                    gagCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACCGGCGCCACGGTGCATGGCACCAGTTTTGCCGGTGTTCGCAGGCACACCCGCTCCCAAAGACt                                                                                            >  2:9820211/1‑137 (MQ=21)
                                    gagCCGGAACGAGGCCTTGGCACCACAGCAGGACTAAGGCTGCCAAACGGTGAGAGCGGGCGTGCCATCGAAACAGGCGCCGCGGTGCATGGCACCGTCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACt                                                                                            <  2:5469515/137‑1 (MQ=18)
                                          gAACGAAGCCTTGGCACCACCGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCACTCCCCCAGAATGATTGc                                                                                      >  2:7079614/1‑137 (MQ=18)
                                          gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGc                                                                                      >  2:10242089/1‑137 (MQ=21)
                                          gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGAATGATTGc                                                                                      >  2:3069670/1‑137 (MQ=21)
                                          gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAGACTGATTGc                                                                                      >  2:12613622/1‑137 (MQ=21)
                                          gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAAAAGGCGCCACGGTGCCTGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAACTGATTGc                                                                                      >  1:14132766/1‑137 (MQ=18)
                                              gaggtcTTGGCACCACAGAAGGACGGAGGCTGCAAAATGGTGGGATTGGGCGTTCCCGCGAAACAGGTTCCGTGGGGCATGGAACTGGTTTGGTTGGTGTTCGTGGGCACGCCCGCTCCCACAGACTGATTGCCAtt                                                                                  <  2:8813186/132‑1 (MQ=255)
                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGTGGGCGTGCCCGCGGAACAGGCGCCGAGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGATCCCACGGGCTGAGTGCCagt                                                                                  >  2:11254793/1‑135 (MQ=18)
                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCAAAACAGGCGCCGCGGGGCAGGGCACCGGCTTTGTCGGTGTTCGCGGGTACGCCCACCCCCACAAACCAATTGCCact                                                                                  >  1:15265554/1‑135 (MQ=255)
                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCACGAAACAAGCGCCGCGGTGCATGGCACCGGATTTGCCGGTGTTCGCGGGGAAGCCCGCTCCCACAGAATGAGTGACAtt                                                                                  >  1:6651840/1‑137 (MQ=255)
                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGTCCCCGAAACAGGAGCCACGGTGTCTGGCCCCAGCTTTGCCGGTGTTTGCGGGCCCACCCACTCCCCCAGACCGGTTGCCact                                                                                  >  1:13507491/1‑135 (MQ=255)
                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGTGCCGCGGTGCATGGCACCAGATTTGCCGGTGTTCGCGGGCACGCCCACTCCCACAGAATGACTGCCact                                                                                  >  2:5352466/1‑135 (MQ=18)
                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCact                                                                                  >  2:3671859/1‑135 (MQ=21)
                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCAtt                                                                                  >  2:7218445/1‑137 (MQ=18)
                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCACCACGGTGCCTGGCACCCGCCTTTCCGGTGTTCTCGGGCACACCCGCCCCCACAAACCGATTTcccct                                                                                  >  2:5264328/1‑134 (MQ=16)
                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAAGCGCCGCAGTGCATGGCACCAGCTTTGCCGGTGTTCGCGGGCACGCCCACCCCCACCGAATGATTGCCact                                                                                  >  1:12233814/1‑135 (MQ=17)
                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCAAAATGGTTGAAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCGGGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCAtt                                                                                  <  2:15150973/137‑1 (MQ=18)
                                              gAAGCCTTGGCACCAAAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGGGCCGCGGTGCCTGGCCCCGGCTTTGCCGGTGTTTGCGGGCCCTCCCCCTCCCACAGACTGATTGCCact                                                                                  >  2:3046908/1‑135 (MQ=21)
                                                  ccTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAGGCGCCGCGGTGCACGGCACCAGGTTTGCCAGTGTTCGCAGGCACGCCCGCTCCCAAAGAATGATTGCCATTAAGc                                                                              >  2:2832541/1‑137 (MQ=255)
                                                  ccTTGGCACCACAACAGGACTAAGACTGCAAAATGGGGGGAGTGGGTGTGCCCGAAAAACAGGCGCCGTGGTGAGTGGAACTGGTTTTGCTGGTGTTCGTGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGc                                                                              <  1:7558198/137‑1 (MQ=255)
                                                    tttccaccacaccaGAACGGAGCCTACCACCCGTTGGGAGCGGGCGTCCCCGCGAAACGGGCTCCGCGGTGCATGTCACCGTCTTTGCCGGTTTTCTGGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                            <  2:4892964/133‑1 (MQ=12)
                                                    ttGGCACCCAACCAGGACTGATCCGGCCAAATGGTGGGAGCGGGCGTACCCGCGACACATGCACCGCAGTCCATAACACCGTCTTTTCCGGTGTTAAAGGACATACCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                            <  2:2015233/137‑1 (MQ=255)
                                                    ttGGCACCACACCAGGACTGAGGCGGCAAAACGGGGGGAGGGGGCGTGCCCGAGAAAGAGGTGCCGTGGGGCATGGCACCGGTTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                            <  2:6652756/137‑1 (MQ=16)
                                                       acaccacAGAGGGATGGAGGCGTCCAAAAGGGGGGAGGGGGGGGGCCCGTAAAACGGGAGCCGTGGGGAATGGCATTGGTTTTGCGGGGGTGCGTGGACAGGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                         <  2:2542474/136‑1 (MQ=255)
                                                        caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGTGTGTCCGCGAAAAAAGCGCCACGGTGCCTGGCACCAGCTTTTTCGGTGTTCGCGGGTCCGCCCCCTCCCACAAAATGATTGCCACTCAACGGGACt                                                                        >  1:7001934/1‑137 (MQ=255)
                                                        caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGTCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGTCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCt                                                                        >  2:6352587/1‑137 (MQ=21)
                                                        caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAGGCGCCGCGGTGCCTGGCACCAGATTTGCCGGTGTTCGCAGGCACGCCCACTCCCACAGAATGATTGCCATTCAACGGGCCt                                                                        >  1:10451558/1‑137 (MQ=255)
                                                        caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGTCCGCGAAACAAGCGCCGCGGTGTACGGCACCGGATTTGCCAGTGTTTTCGGGTACGCCCACCCCCACAAAATGAGTGCCATTGAGAGGGCCt                                                                        >  1:7288381/1‑137 (MQ=255)
                                                          ccacAGCAGGACTGAGGCTGCCAAATGGTGGGAGTGGGTGTACCTGCGAAACAGGCACCGCGGTGCGTGGCACCGGTTTTGCCGGTGTTCGTGGACACTCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTAt                                                                      <  1:13297130/137‑1 (MQ=255)
                                                            acaGCAGGACTGAGTCTTCCAAACGGTGGGGGCGGGCGTGCCCGCGAAACAGTCGCCGCGGTGCAGGGCACTGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                    <  2:11701272/137‑1 (MQ=18)
                                                                    aCTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCGAAACAAGCGCCACGGTGCCTGGCACCAGCCTTTTCAGTGTTCGCAGGCACGCCCGCTCCCACAAAATGATTGCCACTCAGCGGGACCATCGCTCCAATg                                                            >  2:3006359/1‑137 (MQ=12)
                                                                    aCTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAGGCGCCGCGGTGCATGGCACCGGATTTGCCGGTGTTTGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTCAGCGGGACTATCGCTACAGTg                                                            >  2:7547334/1‑137 (MQ=21)
                                                                     cTGAGGCTGCCAAACGGTGGGAGCGGGTGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCACCGGCTTTGCCGGTGTTCGAGGGCACAACCGCTCCCACAGAATGATTGCCACTCAGCGGGACTACCGATAAAGGGa                                                           >  2:5066090/1‑137 (MQ=12)
                                                                      tGAGGCTACCAACCGGGGGGAGTGGGCGTACCCGCGAAACAGGCGCCGCGGTGCATGACACCGGCTTTGCTGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAc                                                          <  2:13798245/137‑1 (MQ=21)
                                                                        aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCCCCTGCTTTGTCGGTGTTCGCGGGCCCGCCCGCTCCCAAAGACTGATTGCCATTAAGCGGGTCTACCGCTACAGTGACgg                                                        >  2:10411353/1‑137 (MQ=21)
                                                                          gttGCCAAACGGGGGGGGCGGGCGTTCCTGTGAAACGTGCGCCGCGGTGCGTGGCACCGGTTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGcc                                                      <  2:11912032/135‑1 (MQ=18)
                                                                           cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACCGGCGCCGCGGTGCATGGCACCGGGTTTGCCAGTGTTTGCGGGACCGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGccc                                                     >  2:1680656/1‑137 (MQ=21)
                                                                           cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGTGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTTGTGGGGACGCCCGCTCCCACAGGCTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAAGGccc                                                     >  2:8942531/1‑137 (MQ=21)
                                                                           cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCACGGTGCACGGCACCAGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGccc                                                     >  2:15911107/1‑137 (MQ=21)
                                                                           cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGAGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAAACGATTGCCATTAAACGGGCCTATCGCTACAGTGACGGccc                                                     >  2:7192611/1‑137 (MQ=18)
                                                                           cTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAAGCACCGCGGTGCATGGCACCGGATTTGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAGAATGATTGCAATTAAGCGGGCCGACCACTCAAGTGACGGacc                                                     >  2:6453055/1‑137 (MQ=9)
                                                                           cTGCCAAACGGTGGGAGCGGGAGTGCCCGTGAAACAGGCGCCGCGGTGCATGGCACCGGTTTTGTCGGTGTTCGCGGGCACGACCGCTCCCAAAGACTGAGTGCCATTAAGCGGGACTCGCACTACAGTGACGGacc                                                     >  2:5531945/1‑137 (MQ=255)
                                                                             gCCAAACGGTGGGAGCGGGCGTGCCCACAAAACAAGCGCCACAGTGCACGGCCCCAGCTTTTCCCGTGTTCTCGGGCACACCCACCCCCCCAAAACGATTGCCACTCAACGGGCCCATCACTCCAATGACAGcccat                                                   >  2:4380961/1‑135 (MQ=255)
                                                                                aaaCGGTGGGAGCGGGCGTGCCCGCAAAACAAGAGCCGCGGTGCACGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCCCCGACTGATTGCCATTCAGGGGGACTATCGCCCCCGTGACGGCCCGtttt                                                >  2:5872676/1‑137 (MQ=9)
                                                                                 aaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGTATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTCAGCGGGCCTCTCGCTACAGTGACGGCCCGttttt                                               >  2:15598995/1‑137 (MQ=18)
                                                                                 aaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGAACGCCCACTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGttttt                                               >  2:8667077/1‑137 (MQ=21)
                                                                                      tGGGAGCGGGTGTGCCCGCGAAACGGGCGCCGTGGTGCAGGGCACCGGCTTTGCTGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgttg                                          <  2:3317264/137‑1 (MQ=21)
                                                                                           gCGGGGGTGCCCGCGGAACAGGCGCCACGGTGCATGGAACCGGCTTTGGCGGTGTTCGCGGGAACGCCCGATCCAACAGGCTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCAATTTTCTGTTGCCTCg                                     >  2:407127/1‑137 (MQ=255)
                                                                                            cGGGCTTGCCCGCGAAACAGGTGCCGTGGTGCGTGGCATCGGCTTTGCCGGTGTTTGAGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc                                    <  2:14577654/137‑1 (MQ=25)
                                                                                            agggCGTACCAGAGAAACAGGCGCCGTGGTGCAGGGCACCGGTTTTGCTGGTGTTCGCGGACACTCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc                                    <  1:9071084/136‑1 (MQ=255)
                                                                                                      ctcgAAAAAGGCGCCGTGTTGCTTTGCATCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                          <  2:9073355/135‑1 (MQ=25)
                                                                                                      cgcgAAACAGGCGCTGCGGTGCATGGCACCGGCTTTGCCGGGGTTCGCGGGAACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCTTTTTTTGTTGCCTCGCTAACGAATCa                          <  2:10674676/137‑1 (MQ=25)
                                                                                                      cgcgAAACAGGCGCAGCGGTGCATGGCACTGGTTTTGCCGGTGTTCGCGGGCACGCCCGCTCCAACAGACTGATTGCCATTAAGCGGGCCTATCGTTATAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                          <  2:13179536/137‑1 (MQ=25)
                                                                                                      cgcaAAACATTCGCCGCGGTGCGTGCCCCCGTCTTTTCCGTTGTTCGCGGACACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                          <  1:4407140/137‑1 (MQ=255)
                                                                                                      cgaaaaaCAGGCGCCGTGGTGCGGGGCACTGGTTTTGCCGGTGTTTGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                          <  1:200421/133‑1 (MQ=37)
                                                                                                      cgaaaaaCAGGCGCCGTGGTGCATGGCACCGTCTTTGCTGGTGTTCGCGGACACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGATAACGAATCa                          <  1:14170624/133‑1 (MQ=255)
                                                                                                           aaCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGACTCATGACt                     >  2:12225911/1‑137 (MQ=25)
                                                                                                                 cgcttcgGTGCATGGCACAAGCTTTGCTGGTGGTCGTGGGCACGCCCGCTCCCACAGACTGATCGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCGCGCTAACGAATCATGACTCCCTTa               <  2:9145510/132‑1 (MQ=25)
                                                                                                                 cgccgtGGTGCGGGGCACCGGTTTGGCCGGTGTTTGTGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa               <  1:3796010/137‑1 (MQ=255)
                                                                                                                 cgccggGTTGCGTGGCACCGGTTTTTCGGGTGTTCGCGGGCATGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa               <  2:6252241/137‑1 (MQ=25)
                                                                                                                 cgccgcGGTGCGTGGCACCGGCTTTGCTGGTTTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa               <  2:1061534/137‑1 (MQ=37)
                                                                                                                 cgcagcGGTGCATGGCACTGTCTTTGCGGTTGTCCGCGGACACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa               <  2:8218430/137‑1 (MQ=255)
                                                                                                                   ccgtGGTGCATGGCGCTGTTTTTGCTGGTGTTTGGGGGCACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa             <  1:10483162/137‑1 (MQ=255)
                                                                                                                      cgGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg          >  1:3694123/1‑137 (MQ=37)
                                                                                                                      cgGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg          >  2:13207512/1‑137 (MQ=37)
                                                                                                                          gCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTAt      >  1:6109434/1‑137 (MQ=37)
                                                                                                                           cATGGCACCGGCTTTTCTGGGGTTCTCGGGCACTCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc     <  2:2669011/137‑1 (MQ=255)
                                                                                                                                 ||                                                                                                                                        
GTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAG  >  AE015451_phaM1/1499351‑1499611

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