Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 288,667 A→G intergenic (‑8/+264) yagK ← / ← yagL inovirus Gp2 family protein/recombinase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273288,6670AG93.3% 195.8 / 4.2 60intergenic (‑8/+264)yagK/yagLinovirus Gp2 family protein/recombinase family protein
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base G (21/35);  total (22/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.54e-01

GTGGAGATCGACACGTACAGCAATTAACCTAGGATACCGGTTAGTCAAATAACATAATAATTTATTAATCTTACTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATACCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCACTGTTGTCGGTATATACTGGTATCAGATTTTAGTACGGACATACTGATCTGATTACTACTCTTA  >  NZ_CP009273/288537‑288801
                                                                                                                                  |                                                                                                                                      
gTGGAGATCGACACGTACAGCAATTAACCTAGGATACCGGTTAGTCAAATAACATAATAATTTATTAATCTTACTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTAtt                                                                                                                                >  2:560623/1‑139 (MQ=255)
   gagaTCGACACGTACAGCAATTAACCTAGGATACCGGTTAGTCAAATAACATAATAATTTATTAATCTTACTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTg                                                                                                                              <  1:560623/138‑1 (MQ=255)
      aTCGACACGTACAGCAATTAACCTAGGATACCGGTTAGTCAAATAACATAATAATTTATTAATCTTACTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTcc                                                                                                                          >  2:430915/1‑139 (MQ=255)
            acGTACAGCAATTAACCTAGGATACCGGTTAGTCAAATAACATAATAATTTATTAATCTTACTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGt                                                                                                                    <  1:1191052/139‑1 (MQ=255)
               tACAGCAATTAACCTAGGATACCGGTTAGTCAAATAACATAATAATTTATTAATCTTACTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCta                                                                                                                 >  2:1179750/1‑139 (MQ=255)
                    cAATTAACCTAGGATACCGGTTAGTCAAATAACATAATAATTTATTAATCTTACTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCa                                                                                                            <  2:184959/139‑1 (MQ=255)
                    cAATTAACCTAGGATACCGGTTAGTCAAATAACATAATAATTTATTAATCTTACTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCa                                                                                                            <  1:846853/139‑1 (MQ=255)
                    cAATTAACCTAGGATACCGGTTAGTCAAATAACATAATAATTTATTAATCTTACTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCa                                                                                                            <  2:146651/139‑1 (MQ=255)
                      aTTAACCTAGGATACCGGTTAGTCAAATAACATAATAATTTATTAATCTTACTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGa                                                                                                          <  1:225567/139‑1 (MQ=255)
                      aTTAACCTAGGATACCGGTTAGTCAAATAACATAATAATTTATTAATCTTACTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGa                                                                                                          <  1:1100776/139‑1 (MQ=255)
                      aTTAACCTAGGATACCGGTTAGTCAAATAACATAATAATTTATTAATCTTACTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGa                                                                                                          <  2:1073020/139‑1 (MQ=255)
                           ccTAGGATACCGGTTAGTCAAATAACATAATAATTTATTAATCTTACTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTa                                                                                                     <  1:1179750/139‑1 (MQ=255)
                               ggATACCGGTTAGTCAAATAACATAATAATTTATTAATCTTACTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAg                                                                                                  <  2:1116355/138‑1 (MQ=255)
                                     cGGTTAGTCAAATAACATAATAATTTATTAATCTTACTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGt                                                                                                                    >  2:830697/1‑114 (MQ=255)
                                     cGGTTAGTCAAATAACATAATAATTTATTAATCTTACTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGt                                                                                                                    <  1:830697/114‑1 (MQ=255)
                                            tCAAATAACATAATAATTTATTAATCTTACTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATAcagg                                                                                    <  1:1260151/139‑1 (MQ=255)
                                                 tAACATAATAATTTATTAATCTTACTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACg                                                                                                                      >  1:118467/1‑100 (MQ=255)
                                                 tAACATAATAATTTATTAATCTTACTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACg                                                                                                                      <  2:118467/100‑1 (MQ=255)
                                                   aCATAATAATTTATTAATCTTACTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGc                                                                              <  1:1098418/138‑1 (MQ=255)
                                                   aCATAATAATTTATTAATCTTACTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGc                                                                              <  1:689189/138‑1 (MQ=255)
                                                              tattaATCTTACTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCa                                                                  >  2:955394/1‑139 (MQ=255)
                                                               attaATCTTACTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCa                                                                  <  2:241866/138‑1 (MQ=255)
                                                               attaATCTTACTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATACCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCa                                                                  <  1:1034220/138‑1 (MQ=255)
                                                                        aCTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATACCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCACtgttgt                                                           <  1:553152/136‑1 (MQ=255)
                                                                        aCTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTGTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATAGAGGCAGGGCTATTCCTGGTCACTGTTGTCg                                                         <  1:1005532/138‑1 (MQ=255)
                                                                         ctctGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCACTGTTGTCGGt                                                       <  2:357679/139‑1 (MQ=255)
                                                                              atatTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATa                                                                                        <  1:1273556/101‑1 (MQ=255)
                                                                              atatTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATa                                                                                        >  2:1273556/1‑101 (MQ=255)
                                                                              atatTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCACTGTTGTCGGtatata                                                  >  1:1155819/1‑139 (MQ=255)
                                                                                  tCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCACTGTTGTCGGTATATACTgg                                              >  1:404048/1‑139 (MQ=255)
                                                                                  tCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATACCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATAcaggc                                                                                   >  1:1145035/1‑102 (MQ=255)
                                                                                  tCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATACCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATAcaggc                                                                                   <  2:1145035/102‑1 (MQ=255)
                                                                                    gAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCACTGTTGTc                                                          >  2:141626/1‑125 (MQ=255)
                                                                                    gAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCACTGTTGTc                                                          <  1:141626/125‑1 (MQ=255)
                                                                                               cTCCATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCACTGTTGTCGGTATATACTGGTATCAGATTTTAg                                 <  2:1155819/139‑1 (MQ=255)
                                                                                               cTCCATGAGTACCTTGATATGTTCTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCACTGTTGTCGGTATATACTGGTATCAGATTTTAg                                 <  2:1248015/139‑1 (MQ=255)
                                                                                                  cATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCACtgttgt                                                           <  1:686817/110‑1 (MQ=255)
                                                                                                  cATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCACTGTTGTCGGTATATACTGGTATCAGATTTTAGTAc                              >  1:476639/1‑139 (MQ=255)
                                                                                                  cATGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGAATGCCATACAGGCAGGGCTATTCCTGGTCACtgttgt                                                           >  2:686817/1‑110 (MQ=255)
                                                                                                    tGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTa                                                                                                                        >  2:1169134/1‑47 (MQ=255)
                                                                                                    tGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTa                                                                                                                        <  1:1169134/47‑1 (MQ=255)
                                                                                                    tGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAggg                                                                                                <  1:708592/71‑1 (MQ=255)
                                                                                                    tGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAggg                                                                                                >  2:708592/1‑71 (MQ=255)
                                                                                                    tGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGc                                                                              >  1:846451/1‑89 (MQ=255)
                                                                                                    tGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGc                                                                              <  2:846451/89‑1 (MQ=255)
                                                                                                    tGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCACTGTTGTCGGTATATACTGGTATCAGATTTTAg                                 <  1:514970/134‑1 (MQ=255)
                                                                                                    tGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCACTGTTGTCGGTATATACTGGTATCAGATTTTAg                                 >  2:514970/1‑134 (MQ=255)
                                                                                                    tGAGTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCACTGTTGTCGGTATATACTGGTATCAGATTTTAGTACg                             >  2:1120961/1‑138 (MQ=255)
                                                                                                       gTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGa                                                                                                          >  1:73229/1‑58 (MQ=255)
                                                                                                       gTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGa                                                                                                          <  2:73229/58‑1 (MQ=255)
                                                                                                       gTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCACtgttgt                                                           >  1:945159/1‑105 (MQ=255)
                                                                                                       gTACCTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCACtgttgt                                                           <  2:945159/105‑1 (MQ=255)
                                                                                                          ccTTGATATGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCACTGTTGTCGGTATATACTGGTATCAGATTTTAGTACGGACATa                       <  2:476639/138‑1 (MQ=255)
                                                                                                               atatGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCACTGTTGTCGGTATATACTGGTATCAGATTTTAGTACGGACATActgat                  <  1:430915/138‑1 (MQ=255)
                                                                                                                tatGTTTTCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAAGGCTATTCCTGGTCACTGTTGTCGGTATATACTGGTATCAGATTTTAGTCCGGACATGctgatct                <  2:403464/139‑1 (MQ=255)
                                                                                                                    ttttCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCa                                                                  >  1:389676/1‑85 (MQ=255)
                                                                                                                    ttttCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCa                                                                  <  2:389676/85‑1 (MQ=255)
                                                                                                                    ttttCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCACTGTTGTCGGTATATACTGGTATCAGATTTTAGTACgg                            <  2:1210561/123‑1 (MQ=255)
                                                                                                                    ttttCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCACTGTTGTCGGTATATACTGGTATCAGATTTTAGTACgg                            >  1:1210561/1‑123 (MQ=255)
                                                                                                                      ttCATAATATATGCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCACTGTTGTCGGTATATACTGGTATCAGATTTTAGTACGGACATACTGATCTGATta           <  2:226388/138‑1 (MQ=255)
                                                                                                                              atataCCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCACTGTTGTCGGTATATACTGGTATCAGATTTTAGTACGGACATACTGATCTGATTACTACTCTTa  >  1:1142791/1‑139 (MQ=255)
                                                                                                                                  |                                                                                                                                      
GTGGAGATCGACACGTACAGCAATTAACCTAGGATACCGGTTAGTCAAATAACATAATAATTTATTAATCTTACTCTGATATTCGAGGATATGTACTCCATGAGTACCTTGATATGTTTTCATAATATATACCTGTATTTGGTCCTACGGTCTATAGCAGACCGTACAGGGATGCCATACAGGCAGGGCTATTCCTGGTCACTGTTGTCGGTATATACTGGTATCAGATTTTAGTACGGACATACTGATCTGATTACTACTCTTA  >  NZ_CP009273/288537‑288801

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