Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BX571857 1,167,168 A→G 9.4% D28G (GAT→GGT)  SAS1093 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BX5718571,167,1680AG 9.4% 209.8 / 8.5 62D28G (GAT→GGT) SAS1093hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/28);  new base G (4/2);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTTCAATACTCGTGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGC  >  BX571857/1167073‑1167261
                                                                                               |                                                                                             
tttttCAATACTCGTGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGAt                                                                                              >  1:770436/1‑97 (MQ=255)
    tCAATACTCGTGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCAcc                                                                                          <  1:741259/97‑1 (MQ=255)
     cAATACTCGTGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCa                                                                                         <  2:1103342/97‑1 (MQ=255)
       ttACTCGGGGGGCCACAACAGAGGGAAATAGGATCACCAATTCCAACAGACCATGCAAGTTGGCGGGGCCGCAACATAGAGAAATTGGATCACCAAt                                                                                       <  2:1071989/96‑1 (MQ=14)
         aCTCGTGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAAttt                                                                                     >  2:1411742/1‑97 (MQ=255)
         aCTCGTGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAACGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAAttt                                                                                     >  2:803937/1‑97 (MQ=255)
          cTCGTGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTc                                                                                    >  2:767955/1‑97 (MQ=255)
          cTCGTGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTc                                                                                    <  2:407506/97‑1 (MQ=255)
          cTCGTGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTc                                                                                    <  2:407381/97‑1 (MQ=255)
          cTCGTGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTc                                                                                    >  1:1253107/1‑97 (MQ=255)
          cTCGTGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTc                                                                                    >  1:814996/1‑97 (MQ=255)
            cGTGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCaa                                                                                  >  1:111476/1‑97 (MQ=255)
            cGTGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCaa                                                                                  >  1:968125/1‑97 (MQ=255)
            cGTGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCaa                                                                                  <  1:1510331/97‑1 (MQ=255)
             gTGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAAc                                                                                 >  1:1507965/1‑97 (MQ=255)
              tGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACa                                                                                <  2:24274/97‑1 (MQ=255)
              tGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACa                                                                                >  1:123278/1‑97 (MQ=255)
              tGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACa                                                                                <  2:861150/97‑1 (MQ=255)
               ggggCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAg                                                                               >  2:630817/1‑97 (MQ=255)
                gggCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGa                                                                              >  1:837243/1‑97 (MQ=255)
                gggCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTGGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGa                                                                              <  2:271980/97‑1 (MQ=255)
                 ggCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGAc                                                                             <  2:1575655/97‑1 (MQ=255)
                       aacaGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCa                                                                       <  1:1465029/97‑1 (MQ=39)
                          agagagAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGt                                                                    >  1:553771/1‑97 (MQ=38)
                          agagagAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGt                                                                    <  1:116568/97‑1 (MQ=39)
                           gagagaAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGtt                                                                   >  1:87217/1‑97 (MQ=38)
                            agagaAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTtg                                                                  <  1:794671/97‑1 (MQ=39)
                             gagaAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTtgg                                                                 >  2:1326447/1‑97 (MQ=38)
                              cgaAATAGGATCACCAATTCCAACCGACAATGCAAGCTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAACGCAAGTtggg                                                                <  2:1191683/96‑1 (MQ=11)
                                  aTAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAAGGCAAGTtggggggg                                                            >  2:493584/1‑97 (MQ=255)
                                   tAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTtggggtggg                                                           <  2:515969/97‑1 (MQ=255)
                                   tAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTtggggtgg                                                            >  1:319535/1‑96 (MQ=40)
                                       aTCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGAccc                                                       >  2:598365/1‑97 (MQ=255)
                                       aTCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACAGAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGccc                                                       <  2:977588/97‑1 (MQ=255)
                                        tCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGcctc                                                      >  2:16174/1‑95 (MQ=255)
                                         cACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGGAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCaa                                                     <  2:446888/97‑1 (MQ=255)
                                          aCCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAac                                                    >  1:96290/1‑97 (MQ=255)
                                                   aaCAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCt                                           <  2:24061/97‑1 (MQ=255)
                                                   aaCAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCt                                           <  1:1251/97‑1 (MQ=255)
                                                   aaCAGACAATGCAAGTTGGCGGGCCCCCAACATAGCGAAATTGGATCACCAATTTCAACAGACAAGGCAAGTTGGGGTGGGGCCCAACACAGAAGCt                                           <  1:2292/97‑1 (MQ=255)
                                                     cAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTgg                                         >  1:1416574/1‑97 (MQ=255)
                                                      aGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGc                                        >  1:404717/1‑97 (MQ=255)
                                                         cAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCCGCACAGAAGCTGGCGaa                                     >  2:864703/1‑97 (MQ=255)
                                                            tGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAg                                  >  1:45439/1‑97 (MQ=255)
                                                               aaGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCa                               <  2:842031/97‑1 (MQ=255)
                                                                 gTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCCACACAGAAGCTGGCGAATAGTCAGc                             >  1:73084/1‑97 (MQ=255)
                                                                 gTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGc                             <  1:226201/97‑1 (MQ=255)
                                                                  ttGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAATAGTCAGCa                            >  1:78045/1‑97 (MQ=255)
                                                                  ttGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCa                            <  2:149083/97‑1 (MQ=255)
                                                                  ttGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCa                            >  1:1094891/1‑97 (MQ=255)
                                                                         ggCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACaaaa                     <  1:1258913/97‑1 (MQ=255)
                                                                          gCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACaaaaa                    <  2:279239/97‑1 (MQ=255)
                                                                          gCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACaaaaa                    <  1:138910/97‑1 (MQ=255)
                                                                            cccAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAAtg                  >  1:1188644/1‑97 (MQ=255)
                                                                               aaCATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATgcgc               >  2:742098/1‑97 (MQ=255)
                                                                               aaCATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGc               >  2:1331768/1‑97 (MQ=255)
                                                                                  aTAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAg            <  1:811251/97‑1 (MQ=255)
                                                                                     gagaAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTg         >  2:1381805/1‑97 (MQ=255)
                                                                                       gaACTTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACGCAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGc       >  1:239913/1‑97 (MQ=255)
                                                                                         aaTTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCgg     <  2:474935/97‑1 (MQ=255)
                                                                                         aaTTGGAGCAACAATTTGAACAGACAAGGCAAGTTGGGGTGAGGCCCCACACAGAAGCTGGAGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCgg     <  2:1425550/97‑1 (MQ=21)
                                                                                            tGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGc  <  1:293352/97‑1 (MQ=255)
                                                                                               |                                                                                             
TTTTTCAATACTCGTGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGC  >  BX571857/1167073‑1167261

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