Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,756,138 Δ1 bp intergenic (+183/‑21) ssrA → / → intA tmRNA/CP4‑57 prophage; integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,756,1370C.88.6% 480.2 / 66.0 123intergenic (+182/‑22)ssrA/intAtmRNA/CP4‑57 prophage; integrase
Reads supporting (aligned to +/- strand):  ref base C (8/6);  new base . (55/54);  total (63/60)
Fisher's exact test for biased strand distribution p-value = 7.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 135 total aligned reads displayed.

TTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACT  >  NC_000913/2756089‑2756185
                                                |                                                
ttAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGc                                                  >  1:1157994/1‑49 (MQ=255)
ttAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGc                                                  >  2:1793842/1‑49 (MQ=255)
ttAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGc                                                  >  2:3475224/1‑49 (MQ=255)
ttAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGc                                                  <  1:1937547/49‑1 (MQ=255)
 tAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGct                                                 >  1:4729146/1‑48 (MQ=255)
 tAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGct                                                 >  2:4909082/1‑48 (MQ=255)
 tAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGct                                                 >  1:669076/1‑48 (MQ=255)
 tAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGct                                                 <  2:1847346/49‑2 (MQ=255)
 tAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGct                                                 >  1:1882005/1‑48 (MQ=255)
  aGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGcta                                                >  2:3265113/1‑47 (MQ=255)
  aGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCt                                                >  2:2559054/1‑49 (MQ=255)
    gCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCTAAAc                                            >  1:5084526/1‑50 (MQ=255)
     cccATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACg                                           >  1:3700899/1‑50 (MQ=255)
     cccATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACg                                           <  1:1969422/50‑1 (MQ=255)
     cccATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACg                                           <  2:232419/50‑1 (MQ=255)
     cccATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACg                                           <  1:421861/50‑1 (MQ=255)
      ccATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACg                                           >  1:179576/1‑49 (MQ=255)
       cATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACg                                           >  1:3303981/1‑48 (MQ=255)
       cATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACg                                           >  2:1887861/1‑48 (MQ=255)
        aTTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACg                                           >  1:4486796/1‑48 (MQ=255)
        aTTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAg                                        <  1:324201/50‑1 (MQ=255)
        aTTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAg                                        >  1:2432360/1‑50 (MQ=255)
        aTTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAg                                        >  2:3863476/1‑50 (MQ=255)
        aTTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAg                                        <  2:1562377/50‑1 (MQ=255)
        aTTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAg                                        <  2:46059/50‑1 (MQ=255)
         ttGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAgg                                       >  2:1021486/1‑50 (MQ=255)
         ttGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAgg                                       <  2:4089556/50‑1 (MQ=255)
         ttGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAgg                                       <  1:5182857/50‑1 (MQ=255)
          tGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAgg                                       >  2:798845/1‑49 (MQ=255)
          tGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAgg                                       >  2:2431220/1‑49 (MQ=255)
           gATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAggg                                      <  1:3638328/49‑1 (MQ=255)
            aTAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAgg                                       >  1:1864727/1‑47 (MQ=255)
            aTAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGga                                     <  2:129805/49‑1 (MQ=255)
              aGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGAc                                  >  1:1579300/1‑50 (MQ=255)
               ggCCCAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGAc                                  <  1:3829087/49‑1 (MQ=255)
                gCCCAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTc                                <  1:1115441/50‑1 (MQ=255)
                  ccAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTccc                              <  2:735958/50‑1 (MQ=255)
                  ccAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTccc                              <  2:2941472/50‑1 (MQ=255)
                  ccAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTccc                              >  1:2470093/1‑50 (MQ=255)
                   cAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTccc                              >  1:483782/1‑50 (MQ=255)
                   cAACGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTccc                              >  1:1849723/1‑49 (MQ=255)
                    aaCGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCA                            >  2:860616‑M1/1‑50 (MQ=255)
                    aaCGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCA                            <  1:1047864‑M1/50‑1 (MQ=255)
                    aaCGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCA                            >  2:4536377‑M1/1‑50 (MQ=255)
                    aaCGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCA                            <  1:14458‑M1/50‑1 (MQ=255)
                    aaCGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCA                            >  2:3395196‑M1/1‑50 (MQ=255)
                     aCGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAc                           <  2:4196361‑M1/50‑2 (MQ=255)
                      cGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcc                          >  1:2094557‑M1/1‑48 (MQ=255)
                      cGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcc                          <  2:1561178‑M1/50‑3 (MQ=255)
                      cGAAAAGCTCTATTGTTTACGTCGGGCTAAACGCAGGGAGACTCCCCAcc                          >  2:3423376‑M1/1‑48 (MQ=255)
                       gAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcct                         >  1:4072192‑M1/1‑47 (MQ=255)
                       gAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcct                         >  1:2969949‑M1/1‑47 (MQ=255)
                       gAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcct                         <  1:2760853‑M1/50‑4 (MQ=255)
                        aaaaGTTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctt                        <  1:4841930‑M1/50‑5 (MQ=255)
                        aaaaGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctt                        <  1:897800‑M1/50‑5 (MQ=255)
                        aaaaGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctt                        <  2:4187221‑M1/50‑5 (MQ=255)
                        aaaaGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctt                        >  1:478730‑M1/1‑46 (MQ=255)
                         aaaGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAccttt                       <  1:679807‑M1/50‑6 (MQ=255)
                         aaaGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAccttt                       >  1:5072478‑M1/1‑45 (MQ=255)
                          aaGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAccttt                       >  2:4598057‑M1/1‑44 (MQ=255)
                           aGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCaa                      <  1:4871058/50‑1 (MQ=255)
                           aGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttga                     <  1:2899495‑M1/50‑8 (MQ=255)
                            gCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaa                    <  1:1243266‑M1/50‑9 (MQ=255)
                            gCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttga                     >  1:4736270‑M1/1‑42 (MQ=255)
                            gCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttga                     >  1:4663328‑M1/1‑42 (MQ=255)
                             ctctATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaa                   <  2:4729923‑M1/50‑10 (MQ=255)
                             ctctATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaa                    >  1:4459740‑M1/1‑41 (MQ=255)
                              tctATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGa                    >  2:4939522/1‑49 (MQ=255)
                              tctATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaa                  >  2:3478096‑M1/1‑40 (MQ=255)
                                tATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaa                  >  1:2104374‑M1/1‑38 (MQ=255)
                                tATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaa                  >  1:1734357‑M1/1‑38 (MQ=255)
                                  ttGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGaaa                  <  2:842642/47‑1 (MQ=255)
                                  ttGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacag               <  2:2157080‑M1/49‑14 (MQ=255)
                                  ttGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacag               >  2:581703‑M1/1‑36 (MQ=255)
                                  ttGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaaca                <  2:2134139‑M1/48‑13 (MQ=255)
                                  ttGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaaca                <  2:2754093‑M1/48‑13 (MQ=255)
                                    gTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggat            >  1:2624056‑M1/1‑34 (MQ=255)
                                     tttACGTTGGGCTAAACGCAGGGAGACTCCCCTcctttgaaaacaggat            <  2:4354265‑M1/49‑17 (MQ=255)
                                     tttACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatg           <  1:4032939‑M1/50‑18 (MQ=255)
                                     tttACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatg           >  1:4264143‑M1/1‑33 (MQ=255)
                                     tttACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatg           >  1:2901833‑M1/1‑33 (MQ=255)
                                     tttACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatg           >  1:1670943‑M1/1‑33 (MQ=255)
                                     tttACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggat            <  2:150967‑M1/49‑17 (MQ=255)
                                     tttACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggat            <  1:3159643‑M1/49‑17 (MQ=255)
                                     tttACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggac            <  2:2876345‑M1/49‑17 (MQ=255)
                                       tACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgt          <  2:3798900‑M1/49‑19 (MQ=255)
                                       tACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatg           >  1:2495381‑M1/1‑31 (MQ=255)
                                        aCGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgta         <  1:89477‑M1/49‑20 (MQ=255)
                                         cGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagc       >  1:1639938‑M1/1‑29 (MQ=255)
                                         cGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagc       <  2:287554‑M1/50‑22 (MQ=255)
                                         cGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagc       >  1:1394479‑M1/1‑29 (MQ=255)
                                         cGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagc       <  1:918749‑M1/50‑22 (MQ=255)
                                         cGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagc       <  2:507646‑M1/50‑22 (MQ=255)
                                         cGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtag        <  1:1507130‑M1/49‑21 (MQ=255)
                                          gTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagc       >  2:5164411‑M1/1‑28 (MQ=255)
                                          gTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagc       >  2:5071393‑M1/1‑28 (MQ=255)
                                           ttGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcga     >  1:360646‑M1/1‑27 (MQ=255)
                                           ttGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcga     <  1:690882‑M1/50‑24 (MQ=255)
                                               gCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACt  <  1:2736790/50‑1 (MQ=255)
                                               gCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACt  >  1:1598521/1‑50 (MQ=255)
                                                |                                                
TTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACT  >  NC_000913/2756089‑2756185

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