Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 1,932,640 +C intergenic (‑107/+43) rluA CDS (23S rRNA (746) and tRNA (32) pseudouridine synthase) ← / ← minE CDS (Cell division topological specificity factor) –/–

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM11,932,6381.C92.9% 1755.0 / 130.2 497intergenic (‑105/+45)rluA CDS (23S rRNA (746) and tRNA (32) pseudouridine synthase)/minE CDS (Cell division topological specificity factor)–/–
Reads supporting (aligned to +/- strand):  ref base . (31/4);  new base C (265/192);  total (301/196)
Fisher's exact test for biased strand distribution p-value = 2.38e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 523 total aligned reads displayed.

TCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑G‑CC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGAT  >  AE015451_phaM1/1932502‑1932769
                                                                                                                                               |                                                                                                                                    
  gAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGGAAGCCCACTCCTCCAAAAGATAGTCGTCCCCTGTAGGA‑G‑A‑GGG‑C‑TT‑T‑cc                                                                                                                                    >  1:4607107/1‑137 (MQ=255)
   aaGAGGATCTGGACATTCGACATCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑G‑ccc                                                                                                                                   <  2:11224683/137‑2 (MQ=255)
   aaGAGAATCTGGACATTCGACAGCGGCATTACAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑G‑ccc                                                                                                                                   <  2:17121540/137‑2 (MQ=255)
      agGATTTGGACGTTCGACATTGGCATTGCAGGTCTAACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGTGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCTTCCCCGGGAGGA‑G‑C‑GGG‑C‑TT‑G‑cccgcg                                                                                                                                <  1:17786108/137‑5 (MQ=37)
           cTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTCCAAAAGATATGCGTCCCCTGTCGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAggg                                                                                                                           >  2:4550995/1‑137 (MQ=38)
           cTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCACTCCCACAAAAGATATGCGTCCCCTGTGAGA‑G‑C‑GGG‑C‑TT‑GCCC‑GAGAAggg                                                                                                                           >  1:9239519/1‑137 (MQ=38)
            tGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGTGGGAAAACCCGCTCCTACAAAGGAGATGCGTCCCCCTTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAaggg                                                                                                                          >  1:3866536/1‑136 (MQ=255)
             ggACATTCGACATCGGCATTGCAGGTCTCACGAAAATAGCGGGAAACGTTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCTTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCt                                                                                                                         <  2:11773777/137‑1 (MQ=37)
             ggACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGTGGGAAAGCCCCCTCCTCCCAAAGAGATGGGTTCCCCGTGGGAGG‑G‑GGG‑T‑TT‑C‑CC‑GCGAAaggat                                                                                                                         >  1:5894129/1‑135 (MQ=255)
             ggACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTCCAAAGGATATGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCt                                                                                                                         >  2:18099402/1‑137 (MQ=255)
                      aCAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTTGCGGGCAAACCCACTCCCCCCAAAGAAATGCATCCCCCGTCGGA‑G‑A‑GGGAT‑TT‑C‑CC‑GCAAAAGGATGCaaagaaa                                                                                                                >  1:11021861/1‑132 (MQ=255)
                      aCAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGGAAAACCACCCCCCCAAAAGAGATGCGTCCCCCGTAGGG‑GAG‑GGG‑T‑TT‑C‑CC‑ACGAAAGGCTGTaaaggag                                                                                                                >  1:7244282/1‑132 (MQ=255)
                      aCAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGTCCGCTCCTACAAAGGATGTGCGTCCCCTGTAGGG‑G‑CGGGG‑T‑TT‑T‑CC‑TAGAAGGGCTGTAAGGGAg                                                                                                                >  2:8882203/1‑137 (MQ=37)
                      aCAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTCCCAAGGATATGCGTCCCCTGTAGGA‑G‑C‑GGG‑CTTT‑C‑CC‑GCGAAGGGCTGCAAAGCAg                                                                                                                >  2:1309393/1‑137 (MQ=37)
                      aCAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCAGGCAAACCCGCTCCCCCAAAAGAGCTGCGTCCCCCTTCGGA‑G‑C‑GGGTT‑TT‑C‑CC‑ACAAAGGGCTTCAAAGCaa                                                                                                                >  2:16671605/1‑136 (MQ=255)
                       cAGCGGAATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAGGCCCGCTCCTACAAAGGATATACGTCCCCTGGAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGc                                                                                                               <  2:12320636/137‑1 (MQ=38)
                          gacagttGCAGGTCCCACAAAAATAGCAGGAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGcccc                                                                                                            <  2:15546497/132‑1 (MQ=37)
                                tGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTAAAAAGGATATGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCCAATTg                                                                                                      >  2:17000982/1‑137 (MQ=255)
                                 gCAGGTCTCACGAAAATAGCAGGAAACGCTGGGGCCGCTTTGCGGCCCTTCGTGGGCAAGCCCGCTCCTAAAAAGGATATGCTTCCCTTGGAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGt                                                                                                     <  1:4066339/137‑1 (MQ=38)
                                      tctcACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACt                                                                                                <  2:10200757/137‑1 (MQ=255)
                                      tctcACGAAAATAGCAGGAAACGCTGGGTCCGCTTTGCGGCCCTTCGCGGACAAACCCGCTCCTACAAAGGATAGGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACt                                                                                                <  1:17570637/137‑1 (MQ=38)
                                      tctaaagAAAGTAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGAAAGCCCGCTCCGAAAAAGGAGAGGCGTCCCTTGTGGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACt                                                                                                <  1:12249130/131‑1 (MQ=25)
                                      tctaaagAAAATAGCAGGAAACGTTGGGGCCGTTTGGCGGCCTTTCGCGGGCAAGCCCGCTCCTACAAAGGAGATGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACt                                                                                                <  1:3707089/131‑1 (MQ=37)
                                        tcACGAAAATAGCAGGAAACGCGGGGGCCGCTTTGTGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATGTGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTAc                                                                                              <  1:16906516/137‑1 (MQ=38)
                                          aCGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAAACCGCTCCCCCAAAAGATATGTGTCCCCCGTAGGA‑G‑C‑GGG‑CTTT‑C‑CC‑GCGAAAGGATGGAAGGCAGCCCCAAACTGTTGACTAcat                                                                                            >  1:2858469/1‑135 (MQ=25)
                                             aaaaTAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCCCCCCCAAAGGATATGCGTCCCCTGTAGGA‑G‑C‑AGG‑C‑TT‑GCCC‑ACGAAGGGCTGCAAGGCAACCCCAAATTGTCGACTCcgaaga                                                                                         >  2:17763409/1‑133 (MQ=37)
                                             aaaaTAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAAACCACTCCCCCAAAAGATATGCGTCCCCCGTAGGA‑GAG‑GGG‑T‑TT‑C‑CC‑ACGAAAGGCTGCAAAGCAACCCCAAACTTTCTCCTACACGGa                                                                                         >  2:1507743/1‑137 (MQ=255)
                                             aaaaTAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGAAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑ACGAAGGGCTGCCAGGCAGCCCCAAATTGTCGACCACACGGa                                                                                         >  2:6408636/1‑137 (MQ=255)
                                                       tAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAACCCGCTCCCCCAAAGGACCTGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGAAGCCCCAAACTGTCGAACACACGGATCCGCGATcc                                                                               >  2:2745385/1‑137 (MQ=25)
                                                       tAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGAAAACCCACCCCCCCCAAAGAAACGCGTCCCCCGTAGGA‑G‑C‑GGG‑T‑TT‑GCCC‑ACGAAGGGATGCAAAGCAGACCCAAACTGTCTACCCCGCGGATCAGCGATcc                                                                               >  1:16535350/1‑137 (MQ=255)
                                                                  gCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTTTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGAt                                                                    <  1:12893434/137‑1 (MQ=255)
                                                                   ccGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATg                                                                   <  2:12817425/137‑1 (MQ=255)
                                                                   ccGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGGGGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATg                                                                   <  1:1910079/137‑1 (MQ=255)
                                                                         tGCGGCACTTCGCGGGAAAGCCCGCTCCTCCAAAAGATATGCGTCCCCTGTCGGA‑G‑C‑GGG‑T‑TT‑GTCC‑GCGGAGGGCTGCAAGGCAGCCCCCAATTTTCGACTCCGCGGATCCGCGATCCGGCAGGGTGATGTTCCGc                                                             >  1:7640787/1‑137 (MQ=255)
                                                                          gCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGAATATGCGTCCCCTGTAGGA‑G‑C‑GGA‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCt                                                            <  1:7978424/137‑1 (MQ=255)
                                                                                     agggAAAGCCCGCTCCCAAAAAGGATAGGAGTCCCCTGGAGGA‑G‑C‑GGG‑C‑TTAC‑CC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGag                                                 <  1:2462113/136‑1 (MQ=38)
                                                                                         cAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGagcagc                                             <  2:14159835/137‑1 (MQ=255)
                                                                                               cGCTCCTACAAAGGATATGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCt                                       >  2:8174062/1‑137 (MQ=255)
                                                                                                   ccTACAAAGGATATGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGtt                                   >  2:9434407/1‑137 (MQ=255)
                                                                                                   ccTACAAAGGATATGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGtt                                   >  1:11379195/1‑137 (MQ=255)
                                                                                                    cTACAAAGGATATGCGTCCCCTGTAGGA‑G‑C‑GGG‑T‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGttt                                  <  2:4809440/137‑1 (MQ=255)
                                                                                                     tACAAAGGATATGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGtttt                                 >  1:15109326/1‑137 (MQ=255)
                                                                                                       cAAAGGATATGCTTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACTCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTcc                               <  2:2112525/137‑1 (MQ=255)
                                                                                                       cAAAGGATATGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑G‑CCTGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTcc                               <  2:5160788/137‑1 (MQ=255)
                                                                                                       cAAAGGAGATGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTcc                               <  1:1597932/137‑1 (MQ=255)
                                                                                                       aaaaGGATATGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑TCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTcc                               <  1:6248077/136‑1 (MQ=255)
                                                                                                         aaGGATATGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCg                                                                              >  1:873742/1‑88 (MQ=255)
                                                                                                         aaGGATATGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAg                             >  1:15541852/1‑137 (MQ=255)
                                                                                                         aaGGATATGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAg                             >  1:10897748/1‑137 (MQ=255)
                                                                                                         aaGGATATGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAg                             >  2:4019794/1‑137 (MQ=255)
                                                                                                          aGGATATGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGc                            <  1:15355795/137‑1 (MQ=255)
                                                                                                            gATGTGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTc                          <  1:5535378/137‑1 (MQ=255)
                                                                                                              tatGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑NCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGa                        <  1:18744042/137‑1 (MQ=255)
                                                                                                               atGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGat                       >  2:5356100/1‑137 (MQ=255)
                                                                                                               atGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGat                       >  2:11490165/1‑137 (MQ=255)
                                                                                                               atGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGat                       >  2:3019789/1‑137 (MQ=255)
                                                                                                               atGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGat                       >  2:3531021/1‑137 (MQ=255)
                                                                                                               atGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCAAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGat                       >  2:13400077/1‑137 (MQ=255)
                                                                                                                 gCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat                     >  2:12940243/1‑137 (MQ=255)
                                                                                                                 gCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat                     >  1:17194717/1‑137 (MQ=255)
                                                                                                                 gCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat                     >  1:14474947/1‑137 (MQ=255)
                                                                                                                 gCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat                     >  2:2556518/1‑137 (MQ=255)
                                                                                                                 gCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat                     >  1:3530363/1‑137 (MQ=255)
                                                                                                                 gCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat                     >  2:10662726/1‑137 (MQ=255)
                                                                                                                 gCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat                     >  2:9209725/1‑137 (MQ=255)
                                                                                                                 gCGTCCCATGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat                     >  2:17323401/1‑137 (MQ=255)
                                                                                                                   gTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGTAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGc                   >  2:13794962/1‑137 (MQ=255)
                                                                                                                   gTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGc                   <  2:10428018/137‑1 (MQ=255)
                                                                                                                    tCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGc                                                                                    >  1:8373822/1‑71 (MQ=255)
                                                                                                                    tCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGca                  >  1:11628183/1‑137 (MQ=255)
                                                                                                                       ccTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTc                          <  2:5844108/126‑1 (MQ=255)
                                                                                                                       ccTGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcat               <  2:14022756/137‑1 (MQ=255)
                                                                                                                         tGTAGGA‑G‑T‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcatca             >  2:10917889/1‑137 (MQ=255)
                                                                                                                         tGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcatca             >  2:8342798/1‑137 (MQ=255)
                                                                                                                         tGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcatca             >  2:7378171/1‑137 (MQ=255)
                                                                                                                         tGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcatca             >  2:5410502/1‑137 (MQ=255)
                                                                                                                         tGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcatca             >  2:12079380/1‑137 (MQ=255)
                                                                                                                         tGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcatca             >  2:9127958/1‑137 (MQ=255)
                                                                                                                         tGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcatca             >  1:10216384/1‑137 (MQ=255)
                                                                                                                         tGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcatca             >  1:14016980/1‑137 (MQ=255)
                                                                                                                         tGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTAGTTTTCCAGCTCGCTATGCAcatca             >  1:5103643/1‑137 (MQ=255)
                                                                                                                         tGTAGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCACCAAATTGTCGACTACGCGGATCAGCGGTCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCCGCTGCCCTGGTTTTCCAGCTCGATATGCAtatca             >  2:1436489/1‑137 (MQ=39)
                                                                                                                            aGGA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCg          >  1:6919489/1‑137 (MQ=255)
                                                                                                                              gA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGtt        >  2:14876242/1‑137 (MQ=255)
                                                                                                                              gA‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGtt        >  2:7149756/1‑137 (MQ=255)
                                                                                                                               a‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACTTCATCGTTg       <  1:9376868/137‑1 (MQ=255)
                                                                                                                               a‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTg       <  1:10561297/137‑1 (MQ=255)
                                                                                                                               a‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTg       <  2:4433290/137‑1 (MQ=255)
                                                                                                                               a‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTg       <  1:3370677/137‑1 (MQ=255)
                                                                                                                               a‑G‑C‑GGG‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTg       <  2:7796768/137‑1 (MQ=255)
                                                                                                                                 cac‑GGG‑C‑TT‑GCCC‑ACAAAGGGCTGCAAAGAAGACCCAAACTTTCGACCCCGCGGATCAACGATCCGGCAGGGTGATGTTCCGCTCCCAAATCGAGCAGCTGCCCTGGTTTTCCCGCTCGATATGCACACCCTCGttt       >  2:6783153/3‑136 (MQ=255)
                                                                                                                                      gg‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGa   >  2:2365445/1‑137 (MQ=255)
                                                                                                                                       g‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCTGAt  >  2:16133172/1‑137 (MQ=255)
                                                                                                                                       g‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGAt  >  2:8559499/1‑137 (MQ=255)
                                                                                                                                       g‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGAt  >  1:4198386/1‑137 (MQ=255)
                                                                                                                                       g‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGAt  >  2:4248616/1‑137 (MQ=255)
                                                                                                                                       g‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGAt  >  1:1824262/1‑137 (MQ=255)
                                                                                                                                       g‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGAt  >  1:13494565/1‑137 (MQ=255)
                                                                                                                                       g‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGAt  >  2:15995127/1‑137 (MQ=255)
                                                                                                                                       g‑C‑TT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCAAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTCTAGCTCCAGAATCAAGCCGCTACACTGGTTTTCCAGCTAGATATGCACATCATCGCTACCGAt  >  2:18689635/1‑137 (MQ=255)
                                                                                                                                               |                                                                                                                                    
TCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGA‑G‑C‑GGG‑C‑TT‑G‑CC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGAT  >  AE015451_phaM1/1932502‑1932769

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