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.

GTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTG‑T‑T‑T‑CGGGCAC‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAG  >  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)
         cGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgtg‑t                                                                                                                          <  2:668486/137‑1 (MQ=255)
         cGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgtg‑t                                                                                                                          <  1:9424266/137‑1 (MQ=255)
         cGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCAGCTTTGCCGgtg‑t                                                                                                                          >  2:4809709/1‑137 (MQ=255)
            tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCCCGAAACAAGCGCCACAGTTCATTGCCCCAGCTTTTCCCGTG‑t‑t‑tt                                                                                                                     >  1:4971436/1‑136 (MQ=255)
            tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCACGCAAAAGGAGCCTCAGGGCATGGCCCCAGATTTTCCGGTG‑t‑tct                                                                                                                      >  1:13026964/1‑135 (MQ=255)
            tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGTCCACGAAACAAGCGCCACGGTGCCTGGCACCGGCTTTGCCTGTG‑t‑tct                                                                                                                      >  1:11701029/1‑135 (MQ=255)
            tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTGCACGGCACCGGCTTTGTCGGTG‑t‑tcg                                                                                                                      >  2:6301126/1‑135 (MQ=37)
            tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCATGCCCGCAAAACAGGCGCCGAGGTGCATGGCACCGGCTTTGCCGGTG‑t‑tcg                                                                                                                      >  2:150466/1‑135 (MQ=38)
            tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAAGCGCCACGGTGCACGGCACCGGATTTTTCGGTG‑t‑t‑tt                                                                                                                     >  1:2853386/1‑136 (MQ=255)
              aaGCCCGTGATGTGGCCATTCAGAGCCTGAACGAAGCATTGGCACCACCGCAGGACTGAGGCTGCCAAACGGTAGAAGCGGGCGTGTCCACGGAACAGGCGCCACCGTGCACGGCACCAACTTTGCCGGTG‑t‑tcgcg                                                                                                                    >  2:6108530/1‑133 (MQ=255)
                      gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGTGTGCCCGAGAAAAAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAa‑g‑cc                                                                                                         >  2:2441653/1‑137 (MQ=25)
                      gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCACCGGATTTGCCGGTG‑T‑TCG‑CGGGCAc‑g‑cc                                                                                                         >  2:16218334/1‑137 (MQ=25)
                      gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCCCTGCTTTGTCGGTG‑T‑TCG‑CGGGCAc‑g‑cc                                                                                                         >  2:10817124/1‑137 (MQ=37)
                      gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAGGCGCCACGGGGTATGGCACCCGCTTTTTCCGTG‑T‑TCT‑CGGGCAcaac                                                                                                           >  1:5258733/1‑134 (MQ=255)
                      gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAAGCGCCACAGTGTATGGCACCAGCTTTGCCGGTG‑T‑TCT‑CGGGCAc‑a‑cc                                                                                                         >  2:14170624/1‑137 (MQ=25)
                      gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAAAAAGAGCCGCGGTGCATGGGACCAGATTTTCCGGTG‑T‑TCT‑TGGGCAc‑a‑cc                                                                                                         >  1:10175655/1‑137 (MQ=255)
                       aTGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTTCCCGCGAAACAGGCGCCGCGGGGCATGGCACCGGCTTTGCCAGTG‑T‑TCG‑CGGGCAc‑g‑ccc                                                                                                        >  2:9600394/1‑136 (MQ=37)
                       aTGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCACCGGCTTTGTCGGTG‑T‑TCG‑CGGGCAc‑g‑ccc                                                                                                        >  2:3984472/1‑136 (MQ=25)
                        tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGGCCGCGAAACAAGAGCCACGGTGCCTGGCCCCGGGTTTTCCGGTGTT‑T‑T‑CGGGCAc‑gaccg                                                                                                        >  1:4777954/1‑133 (MQ=255)
                             ccATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCACCGGCTTTGCCGGTG‑T‑T‑CGCGGGCAC‑GCCCGCTccc                                                                                                   >  2:6751384/1‑137 (MQ=21)
                             ccATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGTCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTccc                                                                                                   >  2:3017836/1‑137 (MQ=21)
                             ccATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTccc                                                                                                   >  2:10071397/1‑137 (MQ=21)
                             ccATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTGCATGGCACCGGGTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTccc                                                                                                   >  2:15297598/1‑137 (MQ=18)
                             ccATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCAAAACAAGCGCCATGGTGCACGGCACCAGCTTTTTCGGTGTT‑T‑T‑CGGGCAC‑GCCCGCTccc                                                                                                   >  1:13872817/1‑137 (MQ=255)
                                    gatccTGAACGAAGCTTTGGCACCACAGCAGGACGGAGGCTGCCAAGCGGTGGGAGCGGGCGTGCCCGTGAAACATGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑T‑CGCGGACAC‑GCCCGCTCCCACAGACt                                                                                            <  2:15394393/134‑1 (MQ=17)
                                    gagCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACCGGCGCCACGGTGCATGGCACCAGTTTTGCCGGTG‑T‑TCG‑CAGGCACAC‑CCGCTCCCAAAGACt                                                                                            >  2:9820211/1‑137 (MQ=21)
                                    gagCCGGAACGAGGCCTTGGCACCACAGCAGGACTAAGGCTGCCAAACGGTGAGAGCGGGCGTGCCATCGAAACAGGCGCCGCGGTGCATGGCACCGTCTTTGCCGGTG‑T‑T‑CGCGGGCAC‑GCCCGCTCCCACAGACt                                                                                            <  2:5469515/137‑1 (MQ=18)
                                          gAACGAAGCCTTGGCACCACCGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCACTCCCCCAGAATGATTGc                                                                                      >  2:7079614/1‑137 (MQ=18)
                                          gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGc                                                                                      >  2:10242089/1‑137 (MQ=21)
                                          gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAGAATGATTGc                                                                                      >  2:3069670/1‑137 (MQ=21)
                                          gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑T‑CGCGGGCAC‑GCCCGCTCCCAAAGACTGATTGc                                                                                      >  2:12613622/1‑137 (MQ=21)
                                          gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAAAAGGCGCCACGGTGCCTGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAAACTGATTGc                                                                                      >  1:14132766/1‑137 (MQ=18)
                                              gaggtcTTGGCACCACAGAAGGACGGAGGCTGCAAAATGGTGGGATTGGGCGTTCCCGCGAAACAGGTTCCGTGGGGCATGGAACTGGTTTGGTTGGTG‑T‑TCG‑TGGGCAC‑GCCCGCTCCCACAGACTGATTGCCAtt                                                                                  <  2:8813186/132‑1 (MQ=255)
                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGTGGGCGTGCCCGCGGAACAGGCGCCGAGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGATCCCACGGGCTGAGTGCCagt                                                                                  >  2:11254793/1‑135 (MQ=18)
                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCAAAACAGGCGCCGCGGGGCAGGGCACCGGCTTTGTCGGTG‑T‑TCG‑CGGGTAC‑GCCCACCCCCACAAACCAATTGCCact                                                                                  >  1:15265554/1‑135 (MQ=255)
                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCACGAAACAAGCGCCGCGGTGCATGGCACCGGATTTGCCGGTG‑T‑TCG‑CGGGGAA‑GCCCGCTCCCACAGAATGAGTGACAtt                                                                                  >  1:6651840/1‑137 (MQ=255)
                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGTCCCCGAAACAGGAGCCACGGTGTCTGGCCCCAGCTTTGCCGGTG‑T‑T‑TGCGGGCCCAC‑CCACTCCCCCAGACCGGTTGCCact                                                                                  >  1:13507491/1‑135 (MQ=255)
                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGTGCCGCGGTGCATGGCACCAGATTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCACTCCCACAGAATGACTGCCact                                                                                  >  2:5352466/1‑135 (MQ=18)
                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCact                                                                                  >  2:3671859/1‑135 (MQ=21)
                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCAtt                                                                                  >  2:7218445/1‑137 (MQ=18)
                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCACCACGGTGCCTGGCACCCGCCTTTCCGGTG‑T‑TCT‑CGGGCACAC‑CCGCCCCCACAAACCGATTTcccct                                                                                  >  2:5264328/1‑134 (MQ=16)
                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAAGCGCCGCAGTGCATGGCACCAGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCACCCCCACCGAATGATTGCCact                                                                                  >  1:12233814/1‑135 (MQ=17)
                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCAAAATGGTTGAAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCGGGGCACCGGCTTTGCCGGTG‑T‑T‑CGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCAtt                                                                                  <  2:15150973/137‑1 (MQ=18)
                                              gAAGCCTTGGCACCAAAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGGGCCGCGGTGCCTGGCCCCGGCTTTGCCGGTG‑T‑T‑TGCGGGCCCTC‑CCCCTCCCACAGACTGATTGCCact                                                                                  >  2:3046908/1‑135 (MQ=21)
                                                  ccTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAGGCGCCGCGGTGCACGGCACCAGGTTTGCCAGTG‑T‑TCG‑CAGGCAC‑GCCCGCTCCCAAAGAATGATTGCCATTAAGc                                                                              >  2:2832541/1‑137 (MQ=255)
                                                  ccTTGGCACCACAACAGGACTAAGACTGCAAAATGGGGGGAGTGGGTGTGCCCGAAAAACAGGCGCCGTGGTGAGTGGAACTGGTTTTGCTGGTG‑T‑T‑CGTGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGc                                                                              <  1:7558198/137‑1 (MQ=255)
                                                    tttccaccacaccaGAACGGAGCCTACCACCCGTTGGGAGCGGGCGTCCCCGCGAAACGGGCTCCGCGGTGCATGTCACCGTCTTTGCCGGTT‑T‑TCT‑GGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                            <  2:4892964/133‑1 (MQ=12)
                                                    ttGGCACCCAACCAGGACTGATCCGGCCAAATGGTGGGAGCGGGCGTACCCGCGACACATGCACCGCAGTCCATAACACCGTCTTTTCCGGTG‑T‑TAA‑AGGACAT‑ACCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                            <  2:2015233/137‑1 (MQ=255)
                                                    ttGGCACCACACCAGGACTGAGGCGGCAAAACGGGGGGAGGGGGCGTGCCCGAGAAAGAGGTGCCGTGGGGCATGGCACCGGTTTTGCCGGTG‑T‑T‑CGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                            <  2:6652756/137‑1 (MQ=16)
                                                       acaccacAGAGGGATGGAGGCGTCCAAAAGGGGGGAGGGGGGGGGCCCGTAAAACGGGAGCCGTGGGGAATGGCATTGGTTTTGCGGGGG‑TGC‑G‑TGGACAG‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                         <  2:2542474/136‑1 (MQ=255)
                                                        caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGTGTGTCCGCGAAAAAAGCGCCACGGTGCCTGGCACCAGCTTTTTCGGTG‑T‑TCG‑CGGGTCC‑GCCCCCTCCCACAAAATGATTGCCACTCAACGGGACt                                                                        >  1:7001934/1‑137 (MQ=255)
                                                        caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGTCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGTCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCt                                                                        >  2:6352587/1‑137 (MQ=21)
                                                        caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAGGCGCCGCGGTGCCTGGCACCAGATTTGCCGGTG‑T‑TCG‑CAGGCAC‑GCCCACTCCCACAGAATGATTGCCATTCAACGGGCCt                                                                        >  1:10451558/1‑137 (MQ=255)
                                                        caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGTCCGCGAAACAAGCGCCGCGGTGTACGGCACCGGATTTGCCAGTGTT‑T‑T‑CGGGTAC‑GCCCACCCCCACAAAATGAGTGCCATTGAGAGGGCCt                                                                        >  1:7288381/1‑137 (MQ=255)
                                                          ccacAGCAGGACTGAGGCTGCCAAATGGTGGGAGTGGGTGTACCTGCGAAACAGGCACCGCGGTGCGTGGCACCGGTTTTGCCGGTG‑T‑TCG‑TGGACAC‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTAt                                                                      <  1:13297130/137‑1 (MQ=255)
                                                            acaGCAGGACTGAGTCTTCCAAACGGTGGGGGCGGGCGTGCCCGCGAAACAGTCGCCGCGGTGCAGGGCACTGGCTTTGCCGGTG‑T‑T‑CGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                    <  2:11701272/137‑1 (MQ=18)
                                                                    aCTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCGAAACAAGCGCCACGGTGCCTGGCACCAGCCTTTTCAGTG‑T‑TCG‑CAGGCAC‑GCCCGCTCCCACAAAATGATTGCCACTCAGCGGGACCATCGCTCCAATg                                                            >  2:3006359/1‑137 (MQ=12)
                                                                    aCTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAGGCGCCGCGGTGCATGGCACCGGATTTGCCGGTG‑T‑T‑TGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTCAGCGGGACTATCGCTACAGTg                                                            >  2:7547334/1‑137 (MQ=21)
                                                                     cTGAGGCTGCCAAACGGTGGGAGCGGGTGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCACCGGCTTTGCCGGTG‑T‑TCG‑AGGGCACAA‑CCGCTCCCACAGAATGATTGCCACTCAGCGGGACTACCGATAAAGGGa                                                           >  2:5066090/1‑137 (MQ=12)
                                                                      tGAGGCTACCAACCGGGGGGAGTGGGCGTACCCGCGAAACAGGCGCCGCGGTGCATGACACCGGCTTTGCTGGTG‑T‑T‑CGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAc                                                          <  2:13798245/137‑1 (MQ=21)
                                                                        aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCCCCTGCTTTGTCGGTG‑T‑TCG‑CGGGCCC‑GCCCGCTCCCAAAGACTGATTGCCATTAAGCGGGTCTACCGCTACAGTGACgg                                                        >  2:10411353/1‑137 (MQ=21)
                                                                          gttGCCAAACGGGGGGGGCGGGCGTTCCTGTGAAACGTGCGCCGCGGTGCGTGGCACCGGTTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGcc                                                      <  2:11912032/135‑1 (MQ=18)
                                                                           cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACCGGCGCCGCGGTGCATGGCACCGGGTTTGCCAGTG‑T‑T‑TGCGGGACC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGccc                                                     >  2:1680656/1‑137 (MQ=21)
                                                                           cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGTGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑T‑TGTGGGGAC‑GCCCGCTCCCACAGGCTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAAGGccc                                                     >  2:8942531/1‑137 (MQ=21)
                                                                           cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCACGGTGCACGGCACCAGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGccc                                                     >  2:15911107/1‑137 (MQ=21)
                                                                           cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGAGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAAAACGATTGCCATTAAACGGGCCTATCGCTACAGTGACGGccc                                                     >  2:7192611/1‑137 (MQ=18)
                                                                           cTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAAGCACCGCGGTGCATGGCACCGGATTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCAAAGAATGATTGCAATTAAGCGGGCCGACCACTCAAGTGACGGacc                                                     >  2:6453055/1‑137 (MQ=9)
                                                                           cTGCCAAACGGTGGGAGCGGGAGTGCCCGTGAAACAGGCGCCGCGGTGCATGGCACCGGTTTTGTCGGTG‑T‑TCG‑CGGGCAC‑GACCGCTCCCAAAGACTGAGTGCCATTAAGCGGGACTCGCACTACAGTGACGGacc                                                     >  2:5531945/1‑137 (MQ=255)
                                                                             gCCAAACGGTGGGAGCGGGCGTGCCCACAAAACAAGCGCCACAGTGCACGGCCCCAGCTTTTCCCGTG‑T‑TCT‑CGGGCACAC‑CCACCCCCCCAAAACGATTGCCACTCAACGGGCCCATCACTCCAATGACAGcccat                                                   >  2:4380961/1‑135 (MQ=255)
                                                                                aaaCGGTGGGAGCGGGCGTGCCCGCAAAACAAGAGCCGCGGTGCACGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCCCCGACTGATTGCCATTCAGGGGGACTATCGCCCCCGTGACGGCCCGtttt                                                >  2:5872676/1‑137 (MQ=9)
                                                                                 aaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGTATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTCAGCGGGCCTCTCGCTACAGTGACGGCCCGttttt                                               >  2:15598995/1‑137 (MQ=18)
                                                                                 aaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑T‑CGCGGGAAC‑GCCCACTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGttttt                                               >  2:8667077/1‑137 (MQ=21)
                                                                                      tGGGAGCGGGTGTGCCCGCGAAACGGGCGCCGTGGTGCAGGGCACCGGCTTTGCTGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgttg                                          <  2:3317264/137‑1 (MQ=21)
                                                                                           gCGGGGGTGCCCGCGGAACAGGCGCCACGGTGCATGGAACCGGCTTTGGCGGTG‑T‑TCG‑CGGGAAC‑GCCCGATCCAACAGGCTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCAATTTTCTGTTGCCTCg                                     >  2:407127/1‑137 (MQ=255)
                                                                                            cGGGCTTGCCCGCGAAACAGGTGCCGTGGTGCGTGGCATCGGCTTTGCCGGTG‑T‑T‑TGAGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc                                    <  2:14577654/137‑1 (MQ=25)
                                                                                            agggCGTACCAGAGAAACAGGCGCCGTGGTGCAGGGCACCGGTTTTGCTGGTG‑T‑TCG‑CGGACAC‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc                                    <  1:9071084/136‑1 (MQ=255)
                                                                                                      ctcgAAAAAGGCGCCGTGTTGCTTTGCATCGGCTTTGCCGGTG‑T‑T‑CGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                          <  2:9073355/135‑1 (MQ=25)
                                                                                                      cgcgAAACAGGCGCTGCGGTGCATGGCACCGGCTTTGCCGGGG‑T‑TCG‑CGGGAAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCTTTTTTTGTTGCCTCGCTAACGAATCa                          <  2:10674676/137‑1 (MQ=25)
                                                                                                      cgcgAAACAGGCGCAGCGGTGCATGGCACTGGTTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCAACAGACTGATTGCCATTAAGCGGGCCTATCGTTATAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                          <  2:13179536/137‑1 (MQ=25)
                                                                                                      cgcaAAACATTCGCCGCGGTGCGTGCCCCCGTCTTTTCCGTTG‑T‑T‑CGCGGACAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                          <  1:4407140/137‑1 (MQ=255)
                                                                                                      cgaaaaaCAGGCGCCGTGGTGCGGGGCACTGGTTTTGCCGGTG‑T‑T‑TGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                          <  1:200421/133‑1 (MQ=37)
                                                                                                      cgaaaaaCAGGCGCCGTGGTGCATGGCACCGTCTTTGCTGGTG‑T‑TCG‑CGGACAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGATAACGAATCa                          <  1:14170624/133‑1 (MQ=255)
                                                                                                           aaCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCAAAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGACTCATGACt                     >  2:12225911/1‑137 (MQ=25)
                                                                                                                 cgcttcgGTGCATGGCACAAGCTTTGCTGGTG‑G‑TCG‑TGGGCAC‑GCCCGCTCCCACAGACTGATCGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCGCGCTAACGAATCATGACTCCCTTa               <  2:9145510/132‑1 (MQ=25)
                                                                                                                 cgccgtGGTGCGGGGCACCGGTTTGGCCGGTG‑T‑T‑TGTGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa               <  1:3796010/137‑1 (MQ=255)
                                                                                                                 cgccggGTTGCGTGGCACCGGTTTTTCGGGTG‑T‑T‑CGCGGGCAT‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa               <  2:6252241/137‑1 (MQ=25)
                                                                                                                 cgccgcGGTGCGTGGCACCGGCTTTGCTGGTT‑T‑T‑CGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa               <  2:1061534/137‑1 (MQ=37)
                                                                                                                 cgcagcGGTGCATGGCACTGTCTTTGCGGTTG‑TCC‑G‑CGGACAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa               <  2:8218430/137‑1 (MQ=255)
                                                                                                                   ccgtGGTGCATGGCGCTGTTTTTGCTGGTG‑T‑T‑TGGGGGCAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa             <  1:10483162/137‑1 (MQ=255)
                                                                                                                      cgGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg          >  1:3694123/1‑137 (MQ=37)
                                                                                                                      cgGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg          >  2:13207512/1‑137 (MQ=37)
                                                                                                                          gCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTAt      >  1:6109434/1‑137 (MQ=37)
                                                                                                                           cATGGCACCGGCTTTTCTGGGG‑T‑TCT‑CGGGCAC‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc     <  2:2669011/137‑1 (MQ=255)
                                                                                                                                 ||                                                                                                                                        
GTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTG‑T‑T‑T‑CGGGCAC‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAG  >  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: ‑