Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 249,538 T→G 28.8% 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,5380TG28.8% 78.3 / 9.1 76H187P (CAC→CCC) flhAflagellar type III secretion system protein FlhA
Reads supporting (aligned to +/- strand):  ref base T (17/31);  new base G (20/5);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 4.48e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGAGCGGGTTACCCTGTGCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCTCAGGTT  >  NC_007779/249304‑249773
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              
gTGAGCGGGTTACCCTGTGCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGCAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCg                                                                                                                                                                                                                                            <  1:824403/238‑1 (MQ=255)
     cGGGTTACCCTGTGCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGCGCGGCAtt                                                                                                                                                                                                                                       >  1:281854/1‑238 (MQ=255)
      gggTTACCCTGTGCTTTTTCCCCCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTcc                                                                                                                                                                                                                                    <  1:1551649/239‑1 (MQ=255)
         ttACCCTGTGCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCCGCATTCCCgg                                                                                                                                                                                                                                 >  1:1961847/1‑239 (MQ=255)
               tgtgCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGGTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCCGGacca                                                                                                                                                                                                                            >  1:2035528/1‑238 (MQ=255)
               tgtgCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGGCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCCGCGGCTGTTTGCTCATACGCCAGCCGGGAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGGGGCAATCCCCGGacca                                                                                                                                                                                                                            >  1:247681/1‑238 (MQ=255)
                  gCTTTGTCCACCAGCGCCACCAGTTTGTAACCGTGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACCGCAGGGGGGGCATTCCCgtgacccccc                                                                                                                                                                                                                         >  2:1832124/1‑229 (MQ=255)
                  gCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGGGGCCTTCCCGGTaccacac                                                                                                                                                                                                                         >  1:454317/1‑236 (MQ=255)
                  gCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCCCCg                                                                                                                                                                                                                        >  1:1854896/1‑238 (MQ=255)
                  gCTTTGTCCACCAGCGCCACCAGATTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGCAAGGTTTCACAACTGACCTGCTGCTAGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCCCCTCCGCCGCCCGCGGCGGATTGACCATCCGACAGCCCGGAAAACCAAGCAAGGCGCTGAACCACCAAAACGGCAGGGGCGGGCAGTCCCgtacccacg                                                                                                                                                                                                                          >  2:1943153/1‑230 (MQ=255)
                     ttGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCTCCAGCCGGTAAAACCCAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCCCCa                                                                                                                                                                                                                     >  1:195341/1‑238 (MQ=255)
                      tGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGGCTTCCCGGGACCACcccccgg                                                                                                                                                                                                                   >  1:1559265/1‑236 (MQ=255)
                      tGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCACTCCCGGTACCACCCCCAGc                                                                                                                                                                                                                   >  1:72616/1‑239 (MQ=255)
                      tGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGCGCGGCATTCCCGGGACCACCCCCAg                                                                                                                                                                                                                    >  1:1676026/1‑238 (MQ=255)
                      tGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCCCATCCGCCAGCCGGTAAAACCAAGCAGGGCCCTGAACAGCAAAAACCGCCGGGGCGGCCCTCCCCCGACCCCCGcccc                                                                                                                                                                                                                    >  2:1474360/1‑236 (MQ=255)
                      tGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGGTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCAAACCCGGTAACCCCGcccg                                                                                                                                                                                                                    >  2:815603/1‑236 (MQ=255)
                      tGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCACCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGGCAAACCAAGCAGGGCCCCGAACAGCACAAACGGGCGGGGGCGCATTCACCGGCACCCACGCCa                                                                                                                                                                                                                     >  2:125934/1‑238 (MQ=255)
                        tCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCAc                                                                                                                                                                                                                 >  1:1246737/1‑239 (MQ=255)
                           accaGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGCGGGGCCTTTCCGGTACCACCGCCCGCACCaa                                                                                                                                                                                                              >  1:1506089/1‑239 (MQ=255)
                                cgcCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCCATCCCGGGACCACCCCCAGCACGCACATaa                                                                                                                                                                                                         >  1:857398/1‑239 (MQ=255)
                                  ccaccaGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGGGGGCTCCCCGGGACCACCGGCCAGCagaaaaaaat                                                                                                                                                                                                         >  1:1842465/1‑229 (MQ=255)
                                  ccaccaGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGGGGGAATCCCGGTACCCCCGCCAACACGaaaaaaaa                                                                                                                                                                                                        >  1:775484/1‑232 (MQ=255)
                                           tGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCgg                                                                                                                                                                                                                                 >  1:92949/1‑205 (MQ=255)
                                                   aGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGGTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGTAAACCCGCCAGACCGAACATAATCCCGGTGGGGGTGGaggg                                                                                                                                                                                      >  1:1330539/1‑237 (MQ=255)
                                                   aGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGGGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTTCTCATCCGCCAGCCGGTAAAACCCAACAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGGCATTCCCCGAACCCCCCCCCCACCGAACAAAATCCCGGGGCGgcgtgaga                                                                                                                                                                                      >  2:1308754/1‑235 (MQ=255)
                                                     gCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGGACCACCGCCAGCACGAACATAATCCCCGGGGCgggggtgagc                                                                                                                                                                                     >  1:187425/1‑230 (MQ=255)
                                                              gATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGAGGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCCCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGcc                                                                                                                                                                                                                      <  1:154924/196‑1 (MQ=255)
                                                               aTAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGGACCCACGCCCGCACGCACAAAAACCCGGGGGCGGGGGAGAGCACCCACGGGCt                                                                                                                                                                          >  1:607735/1‑239 (MQ=255)
                                                               aTAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGTACCACCGCCCGCACGAACATAAACCCGGTGGCGGGGGAGAGGACCGACGGGcc                                                                                                                                                                          >  1:973135/1‑238 (MQ=255)
                                                                    cTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCCGCACGAACATAATCCCGGGGGCGGGGTAGAGCa                                                                                                                                                                                    >  1:553749/1‑224 (MQ=255)
                                                                       gATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCCTTCCCGGTACCACCGCCAGCACGAAAATAATCCCGGGGGCGGGGGAGAGCACCCACGGGCTTGCcagca                                                                                                                                                                  >  1:216165/1‑239 (MQ=255)
                                                                             cGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGGATTCCCGGTAACACCGCCAGCACGAACATAATCCCGGTGGCGGTGGAGAGCACCGACGGGGTTGCCAGCAGCTggg                                                                                                                                                            >  1:1765312/1‑238 (MQ=255)
                                                                             cGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGTACCACCGCCCGCACCAACAAAAACCCCGGGGGGGGGTATAGCACCCACCGGCTCGCCCCCCGCTGGt                                                                                                                                                            >  2:1098697/1‑239 (MQ=255)
                                                                                aaTGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGTTCAGCGTTTCGAGGCTTTTCTCTTCCCCCTCCGCCGCCAGCGCCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTa                                                                                                                                                          <  1:472351/238‑1 (MQ=255)
                                                                                 aTGGTTTCCCACCTGCCCTGCTGCTCGCTGTTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTCCCGCCCCCGCCGCCAGCGGCTGTTTGCTCACCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACg                                                                                                                                                        <  1:1585608/239‑1 (MQ=255)
                                                                                       tCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGGGGCACTCCCGGGACCCCCCCCAGCACCAACAAAATCCCCGGGGCGGGGTGGAGGCCCCCCACGCTTGCCAGCAGCCGGGTAcgcgcgg                                                                                                                                                   >  2:585294/1‑237 (MQ=255)
                                                                                               gACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGGAGAGCACCCACGGGCTTGCCAGCAGCTGGGTACGGCCGTCATGGGCGa                                                                                                                                          >  1:1235891/1‑239 (MQ=255)
                                                                                               gACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGGGGCGGGGTAGAGCACCCACGGGCTTGCCAGCAGCTGGTTACGCACGGCATGGGCGa                                                                                                                                          >  1:1375497/1‑239 (MQ=255)
                                                                                                                     gTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGGGGCGGGGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGGTATCGCCGTTGTCGCCGAcgcg                                                                                                                    >  1:431859/1‑239 (MQ=255)
                                                                                                                                           gtggttgttcctctGCCCCCCCCGGCGCGGGGGGGTGTTGGCCCCTCCGCCGGCCGGAAAACCCAAGCGGGGCGCTGACAGGAAAAAACGGGAGGGGCGGCATCCCCGGTCCCGCCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCcgcg                                                                                              <  2:2116240/229‑1 (MQ=255)
                                                                                                                                               gCTTTCCCCTCCCGCCCCCGCCGCCAGCGGCTGTTTGCTCACCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTcc                                                                                                                                                                                                                                    <  2:1778938/102‑1 (MQ=255)
                                                                                                                                                   tcctctCCCGCCCCCGCCGCCAGCGGCTTTTTGCTCACCCCCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACCGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGagag                                                                                      <  1:1726709/237‑1 (MQ=255)
                                                                                                                                                    tctcTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGGGGAGagc                                                                                     >  1:1824225/1‑239 (MQ=255)
                                                                                                                                                      tctcCCGCCCCGGCCGCAGGGGGTTGTTTCCTCTCCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGagcag                                                                                   <  1:34480/239‑1 (MQ=255)
                                                                                                                                                      cctcccgccccggccGCCGGCGCCTGTTGCCCCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGagcag                                                                                   <  1:1796471/235‑1 (MQ=255)
                                                                                                                                                        tcccccccccgcggCCGGCGGCTGTTGGCCCAGCCGCCAGCCGGTAAACCCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAgcagca                                                                                 <  1:1474360/231‑1 (MQ=255)
                                                                                                                                                        taccgccccccccCCCAGCGGCTTTTGGCCCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAgcagca                                                                                 <  1:286896/237‑1 (MQ=255)
                                                                                                                                                         tccgcctccgccGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTaccacc                                                                                                                                                                                                                          >  1:1778938/1‑102 (MQ=255)
                                                                                                                                                                     gccAGCGCTTGTTTGCTCACCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTg                                                                    <  1:454764/239‑1 (MQ=255)
                                                                                                                                                                                 atgtcggcccgccgGCGGGAAACACCAAGCAGGGGCAGAAACAGAAAAAACGACGGGGGCGGCTTTCCCGGCACCCCCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGcc                                                        <  2:1842465/231‑1 (MQ=255)
                                                                                                                                                                                         gcgtcggCCGGTAAATCCACGCAGGGTGCTGAGAATCAAAACGGGCGGGTGGGGCGTCCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGGGCACCGACGGGCTGGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGa                                                <  2:1297540/233‑1 (MQ=255)
                                                                                                                                                                                         cccccAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGa                                                <  1:1617083/236‑1 (MQ=255)
                                                                                                                                                                                           cccgccgggTAAACCCGAGAGGGGCGCGGAACAGCAAAACGGGCAGGTGGGGCTTCCCGGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATg                                              <  2:187425/232‑1 (MQ=255)
                                                                                                                                                                                            ccAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGGGGGGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATgg                                             >  1:583746/1‑239 (MQ=255)
                                                                                                                                                                                               gCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCCTTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGtca                                          <  1:199347/239‑1 (MQ=255)
                                                                                                                                                                                                 ggggaaaaCCAAGCAGGGCGCTGAACGACAAAACCGGCGGGGGCGGCCTTCCCGGTCCCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGtcat                                         <  2:775484/234‑1 (MQ=255)
                                                                                                                                                                                                 gggTAAAACCAAGGAGGGCGCGCAACAGCAAAAACGGGAGGTGCGGCGTTCCCGGTCCCACCGCCAGCACGAACATAAGCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGtcat                                         <  1:1098697/237‑1 (MQ=255)
                                                                                                                                                                                                 cGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGtcat                                         <  1:519516/238‑1 (MQ=255)
                                                                                                                                                                                                    tAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCg                                                                                                                                                                                              <  2:1796627/86‑1 (MQ=255)
                                                                                                                                                                                                    tAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCg                                                                                                                                                                                              >  1:1796627/1‑86 (MQ=255)
                                                                                                                                                                                                     aaaaCCACGCAGGGCGCGGAACAGCAAAAGCGGCAGGTGCGGCATTCCCGGTACCCCCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGc                                     <  1:1832124/238‑1 (MQ=255)
                                                                                                                                                                                                       caCCAAGCAGGGCACGAAACAGAAAAAGCGGCGGTGGCGGCGTTCCGGGAGCCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTGGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGcaca                                  <  2:379531/238‑1 (MQ=255)
                                                                                                                                                                                                       aaCCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGcaca                                  <  1:1620813/239‑1 (MQ=255)
                                                                                                                                                                                                        cccaggaggggCGCGGAACAGAAAAACCGGCGGGTGCGGCATCCCCGGTCCCCCCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACAt                                 <  2:72616/231‑1 (MQ=255)
                                                                                                                                                                                                         ccAAGCAGGGCCGTAAACAGCAAAAGCGGCAGGTGCGGCATTCCCGGTCCCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGcaca                                  <  2:413699/237‑1 (MQ=255)
                                                                                                                                                                                                              cAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGc                                     <  1:85880/229‑1 (MQ=255)
                                                                                                                                                                                                                   cgcTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCCCGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGGGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTggg                        >  1:413699/1‑236 (MQ=255)
                                                                                                                                                                                                                   cgcTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGcaca                                  <  1:1301160/227‑1 (MQ=255)
                                                                                                                                                                                                                         aaCAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGGGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGGTACGCACGGCCTGGGCGATATCCCCGTTGTCGCTGACGCGGGGGACGATAATCGCCGCCGCGGGGGAGAGCAGCAGGGGAGGGACGTGCGCCACCAGGGCGGCGCCGACGGGCATCAGCCCATACCGCTggctgccctt                 >  2:85880/1‑231 (MQ=255)
                                                                                                                                                                                                                               aaaagCGGCGGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATc              <  1:1695013/236‑1 (MQ=255)
                                                                                                                                                                                                                               aaaaaGGGGACGTGCGGCATCCCCGGTACACCCGCCAGCACGAACATAATCCCGGTGTCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCAt               <  2:583746/235‑1 (MQ=255)
                                                                                                                                                                                                                                aaaaCGGCAGGTGCGGCATTCCCGGTCCCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGgcg          <  1:585294/239‑1 (MQ=255)
                                                                                                                                                                                                                                  acggcgaGGTGGGGCTTCCCCGTGCCCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCt        <  2:216165/233‑1 (MQ=255)
                                                                                                                                                                                                                                  accggcggGTGCGGCATCCCCGGCCCCACCGCCAGCGCGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCt        <  1:1086796/237‑1 (MQ=255)
                                                                                                                                                                                                                                  accGGCCGGTGCGGAGTTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGGGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCt        <  2:1573407/237‑1 (MQ=255)
                                                                                                                                                                                                                                        aGGTGCGGCCCTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGGCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCTCAGGt   <  1:87053/238‑1 (MQ=255)
                                                                                                                                                                                                                                           tGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGGCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCGCCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCTCAGGtt  <  1:1943153/236‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              
GTGAGCGGGTTACCCTGTGCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCTCAGGTT  >  NC_007779/249304‑249773

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