Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 249,538 T→G 42.3% H187P (CAC→CCC)  flhA ← flagellar type III secretion system protein FlhA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007779249,5380TG42.3% 37.1 / 23.6 64H187P (CAC→CCC) flhAflagellar type III secretion system protein FlhA
Reads supporting (aligned to +/- strand):  ref base T (13/21);  new base G (21/7);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 5.04e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGAGCGGGTTACCCTGTGCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCTCAGGT  >  NC_007779/249304‑249772
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           
gTGGGCGGGTTACCCTGTGCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCg                                                                                                                                                                                                                                          <  1:553146/238‑1 (MQ=255)
  gggCGGGTTACCCTGTGCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGc                                                                                                                                                                                                                                        <  1:876702/236‑1 (MQ=255)
   aGCGGGTTACCCTGTGCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGCGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCCGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCa                                                                                                                                                                                                                                       >  1:552040/1‑238 (MQ=255)
           aCCCTGTGCTTTGTCCACCAGCGCCACCAGTTTGTTACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCCGCCGGTAAAACCAAGCAGGGCGCTGCACAGCCAAAACGGCAGGGGGGGGCATCCCCGGt                                                                                                                                                                                                                               >  2:644609/1‑239 (MQ=255)
                gtgCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTaccac                                                                                                                                                                                                                          >  1:836597/1‑238 (MQ=255)
                 tgatttGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGACTGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTTAACCGCATAAACGGCCGGGGGCGCCAGCCCGGTaccagc                                                                                                                                                                                                                         >  1:861582/4‑236 (MQ=255)
                 tgCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCCGACTTTCCGGGaccacc                                                                                                                                                                                                                         >  1:1300072/1‑238 (MQ=255)
                 tgCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAACCAGGGCGCGGACCAGCAAAAACGGCGGGGGCGCGATTCCCgggcccccc                                                                                                                                                                                                                         >  2:559685/1‑231 (MQ=255)
                  gCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCCGGTGCGGGCTTCCCCGGagcaccc                                                                                                                                                                                                                        >  2:1353102/1‑237 (MQ=255)
                  gCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCCATCCCGGGACCACCg                                                                                                                                                                                                                        >  1:1472239/1‑238 (MQ=255)
                  gCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAAAGGCAGGGGGGGGAATCCCCGGAccgccg                                                                                                                                                                                                                        >  2:193048/1‑238 (MQ=255)
                      tGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGGATTCCCGGTACCACCGCCAGc                                                                                                                                                                                                                   >  1:1584740/1‑239 (MQ=255)
                                cgcCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGGAATCCCGGGACCACCCCCACCACGAACATaa                                                                                                                                                                                                         >  1:377889/1‑239 (MQ=255)
                                          ttGTCACCCGGGCTTAAATCTAGTGGGCTTGCTCAGCAGAGGGTTTCCCCAGCTCCCCTGCTGCTCGGGGGTTTGTGTGGGGGGCGGGGTCGGTGTTCCGGGTCTTTCCCCCCCCCCCCCCCCCCACGGCGGTTGTTGCCTCACCCCCCAGCCGGTAAACCCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCGGGTACCACCGCCAGCACGAACATAATCCCGGTgg                                                                                                                                                                                                <  2:280057/239‑1 (MQ=255)
                                          ttGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGGGGCAATACCGGTACCACCCCCAGCACGAACATAATACCGGTggg                                                                                                                                                                                               >  1:2256149/1‑238 (MQ=255)
                                                  gAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGGGGCAATCCCGGGACCACCGCCAGCCCGAACAAAAACCCCGtctcctattcaa                                                                                                                                                                                       >  1:472471/1‑228 (MQ=255)
                                                   aGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCCTTCCCGGTACCACCGCCAGCACGAACATAATCCCGGGGGCgggggcgag                                                                                                                                                                                      >  1:393904/1‑232 (MQ=255)
                                                                aTAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCa                                                                                                                                                                                                                                       <  1:1850265/178‑1 (MQ=255)
                                                                aTAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCa                                                                                                                                                                                                                                       >  2:1850265/1‑178 (MQ=255)
                                                                aTAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCCGTACCCCCGCCAGCACGGACAAAATCCCGGGGGGGGGGTAGAGccctttcctttt                                                                                                                                                                          >  1:976043/1‑228 (MQ=255)
                                                                     cTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGTACCACCCCCAGAACGaa                                                                                                                                                                                                              >  1:1719241/1‑198 (MQ=255)
                                                                      tCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGTACCACCGCCAGCACGAACATAAACCCGGTGGCGGGGGAGAGGACCCACGGGGCTGCca                                                                                                                                                                     >  1:1087641/1‑238 (MQ=255)
                                                                          tCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGGGGCGGGGTAGAGCACCGCCCGGCCTGCcagcag                                                                                                                                                                 >  1:687724/1‑238 (MQ=255)
                                                                              cGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCGATCCCCGTACCCCCCCCCCCCCGAAACAAAACCCGGTGGCGGGTGAGAGGACCGACCGGGCTGCCCGCAGCAgg                                                                                                                                                             >  2:1500864/1‑238 (MQ=255)
                                                                                        tCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCCGCCGGTAAAACCAAGCAGGGCGCCGAACAGCAAAAACGGCCGGGGCGGCCTTCCCCGGGCCCCCCCCCCCCCCAACACAAACACCGGGGCGGGGGAGAGCACCCACCGGGCTGCCACCAGCCTGGTACGccgca                                                                                                                                                   >  2:1851615/1‑234 (MQ=255)
                                                                                           cAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCCTTCCCGGTACCCCCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCCACCGGCTTGCCAGCAGCCGGGTACGCACGTCATg                                                                                                                                               >  1:1541780/1‑239 (MQ=255)
                                                                                             aCTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGCGCGGCATTCCCGGGACCACCCCCAGCCCGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCCGCAGCTGGTTACGCACGTCATggg                                                                                                                                             >  1:1921830/1‑239 (MQ=255)
                                                                                                gACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGCCGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGa                                                                                                                                          >  1:1133701/1‑239 (MQ=255)
                                                                                                                          gggtcggggTCGCCGTTCCGGGGCTTTTCTCTCCCGCCCCGGCCCCCGGCGGCTGTTGGCCCACCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGt                                                                                                                <  1:143919/232‑1 (MQ=255)
                                                                                                                               cGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGGACCCCCCCCAGCACGAACATAATCCCGGGGGGGGGGTAGAGCACCCGCGGGCTTTCCCGCAGCCGGGTACCCCCGGCCTGGGCGGTTACGCCCGTGTCCGCGGCCGCGGGTgcgat                                                                                                           >  2:751707/1‑235 (MQ=255)
                                                                                                                                 ggTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGTACCACCCCCAGCACGAACAAAATCCCGGTGGGGGTGTAGAGCACCCACCGGCGTGCCCGCAGCTGGTTAAGCCCGTCATGGGCGATATCGCCGGTGGCGCGCGCGCgggggtggagaa                                                                                                        >  2:1246040/1‑230 (MQ=255)
                                                                                                                                     aGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCCGTACCACCGCCCACACGAACAAAAACCCGGGGGCGGGGGAGAGCACCCACGGGCGTGCCAGCCGCTGGTTTAGCCCGTCCTGGGCGATATCGCCGTTGTCCCCGACGCGGGGGGCGATAATcg                                                                                                     >  2:2033792/1‑238 (MQ=255)
                                                                                                                                         tttCGAGGTTTTTCTCTCCCGCCCCCGCCGCGGGCGGCGGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATcgccgc                                                                                                 <  1:193048/238‑1 (MQ=255)
                                                                                                                                            cGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGGGGCGGTGGAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGGCATGGGCGATATCGCCGTTGTCGCCGACGCGGGGGACGGTAATCGCCGGcgcg                                                                                             >  1:116377/1‑239 (MQ=255)
                                                                                                                                                       tctgcgccccccgacgCCAGCGCTTGTTTGCTCACCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCGGGGGCGGCATTCCCGGTCCCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGagcag                                                                                  <  1:384570/229‑1 (MQ=255)
                                                                                                                                                              cctccgccGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGAGCGGCATTCCCGGTACCACCGCCAGCACGAAAATAATCCCGGGGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATgggg                                                                                                                                            >  1:1743950/1‑174 (MQ=255)
                                                                                                                                                               ctccgccGCCGGCGGCGGTTCGCCAACCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGACCAGCAAAAACGGCGGGGGCGGCATTCCCGGTCCCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCagggaa                                                                           <  1:2033792/238‑1 (MQ=255)
                                                                                                                                                                               gttgtGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCagggaa                                                                           <  1:572864/218‑1 (MQ=255)
                                                                                                                                                                                        aTCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGTACCACCGCCAGCACCAACAAAAACCCCGTGGCGGGGTAGAGCACCCGCGGGCCTGCCCGCAGCCGGGTACCCCCGGCATGGGCGATATCCCCGTTTTCGCTGCCGCGGGGGACGATAATCTCCCCCGCGGGGGGGCGCAGCAGGGGGGGGGGACGCGCCACCCAGCCGgcgc                                                  >  2:1621496/1‑238 (MQ=255)
                                                                                                                                                                                          cgcgcAGCGGGTAAATCACAGCGGGGCGCAGAACAGCAAAAACGTCAGGGGCGGCGTTCCCGGTCCCCCCTCCAGCCCGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGa                                               <  2:343064/235‑1 (MQ=255)
                                                                                                                                                                                               aGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGGGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGtc                                          >  1:1312450/1‑239 (MQ=255)
                                                                                                                                                                                                 ccGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGtca                                         <  1:1692562/238‑1 (MQ=255)
                                                                                                                                                                                                  cGGTAAAACCAAGCAGGGCGCTAAAAAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGtcatc                                       <  1:1246040/239‑1 (MQ=255)
                                                                                                                                                                                                   ggTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGAATTCCCGGTCCCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATcgccgccg                                                                                               <  1:1394997/182‑1 (MQ=255)
                                                                                                                                                                                                        acccAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTCCCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGcaca                                 <  1:486796/237‑1 (MQ=255)
                                                                                                                                                                                                           cAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATAc                              <  1:751707/239‑1 (MQ=255)
                                                                                                                                                                                                              gCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGGGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGCTGGTCATCAGCACATActgc                           >  1:240901/1‑239 (MQ=255)
                                                                                                                                                                                                                    cgcTGAACAGCAAAACCGGCAGGTGGGGCATTCCCGGTACCCCCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGaa                      <  1:255834/238‑1 (MQ=255)
                                                                                                                                                                                                                    cgcTGAACAGCAAAAACGGCAGGTGGGGCATTCCCGTTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGaa                      <  1:644609/238‑1 (MQ=255)
                                                                                                                                                                                                                         aaCAGCAAAAACGGCAGGTGCGGCCTTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCCTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTTTCGCTGGCGCGGGGGACGATAATCGCCCCCCCGGTGGGGCGCAGGAGGGGGAGGATCTGCGCCACCAGGCCCTCGCCGCTGGGCATCATGACAAACTCCTGGGAGgcagc                >  2:543185/1‑239 (MQ=255)
                                                                                                                                                                                                                           ccgCAAAAACGGCAGGGGCGGCATTCCCGGTCCCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCAt              <  2:1133701/237‑1 (MQ=255)
                                                                                                                                                                                                                            aGCAAAAACGGCAGGTGCGGCATTCCCGGTACCCCCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGagag                                                                                     <  1:1581176/167‑1 (MQ=255)
                                                                                                                                                                                                                              cAAAACCGGCAGGTGGGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATcgccgccg                                                                                               <  2:15440/155‑1 (MQ=255)
                                                                                                                                                                                                                              cAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGGGTAGAGCACCGACGGGCTTGCCAGCAGCCGGGTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGGGACGATAATCGCCGCCGCGGGGGAGAGCCGCCGGGGAGGGAGCTGCGCCCCCAGGCCCGCGCCGGTGGGCATCAGCACAAACCGCTGGAAGGCAGCATccg           >  2:192605/1‑237 (MQ=255)
                                                                                                                                                                                                                              cAAAAACGGCAGGGGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTg                                                                                          <  1:94302/160‑1 (MQ=255)
                                                                                                                                                                                                                                 aaaCGGCAGGTGCGGCATTCCCGGTACCACCGCCATCACGAACATAATCCCGGGGGCGGGGTAGAGCACCCACGGGCTTGCCCGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGGGACGATAATCGCCGCCGCGGGGGAGAGCAGCAGGGGAGGGAGCTGCGCCACCCGGCCCTCCCCCAGGGTCATCAGCACATACCGCTGGAAGGGAGCACCGgcg         >  2:648888/1‑238 (MQ=255)
                                                                                                                                                                                                                                  accGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCt       <  1:1621496/237‑1 (MQ=255)
                                                                                                                                                                                                                                  accGGCAGGTGCGGCATCCCCGGTCCCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCt       <  2:240901/237‑1 (MQ=255)
                                                                                                                                                                                                                                  aaCGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCt       <  1:1586026/239‑1 (MQ=255)
                                                                                                                                                                                                                                  aaCGGCAGGTGCGGCATCCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCt       <  1:559685/239‑1 (MQ=255)
                                                                                                                                                                                                                                  aaCGGCAGGGGGTGTATTCCCGGTCCCCCCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATTCTGCTGGAAGGCAGCATCGGCGCt       <  1:543185/239‑1 (MQ=255)
                                                                                                                                                                                                                                       cAGGGGCTGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCTCAgg   <  1:1089156/238‑1 (MQ=255)
                                                                                                                                                                                                                                        aGGTGCGGCATTCCCGGTACCACCGCCAGCCAGAACATAATCCCGGTGGCGGTGTAGAGCACCGACCGGCTTGCCAGCAGCTGGTTACGCCCGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGagag                                                                                     >  1:15440/1‑155 (MQ=255)
                                                                                                                                                                                                                                        aGGTGCGGCATTCCCGGCACCCCCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGACGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCTCAGGt  <  1:1691765/238‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           
GTGAGCGGGTTACCCTGTGCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCTCAGGT  >  NC_007779/249304‑249772

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