Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BX571857 1,222,841 T→A 20.3% I360N (ATT→AAT)  SAS1146 → primosomal protein n'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BX5718571,222,8410TA20.3% 185.8 / 44.5 64I360N (ATT→AAT) SAS1146primosomal protein n'
Reads supporting (aligned to +/- strand):  ref base T (19/32);  new base A (5/8);  total (24/40)
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.24e-01

AATGGTTTTACGCTTCAAACGTCGATTTGGTGATGACGTTGCTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCAAAAATTTAGGGTTAATCATCATTGATGAAGAACATG  >  BX571857/1222746‑1222936
                                                                                               |                                                                                               
aaTGGTTTTACGCTTCAAACGTCGATTTGGTGATGACGTTGCTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAAtt                                                                                                <  1:577803/97‑1 (MQ=255)
  tGGTTTTACGCTTCAAACGTCGATTTGGTGATGACGTTGCTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAAATAg                                                                                              >  2:889350/1‑97 (MQ=255)
      tttACGCTTCAAACGTCGATTTGGTGATGACGTTGCTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGAt                                                                                          <  1:1061249/97‑1 (MQ=255)
       ttACGCTTCAAACGTCGATTTGGTGATGACGTTGCTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATg                                                                                         <  1:1282925/97‑1 (MQ=255)
       ttACGCTTCAAACGTCGATTTGGTGATGACGTTGCTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAAAGAGGGATg                                                                                         >  1:508453/1‑97 (MQ=255)
         aCGCTTCAAACGTCGATTTGGTGATGACGTTGCTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGt                                                                                       <  1:668362/97‑1 (MQ=255)
          cGCTTCAAACGTCGATGTGGGGATGACGTTGCGGGATTACAGTCTGGCTTGTCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTc                                                                                      <  2:541103/97‑1 (MQ=255)
              tgaaaCGTCGATTTGGTGATGACGTTGCTGTATTACATTCTGGCTTAGCTAATGGGGAACGTTATGATGAGTGGCAAAAAAATAGGGATGGTCGTGc                                                                                  <  2:1033581/95‑1 (MQ=255)
                 aaCGTCGATTTGGTGATGACGTTGCTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCgag                                                                               <  1:32442/97‑1 (MQ=255)
                 aaCGTCGATTTGGTGATGACGTTGCTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGAAAAAAAAAAAAAAATTAGGGATGGTCGTGCgag                                                                               <  2:347064/97‑1 (MQ=255)
                     tCGATTTGGTGATGACGTTGCTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAAATAGGGATGGTCGTGCGAGAGTa                                                                           <  2:437707/97‑1 (MQ=255)
                       gATTTGGTGATGACGTTGCTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAg                                                                         >  1:657637/1‑97 (MQ=255)
                       gATTTGGTGATGACGTTGCTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAg                                                                         >  1:1084919/1‑97 (MQ=255)
                          ttGGTGATGACGTTGCTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAgtgt                                                                      >  1:650005/1‑97 (MQ=255)
                          ttGGTGATGACGTTGCTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAgtgt                                                                      >  1:206436/1‑97 (MQ=255)
                           tGGTGATGACGTTGCTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGtt                                                                     <  1:1211459/97‑1 (MQ=255)
                            ggTGATGACGTTGCTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTg                                                                    <  1:1377859/97‑1 (MQ=255)
                              tgatgaCGTTGCTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGt                                                                  >  2:338073/1‑97 (MQ=255)
                                 tgaCGTTGCTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCa                                                               <  2:689151/97‑1 (MQ=255)
                                   aCGTTGCTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAg                                                             <  2:894835/97‑1 (MQ=255)
                                   aCGTTGCTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAg                                                             <  2:423839/97‑1 (MQ=255)
                                   aCGTTGCTGTATTACATTCTGGCCTATCTACTGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGGCGTGAGAGAGTAAGTGTTGGTGcgcg                                                             >  2:865172/1‑94 (MQ=255)
                                    cGTTGCTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAgg                                                            >  1:1035691/1‑97 (MQ=255)
                                     gTTGCTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCCAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAgga                                                           >  2:621073/1‑96 (MQ=255)
                                     gTTGCTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAAATAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGt                                                           >  2:463724/1‑97 (MQ=255)
                                         cTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAg                                                       <  2:65740/97‑1 (MQ=255)
                                         cTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAAATAGGGATGGTCGTGCGTGAGTAAGTGTTGGTGCAAGGTCAAg                                                       <  2:258547/97‑1 (MQ=255)
                                           gTATTACATTCTGGCTTATCTAATGGGGAGCGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAgtg                                                     <  1:646264/97‑1 (MQ=255)
                                           gTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAgtg                                                     <  1:1008942/97‑1 (MQ=255)
                                           gTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAgtg                                                     >  1:1229299/1‑97 (MQ=255)
                                             aTTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAgtgtg                                                   <  1:778966/97‑1 (MQ=255)
                                             aTTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAgtgtg                                                   >  2:153823/1‑97 (MQ=255)
                                              ttACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAgtgtgt                                                  >  2:231398/1‑97 (MQ=255)
                                                    tCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCAc                                            <  1:849312/97‑1 (MQ=255)
                                                     cTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCAcc                                           <  1:106792/97‑1 (MQ=255)
                                                     cTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAAATAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCAcc                                           <  2:257795/97‑1 (MQ=255)
                                                      tGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCt                                          <  1:1406387/97‑1 (MQ=255)
                                                       ggCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCtt                                         <  1:1022556/97‑1 (MQ=255)
                                                         cTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTc                                       <  1:792613/97‑1 (MQ=255)
                                                             tCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCaaaa                                   <  2:54368/97‑1 (MQ=255)
                                                             tCTAATGGGGAACGTTATGATGAGTGGCAAAAAAATAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCaaaa                                   >  1:199088/1‑97 (MQ=255)
                                                             tCTAATGGGGAACGTTATGATGAGTGGCAAAAAAATAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCaaaa                                   <  2:1398806/97‑1 (MQ=255)
                                                              cTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCaaaaa                                  <  1:1257122/97‑1 (MQ=255)
                                                                 aTGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCAAAAAttt                               <  1:1317856/97‑1 (MQ=255)
                                                                  tGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCAAAAATTTa                              >  1:70044/1‑97 (MQ=255)
                                                                   ggggAACGTTATGATGAGTGTCAAAAAATTAGGGATGGTCGTGTGAGGGTAAGTGTTGGTGCAAGGTCAAGTGTGGTCGCACCGTGCAAAAATTCAg                             >  1:818207/1‑97 (MQ=255)
                                                                   ggggAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCAAAAATTTAg                             >  1:908686/1‑97 (MQ=255)
                                                                   ggggAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCAAAAATTTAg                             <  2:219924/97‑1 (MQ=255)
                                                                   ggggAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCAAAAATTTAg                             <  2:219908/97‑1 (MQ=255)
                                                                    gggAACGTTATGATGAGTGGCAAAAAAATAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCAAAAATTTAgg                            <  2:295358/97‑1 (MQ=255)
                                                                           ttATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCAAAAATTTAGGGTTAatc                     <  2:611239/97‑1 (MQ=255)
                                                                             atgatgAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCAAAAATTTAGGGTTAatcat                   <  2:1045424/97‑1 (MQ=255)
                                                                               tatgaGTGGCAAAAAAATAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCAAAAATTTAGGGTTAAtcatca                 <  1:992860/96‑1 (MQ=255)
                                                                               gatgaGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCAAAAATTTAGGGTTAAtcatca                 <  1:1087063/97‑1 (MQ=255)
                                                                                atgaGTGGCAAAAAAATAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCAAAAATTTAGGGTTAATcatcat                >  2:300370/1‑97 (MQ=255)
                                                                                  gaGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCAAAAATTTAGGGTTAATCATCATtg              <  1:524983/97‑1 (MQ=255)
                                                                                   aGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCAAAAATTTAGGGTTAATCATCATtga             >  1:93781/1‑97 (MQ=255)
                                                                                   aGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCAAAAATTTAGGGTTAATCATCATtga             >  1:209874/1‑97 (MQ=255)
                                                                                       gCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCAAAAATTTAGGGTTAATCATCATTGATgaa         <  1:248820/97‑1 (MQ=255)
                                                                                           aaaaTTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCAAAACTTGAGGGTTAATCATCATCGATGAAGAAc     >  2:359084/1‑97 (MQ=255)
                                                                                           aaaaTTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCAAAAATTTAGGGTTAATCATCATTGATGAAGAAc     <  1:887074/97‑1 (MQ=255)
                                                                                           aaaaATAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCAAAAATTTAGGGTTAATCATCATTGATGAAGAAc     <  2:974968/97‑1 (MQ=255)
                                                                                            aaaTTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCAAAAATTTAGGGTTAATCATCATTGATGAAGAACa    >  1:388353/1‑97 (MQ=255)
                                                                                             aaTTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCAAAAATTTAGGGTTAATCATCATTGATGAAGAACAt   <  1:65233/97‑1 (MQ=255)
                                                                                              aTTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCAAAAATTTAGGGTTAATCATCATTGATGAAGAACATg  >  2:1342832/1‑97 (MQ=255)
                                                                                               |                                                                                               
AATGGTTTTACGCTTCAAACGTCGATTTGGTGATGACGTTGCTGTATTACATTCTGGCTTATCTAATGGGGAACGTTATGATGAGTGGCAAAAAATTAGGGATGGTCGTGCGAGAGTAAGTGTTGGTGCAAGGTCAAGTGTGTTCGCACCTTTCAAAAATTTAGGGTTAATCATCATTGATGAAGAACATG  >  BX571857/1222746‑1222936

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: