Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,615,078 A→G 17.0% G919G (GGT→GGC USA300HOU_2480 ← helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,615,0780AG17.0% 94.3 / 15.6 47G919G (GGT→GGCUSA300HOU_2480helicase
Reads supporting (aligned to +/- strand):  ref base A (27/12);  new base G (4/4);  total (31/16)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

CAACTTTACTAACATTCATTTTTCGGTTTACATTAATTTGTTAGATAACGATATATATCATCTCTTACAGCTTTATCCAGTGCTAAATCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTCCTTTTGCTCTATGGGACAAAATTTTCTGAACTTCTTTTGATTCTAGCGTTCTA  >  CP000730/2614941‑2615216
                                                                                                                                         |                                                                                                                                          
cAACTTTACTAACATTCATTTTTCGGTTTACATTAATTTGTTAGATAACGATATATATCATCTCTTACAGCTTTATCCAGTGCTAAATCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAGCCt                                                                                                                                         >  1:186040/1‑141 (MQ=255)
cAACTTTACTAACATTCATTTTTCGGTTTACATTAATTTGTTAGATAACGATATATATCATCTCTTACAGCTTTATCCAGTGCTAAATCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAGCCt                                                                                                                                         >  1:451178/1‑141 (MQ=255)
cAACTTTACTAACATTCATTTTTCGGTTTACATTAATTTGTTAGATAACGATATATATCATCTCTTACAGCTTTATCCAGTGCTAAATCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAGCCt                                                                                                                                         >  1:343131/1‑141 (MQ=255)
cAACTTTACTAACATTCATTTTTCGGTTTACATTAATTTGTTAGATAACGATATATATCATCTCTTACAGCTTTATCCAGTGCTAAATCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAGCCt                                                                                                                                         >  1:87115/1‑141 (MQ=255)
             attcattTTTCGGTTTACATTAATTTGTTAGATAACGATATATATCATCTCTTACAGCTTTATCCAGTGCTAAATCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCcag                                                                                                                            >  1:148289/1‑141 (MQ=255)
                   ttttCGGTTTACATTAATTTGTTAGATAACGATATATATCATCTCTTACAGCTTTATCCAGTGCTAAATCCATCGTTCCTATTTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTc                                                                                                                      >  1:354797/1‑141 (MQ=255)
                   ttttCGGTTTACATTAATTTGTTAGATAACGATATATATCATCTCTTACAGCTTTATCCAGTGCTAAATCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTc                                                                                                                      >  2:414977/1‑141 (MQ=255)
                   ttttCGGTTTACATTAATTTGTTAGATAACGATATATATCATCTCTTACAGCTTTATCCAGTGCTAAATCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTc                                                                                                                      >  1:154572/1‑141 (MQ=255)
                        ggTTTACATTAATTTGTTAGATAACGATATATATCATCTCTTACAGCTTTATCCAGTGCTAAATCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTaaa                                                                                                                 >  1:386319/1‑141 (MQ=255)
                        ggTTTACATTAATTTGTTAGATAACGATATATATCATCTCTTACAGCTTTATCCAGTGCTAAATCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTaaa                                                                                                                 >  1:98695/1‑141 (MQ=255)
                        ggTTTACATTAATTTGTTAGATAACGATATATATCATCTCTTACAGCTTTATCCAGTGCTAAATCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTaaa                                                                                                                 >  1:361528/1‑141 (MQ=255)
                           ttACATTAATTTGTTAGATAACTATATATATCATCTCTTACAGCTTTATCCAGTGCTAAATCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAaataa                                                                                                              <  2:32423/141‑1 (MQ=255)
                           ttACATTAATTTGTTAGATAACGATATATATCATCTCTTACAGCTTTATCCAGTGCTAAATCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAaataa                                                                                                              <  2:323079/141‑1 (MQ=255)
                                         tAGATAACGATATATATCATCTCTTACAGCTTTATCCAGTGCTAAATCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATAcc                                                                                                >  2:182865/1‑141 (MQ=255)
                                         tAGATAACGATATATATCATCTCTTACAGCTTTATCCAGTGCTAAATCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATAcc                                                                                                >  2:352666/1‑141 (MQ=255)
                                         tAGATAACGATATATATCATCTCTTACAGCTTTATCCAGTGCTAAATCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATAcc                                                                                                >  2:367546/1‑141 (MQ=255)
                                          aGATAACGATATATATCATCTCTTACAGCTTTATCCAGTGCTAAATCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACca                                                                                               >  1:438786/1‑141 (MQ=255)
                                            aTAACGATATATATCATCTCTTACAGCTTTATCCAGTGCTAAATCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACcatc                                                                                             >  2:37357/1‑141 (MQ=255)
                                                        atcatcTCTTACAGCTTTATCCAGTGCTAAATCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCtttttt                                                                                 >  1:161785/1‑141 (MQ=255)
                                                        atcatcTCTTACAGCTTTATCCAGTGCTAAATCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCtttttt                                                                                 >  1:161813/1‑141 (MQ=255)
                                                               ctTACAGCTTTATCCAGTGCTAAATCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACa                                                                          <  2:291416/141‑1 (MQ=255)
                                                               ctTACAGCTTTATCCAGTGCTAAATCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACa                                                                          >  2:296543/1‑141 (MQ=255)
                                                                                tGCTAAATCCATAGTTACTACGTTTTAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTAAAAAGTAAAAAAATATACATTTTAAt                                                         <  2:397774/141‑1 (MQ=255)
                                                                                  cTAAATCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTc                                                       <  1:137403/141‑1 (MQ=255)
                                                                                     aaTCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTCCtt                                                    >  1:52744/1‑141 (MQ=255)
                                                                                      aTCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTCCttt                                                   >  1:439188/1‑141 (MQ=255)
                                                                                      aTCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTCCttt                                                   >  2:34219/1‑141 (MQ=255)
                                                                                       tCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTCCtttt                                                  >  1:180350/1‑141 (MQ=255)
                                                                                              ttACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTCCTTTTGCTCTAt                                           >  1:272249/1‑141 (MQ=255)
                                                                                                                  ggCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTCCTTTTGCTCTATGGGACAAAATTTTCTGAAct                       >  2:397907/1‑141 (MQ=255)
                                                                                                                   gCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTCCTTTTGCTCTATGGGACAAAATTTTCTGAActt                      <  2:161813/141‑1 (MQ=255)
                                                                                                                   gCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTCCTTTTGCTCTATGGGACAAAATTTTCTGAActt                      <  2:161785/141‑1 (MQ=255)
                                                                                                                     aTTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTCCTTTTGCTCTATGGGACAAAATTTTCTGAActtct                    >  2:157561/1‑141 (MQ=255)
                                                                                                                     aTTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTCCTTTTGCTCTATGGGACAAAATTTTCTGAActtct                    >  2:33860/1‑141 (MQ=255)
                                                                                                                       tttATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTCCTTTTGCTCTATGGGACAAAATTTTCTGAACTTCttt                  <  1:182865/141‑1 (MQ=255)
                                                                                                                       tttATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTCCTTTTGCTCTATGGGACAAAATTTTCTGAACTTCttt                  <  1:352666/141‑1 (MQ=255)
                                                                                                                       tttATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTCCTTTTGCTCTATGGGACAAAATTTTCTGAACTTCttt                  <  1:367546/141‑1 (MQ=255)
                                                                                                                          aTCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTCCTTTTGCTCTATGGGACAAAATTTTCTGAACTTCTTTTGa               >  1:83742/1‑141 (MQ=255)
                                                                                                                          aTCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTCCTTTTGCTCTATGGGACAAAATTTTCTGAACTTCTTTTGa               >  1:10222/1‑141 (MQ=255)
                                                                                                                          aTCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTCCTTTTGCTCTATGGGACAAAATTTTCTGAACTTCTTTTGa               >  1:384329/1‑141 (MQ=255)
                                                                                                                            cTTGCTTTTCTGAGCCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTCCTTTTGCTCTATGGGACAAAATTTTCTGAACTTCTTTTGAtt             <  2:186040/141‑1 (MQ=255)
                                                                                                                            cTTGCTTTTCTGAGCCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTCCTTTTGCTCTATGGGACAAAATTTTCTGAACTTCTTTTGAtt             <  2:87115/141‑1 (MQ=255)
                                                                                                                            cTTGCTTTTCTGAGCCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTCCTTTTGCTCTATGGGACAAAATTTTCTGAACTTCTTTTGAtt             <  2:343131/141‑1 (MQ=255)
                                                                                                                            cTTGCTTTTCTGAGCCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTCCTTTTGCTCTATGGGACAAAATTTTCTGAACTTCTTTTGAtt             <  2:451178/141‑1 (MQ=255)
                                                                                                                             ttGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTCCTTTTGCTCTATGGGACAAAATTTTCTGAACTTCTTTTGATTc            <  2:180350/141‑1 (MQ=255)
                                                                                                                                    tCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTCCTTTTGCTCTATGGGACAAAATTTTCTGAACTTCTTTTGATTCTAGCGtt     >  1:82459/1‑141 (MQ=255)
                                                                                                                                       gAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTCCTTTTGCTCTATGGGACAAAATTTTCTGAACTTCTTTTGATTCTAGCGTTCTa  <  1:397907/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                          
CAACTTTACTAACATTCATTTTTCGGTTTACATTAATTTGTTAGATAACGATATATATCATCTCTTACAGCTTTATCCAGTGCTAAATCCATCGTTACTACGTTTGAACCATTCGGCATTTTATCTTGCTTTTCTGAACCTTCTATATATCCAGCAGTTCCTAAATAATAAAAATATATACCATCATCATCTTTTTTCTGTACAAAAATATACATTTTAATTCCTTTTGCTCTATGGGACAAAATTTTCTGAACTTCTTTTGATTCTAGCGTTCTA  >  CP000730/2614941‑2615216

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