Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 282,361 C→T G28G (GGC→GGT yagE → 2‑keto‑3‑deoxy gluconate (KDG) aldolase; CP4‑6 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913282,3610CT98.6% 274.1 / ‑5.4 71G28G (GGC→GGTyagE2‑keto‑3‑deoxy gluconate (KDG) aldolase; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (42/28);  minor base G (1/0);  total (43/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCG  >  NC_000913/282214‑282503
                                                                                                                                                   |                                                                                                                                              
tGATCGTTAAATAAACGCACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTAc                                                                                                                                              >  1:139525/1‑150 (MQ=255)
    cGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTAccgccg                                                                                                                                         >  2:2063/1‑151 (MQ=255)
        aaataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCt                                                                                                                                     >  1:961829/1‑151 (MQ=255)
         aataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTg                                                                                                                                    >  1:210875/1‑151 (MQ=255)
         aataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTg                                                                                                                                    <  2:110104/151‑1 (MQ=255)
          ataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGa                                                                                                                                   <  2:1039073/151‑1 (MQ=255)
             aacgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATcg                                                                                                                                >  2:462426/1‑151 (MQ=255)
               cgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATcgac                                                                                                                              >  2:122879/1‑151 (MQ=255)
                 aacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATcgacg                                                                                                                             >  2:282323/1‑150 (MQ=255)
                       gTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTg                                                                                                                        <  1:479082/149‑1 (MQ=255)
                       gTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGAt                                                                                                                      <  2:598001/151‑1 (MQ=255)
                           tataATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCaa                                                                                                                   >  1:393595/1‑150 (MQ=255)
                              aaTGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCa                                                                                                               <  1:1042627/151‑1 (MQ=255)
                              aaTGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCa                                                                                                               <  1:891144/151‑1 (MQ=255)
                                  tAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCg                                                                                                           >  1:707815/1‑151 (MQ=255)
                                   aGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCg                                                                                                           >  2:349723/1‑150 (MQ=255)
                                   aGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGCCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCCCCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCg                                                                                                           <  2:914564/150‑1 (MQ=255)
                                     aaCAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTg                                                                                                        <  1:392447/151‑1 (MQ=255)
                                     aaCAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTg                                                                                                        <  1:676378/151‑1 (MQ=255)
                                       cAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGAc                                                                                                      >  1:692649/1‑151 (MQ=255)
                                       cAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGAc                                                                                                      >  1:71503/1‑151 (MQ=255)
                                        aaaaTGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGAc                                                                                                      <  2:847895/150‑1 (MQ=255)
                                        aaaaTGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACg                                                                                                     >  2:325118/1‑151 (MQ=255)
                                         aaaTGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACgg                                                                                                    >  1:108368/1‑151 (MQ=255)
                                          aaTGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGc                                                                                                   >  1:886359/1‑151 (MQ=255)
                                           aTGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGcc                                                                                                  >  2:313191/1‑151 (MQ=255)
                                           aTGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGcc                                                                                                  >  1:953621/1‑151 (MQ=255)
                                           aTGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGcc                                                                                                  >  1:85159/1‑151 (MQ=255)
                                               ttCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGtt                                                                                              >  2:48666/1‑151 (MQ=255)
                                               ttCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGtt                                                                                              >  2:185076/1‑151 (MQ=255)
                                                    cAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTcc                                                                                         >  1:159915/1‑151 (MQ=255)
                                                    cAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTcc                                                                                         >  2:920872/1‑151 (MQ=255)
                                                     agggATATCTCAGGCCGCAGTCCGCTTTGTTCACGGGAATCATTCCCCCGGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCt                                                                                        <  2:1011253/149‑1 (MQ=255)
                                                          gaTCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGc                                                                                    >  2:486100/1‑150 (MQ=255)
                                                            tctcATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGc                                                                                 >  1:389331/1‑151 (MQ=255)
                                                             ctcATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGc                                                                                    <  2:990129/147‑1 (MQ=255)
                                                              tcATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCg                                                                                >  2:500455/1‑150 (MQ=255)
                                                               cATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGt                                                                              >  1:435576/1‑151 (MQ=255)
                                                                       aGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGt                                                                       >  2:658349/1‑150 (MQ=255)
                                                                         tCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTtc                                                                     <  2:716960/150‑1 (MQ=255)
                                                                             cgTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGAGGTGGCGAGTTCTccc                                                                 >  1:330940/1‑150 (MQ=255)
                                                                                 gttCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCt                                                             <  2:177284/150‑1 (MQ=255)
                                                                                 gttCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCt                                                             <  1:546933/150‑1 (MQ=255)
                                                                                      cGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGgcgc                                                       >  2:227550/1‑151 (MQ=255)
                                                                                      cGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGgcgc                                                       <  1:26434/151‑1 (MQ=255)
                                                                                        ggAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCg                                                     >  2:648725/1‑151 (MQ=255)
                                                                                           aTCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGa                                                    >  2:515135/1‑149 (MQ=255)
                                                                                                 cccccTGTCTCCCCCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTa                                            <  2:954517/151‑1 (MQ=255)
                                                                                                  ccccTGTCTCCCCCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTaa                                           <  2:608695/151‑1 (MQ=255)
                                                                                                      tGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGcc                                       >  1:639993/1‑151 (MQ=255)
                                                                                                      tGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGcc                                       >  2:123115/1‑151 (MQ=255)
                                                                                                      tGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGcc                                       >  2:831871/1‑151 (MQ=255)
                                                                                                        tctcCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCg                                              >  2:67586/1‑142 (MQ=255)
                                                                                                         ctcCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCg                                              <  1:67587/141‑1 (MQ=255)
                                                                                                         ctcCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCAtt                                    >  2:159808/1‑151 (MQ=255)
                                                                                                           cccccaTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTg                                   <  2:761976/147‑1 (MQ=255)
                                                                                                           ccaccaTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCAtt                                    <  2:393585/149‑1 (MQ=255)
                                                                                                                 tttttACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCtt                            >  2:833640/1‑151 (MQ=255)
                                                                                                                 tttttACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCt                             <  2:335583/150‑1 (MQ=255)
                                                                                                                   tttACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTg                          >  1:1007498/1‑151 (MQ=255)
                                                                                                                      aCCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTa                       <  1:654308/151‑1 (MQ=255)
                                                                                                                        cgccgACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATc                     <  1:103037/151‑1 (MQ=255)
                                                                                                                            gACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATc                 >  2:724247/1‑151 (MQ=255)
                                                                                                                            gACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATc                 >  2:193152/1‑151 (MQ=255)
                                                                                                                             aCGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCa                <  1:123119/151‑1 (MQ=255)
                                                                                                                                 ccAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGt             <  2:35451/150‑1 (MQ=255)
                                                                                                                                 ccAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTc            <  1:220085/151‑1 (MQ=255)
                                                                                                                                      tccaTAAACCGGGGACCCCCCCGCTGATCGACGATCTGATCAAAGCAAGCGGTGGCGGCCCGTTCTTCCTGGGCAGGGGGGGCGAGTTCTCCCCGCTCGGCGCCGAAGAGCGGAAAGCCATTTCCCGCTTTTCTATCGATCCTGTCGAtc        >  2:882874/4‑150 (MQ=255)
                                                                                                                                         aTAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGAtcgt      <  2:754085/149‑1 (MQ=255)
                                                                                                                                           aaGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTcgcg  >  2:743670/1‑151 (MQ=255)
                                                                                                                                           aaGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTcgcg  >  2:964471/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
TGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCG  >  NC_000913/282214‑282503

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: