Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,768,100 C→T A699A (GCC→GCT)  PP_1577 → lambda family tail tape meausure protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,768,1000CT97.9% 142.7 / NA 48A699A (GCC→GCT) PP_1577lambda family tail tape meausure protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (21/26);  minor base G (0/1);  total (21/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGCGGCACAGGTGGCCATGGCCGGCCTGAACGCATTCACCGCGATGATCGGGGTTCCCTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGC  >  NC_002947/1767957‑1768218
                                                                                                                                               |                                                                                                                      
gcggcgGCACAGGTGGCCATGGCCGGCCTGAACGCATTCACCGCGATGATCGGGGTTCCCTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCt                                                                                                                  <  1:190787/150‑1 (MQ=255)
            gTGGCCATGGCCGGCCTGAACGCATTCACCGCGATGATCGGGGTTCCCTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCAc                                                                                                      >  1:944446/1‑150 (MQ=255)
               gCCATGGCCGGCCTGAACGCATTCACCGCGATGATCGGGGTTCCCTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACggt                                                                                                   >  2:1450091/1‑150 (MQ=255)
                 cATGGCCGGCCTGAACGCATTCACCGCGATGATCGGGGTTCCCTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACggtgg                                                                                                 <  1:207085/150‑1 (MQ=255)
                  aTGGCCGGCCTGAACGCATTCACCGCGATGATCGGGGTTCCCTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACggtggt                                                                                                >  1:311696/1‑150 (MQ=255)
                    ggccggccTGAACGCATTCACCGCGATGATCGGGGTTCCCTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACggtctcag                                                                                              <  2:50843/150‑6 (MQ=255)
                          ccTGAACGCATTCACCGCGATGATCGGGGTTCCCTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGaaaa                                                                                        <  2:311696/150‑1 (MQ=255)
                              aaCGCATTCACCGCGATGATCGGGGTTCCCTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCttt                                                                                                                    >  2:369662/1‑118 (MQ=255)
                              aaCGCATTCACCGCGATGATCGGGGTTCCCTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCttt                                                                                                                    <  1:369662/118‑1 (MQ=255)
                                cGCATTCACCGCGATGATCGGGGTTCCCTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATccc                                                                                  >  2:285414/1‑150 (MQ=255)
                                cGCATTCACCGCGATGATCGGGGTTCCCTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGAAGTCCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATccc                                                                                  >  2:1238165/1‑150 (MQ=255)
                                    ttCGCCGCGATGGTCGGGTTTCGCTGCGTGGGGCCGCCAACTGGTCCTGTCGCTGCCATTGCGGCGGCAGCCGCAACAGGCGTGTTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCga                                                                              <  2:577942/150‑1 (MQ=255)
                                                     ggTTCCCTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGagag                                                                           >  2:1119275/1‑136 (MQ=255)
                                                     ggTTCCCTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGagag                                                                           <  1:1119275/136‑1 (MQ=255)
                                                     ggTTCCCTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGgctgct                                                             <  1:285414/150‑1 (MQ=255)
                                                        tgcggtCGTGGGGCGGGCAAGTGCAGCGGGCGCTGCGATTGCAGCGGCAGCCGCGGCAGGCGTGCTGAGGACGGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGa                                                          <  2:534620/145‑1 (MQ=255)
                                                         cccTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGAt                                                         <  1:701276/150‑1 (MQ=255)
                                                           cTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATgg                                                       >  2:1237580/1‑150 (MQ=255)
                                                                 ggggCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGgcga                                                 <  1:1011915/150‑1 (MQ=255)
                                                                  gggCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGgcgag                                                >  1:798720/1‑150 (MQ=255)
                                                                   ggCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGgcgagc                                               <  2:798720/150‑1 (MQ=255)
                                                                     ccGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGgcgagcga                                             >  2:470620/1‑150 (MQ=255)
                                                                     ccGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGCCGGGCCGATGACTGGGGTTGACTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGAGTACCTGGCTGCTGGATGGCGgcgagcga                                             <  1:1218234/150‑1 (MQ=255)
                                                                            ttGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTcc                                         <  1:240883/147‑1 (MQ=255)
                                                                            ttGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTcc                                         >  2:240883/1‑147 (MQ=255)
                                                                             tGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAg                                     <  1:228689/150‑1 (MQ=255)
                                                                              gCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGc                                    <  1:470620/150‑1 (MQ=255)
                                                                               cACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGgcgagc                                               >  1:1010685/1‑138 (MQ=255)
                                                                               cACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGgcgagc                                               <  2:1010685/138‑1 (MQ=255)
                                                                                  cTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTg                                >  1:387297/1‑150 (MQ=255)
                                                                                      cgcTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCTGCGAGAggg                                                                         >  1:1365296/1‑105 (MQ=255)
                                                                                      cgcTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAggg                                                                         <  2:1365296/105‑1 (MQ=255)
                                                                                       gcTGCCATTGCAGCCACAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAg                           >  1:1572147/1‑150 (MQ=255)
                                                                                             aTTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCgg                     <  1:569231/150‑1 (MQ=255)
                                                                                             aTTGCAGCAGCAGACGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCgg                     >  2:468897/1‑150 (MQ=255)
                                                                                                gcagcggcagccGGAGCAGGCGTGCGGAGGACTGGGGTGGGCTCGGCGCTTTCTGGGCAGGCTCACGGTGGTATGGAAAGCATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGAc                  <  1:513659/150‑1 (MQ=255)
                                                                                                 cagccgcagccgCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCCGGCCTTCAGCCCTGCGCAGAACCGGGAcc                 >  1:402537/1‑150 (MQ=255)
                                                                                                   gccgcagccgcaGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTc               <  1:429275/150‑1 (MQ=255)
                                                                                                     cgcagccgcagcAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCAc             <  2:1555762/150‑1 (MQ=255)
                                                                                                      gcagccgcagcAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCAcc            >  2:850527/1‑150 (MQ=255)
                                                                                                       cagccgcagcAGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCg           <  2:944446/150‑1 (MQ=255)
                                                                                                            gcagcaGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTAc      <  1:850527/150‑1 (MQ=255)
                                                                                                            gcagcaGGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTAc      <  1:468897/150‑1 (MQ=255)
                                                                                                                cagGCGTGCTGATGACTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGc  >  1:488516/1‑150 (MQ=255)
                                                                                                                         tgatgaCTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGAc                  >  2:542327/1‑125 (MQ=255)
                                                                                                                         tgatgaCTGGGGTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGAc                  <  1:542327/125‑1 (MQ=255)
                                                                                                                                    gTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGctgctg                                                            <  1:924009/72‑1 (MQ=255)
                                                                                                                                    gTTGGCTCGGCTCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGctgctg                                                            >  2:924009/1‑72 (MQ=255)
                                                                                                                                               |                                                                                                                      
GCGGCGGCACAGGTGGCCATGGCCGGCCTGAACGCATTCACCGCGATGATCGGGGTTCCCTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGC  >  NC_002947/1767957‑1768218

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