Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,201,277 A→C 19.0% G137G (GGT→GGG zraP ← Zn‑dependent periplasmic chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,201,2770AC19.0% 158.8 / 3.8 71G137G (GGT→GGGzraPZn‑dependent periplasmic chaperone
Reads supporting (aligned to +/- strand):  ref base A (49/6);  new base C (3/13);  total (52/19)
Fisher's exact test for biased strand distribution p-value = 1.92e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGTTCGCGGCTAACAGGGCATTGTATTCATAACGCTTCGTCACCAGTTGCTGTTGCAGTGCGCTGCTT  >  NC_000913/4201046‑4201514
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             
gcTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCCTGGCTGGAACCGCCCGAAAGTTCGCAATTACTGTTAGTGGCAAACCCTGTTTTCTGTCGCTGGCAACCAAACGGACAAACGGCTTGATTTCCACGTGGcccacccccataa                                                                                                                                                                                                                                        >  1:1312075/1‑238 (MQ=255)
gcTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAATGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGccacaacccatag                                                                                                                                                                                                                                        >  2:892099/1‑238 (MQ=255)
 cTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGcccataccccgtt                                                                                                                                                                                                                                        >  2:203410/1‑234 (MQ=255)
 cTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATg                                                                                                                                                                                                                                       >  2:58146/1‑239 (MQ=255)
  tGGCATGCGATGGAACAGGCCGTGACCACCTGGGAAGGGCAACCCATTAAACCTGATTTAGATTCGGAGGCGATAGCGCATTTCGAGCGCCGCCAGACCCGCTTTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATg                                                                                                                                                                                                                                       <  2:963658/238‑1 (MQ=255)
       gggggtgaaactgcacggcacacccTGCGAAGGGCAAACGATTAACCCTAATTGAGATGCGCAGGCAATAGGGCACAGCGGGGGCCGCCAGGGCAGCTTTTCTGGTGATTTCAGCGTGGCATGGCTGGAAGCGACCGAAGGTTCGCATTTACGGTTAGGGGAAACCTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCGCTACCCATATGAgcgc                                                                                                                                                                                                                                  <  1:524423/215‑1 (MQ=255)
           gggtaaCAGTCCGGCACCCCCGGGGAGGGGAAACCCATCAACCCTGATTTAGATGCGCAGGCGTTAGGCAATTGCGGGCGCCCGCAGGGCCGCTCTTCTGTTGAGGTCAGCGGGGCATGGCGGGAAGCGCCCGAAGGTCCGAAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGccacc                                                                                                                                                                                                                              <  1:772623/234‑1 (MQ=255)
                    aCCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAgcc                                                                                                                                                                                                                     >  2:633352/1‑238 (MQ=255)
                    aCCGTGACCACCTGCGAAGGGCAAACCATTAAACCGGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGACGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCGGTCAACCCAACCAGAAAAGGGTCTGTTTTACAAATGGcccctacccaacagacacccccaaacccaagc                                                                                                                                                                                                                     >  1:1104878/1‑216 (MQ=255)
                      cGTGACCACCTGCGAAGGGCCAACCATTAAACCTGATTTAGATTCGCAGGCGAGAGCGCATATCGAGCGCCGCAAGAGCCGCTCTTCTGTTGATGTCAGCGGGGCATGGCGGGAACCGCCGGAAGGCTCGCACTTACTGTTAGTTGCACACGCTGCTTTCTGTCCCAGGCACGCAAACAAGGAACGCTTGTGGTTTACCAGTGGGCCCCCCCCATATGaacacaaccccagcaccacgg                                                                                                                                                                                                                  >  1:1058289/1‑219 (MQ=255)
                      cGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCCCCTCAgccgcc                                                                                                                                                                                                                  >  2:882541/1‑239 (MQ=255)
                      cGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTCATTTACCAGTGGCCAATACCCCTTTTACCGCccccactcacaccccc                                                                                                                                                                                                                  >  2:1141487/1‑229 (MQ=255)
                       ggcacccatgtagaatgcaatagcATTAAACGTTATTTAGGTGCGGAGGTGATCGCGATAATGGAGCGCCTCCAGACCCGCTCTGCTTTTGATGTCAGTGGGGCATGGCTGGAAGCGCCGAAAGGTCCGAAATTACTTTTGTTGCAAATCTCTGATTTCTGTCCCTGGCAACCCAACGAGAAAACGTTCTGTTTTACCATTGCCCCATCCCCATTTGACCGCCACCTCCGCAGccgccg                                                                                                                                                                                                                 <  1:1196086/216‑1 (MQ=255)
                        tGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCACACTCTGATTTCTGTCGCTGGCAACCCACCGAGAAAACGTTCTGATTTACCACTTGCCACTACCCATATGACCCCCaccacccaccccccccc                                                                                                                                                                                                                >  1:292411/1‑228 (MQ=255)
                        tGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGt                                                                                                                                                                                                                >  2:229338/1‑239 (MQ=255)
                        tGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCACTGGCCCATACCCATATGACCGCCACAACCGCAACCGCCGt                                                                                                                                                                                                                >  2:967991/1‑239 (MQ=255)
                         gACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGGTTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGccgccgga                                                                                                                                                                                                               >  2:666237/1‑237 (MQ=255)
                         gACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTCTTGATTTTCCAGTGGCCCCTCCCCATATGACGCCCACCACCACGGCAACCGTa                                                                                                                                                                                                               >  2:194377/1‑239 (MQ=255)
                         gACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAAAGTTCTGATTTACCAGTGGCCCATACCCATAGGACCGCCACCACCGGACCCGCCGTa                                                                                                                                                                                                               >  2:987011/1‑239 (MQ=255)
                                  cGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGCAAGCGACCGAAGGTTCGCAATTACTGTTACTGGCCAACTCCGAGTCCTGCCCCTGGCAACCCAACGGGAAAACGTTTTGATTTACCAGAGGCCCCAACCCCAAACAAGCCCACCACCGCAgcccaccaagccccatc                                                                                                                                                                                                       >  1:231515/1‑224 (MQ=255)
                                   gAAGGGAAAACCATTAACCCTGATTTAGTTTCGCGGGCGATAGCGCATATCGAGCGCCGCAAGGGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCCCCCGAAGGTTCGAAATTACTGTTAGTGGAAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGGGGCCCATCCCCATATGACCGCCACCACCGCAGCCGCCGTAgcccatgcc                                                                                                                                                                                                      <  2:1008630/238‑1 (MQ=255)
                                         cAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGCCCGCCACCACCGCAGCCGCCGTAACCCATGCCCATTcc                                                                                                                                                                                                >  2:757656/1‑238 (MQ=255)
                                                 tAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACCCTGATTTCTGTCGCTGGCAACCCAACGAGCAAACGATCTGATTTACCAGTGGCCCATCCCCATATGACCCCCCCCACCGCCTCCGCCGCAACCCACTCCCATTCCCGCCCcccg                                                                                                                                                                                       >  1:1017887/1‑239 (MQ=255)
                                                        gATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGACGCGCCCGAAGGTTCGCACTTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTCCCAGTGGCCCATACCCATATGACCCCCACACCCCCACCCCCCGAACCCCCTCCCCCTTCCgcgtcccccccgaaa                                                                                                                                                                                 >  1:1194786/1‑224 (MQ=255)
                                                        gATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAAt                                                                                                                                                                                 >  2:640485/1‑238 (MQ=255)
                                                         aTTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATa                                                                                                                                                                                >  2:561122/1‑238 (MQ=255)
                                                               aTTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCCCCGCAGCCGCCGTACCCCACGCCCATTCCgcgcccccccggaaaaacgcct                                                                                                                                                                          >  1:150213/1‑217 (MQ=255)
                                                                   gCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATTCCCATTCCGGCTCCGCCCGGAATACCCGCCTCAg                                                                                                                                                                      >  2:254367/1‑238 (MQ=255)
                                                                          aTAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAACCGTTCTGATTTACCACTGGCCCCTACCCATATGACCGCCACCACCGCAGCCGCCGGAGCCCATACCCCTTCCGCCTCCCCACGGCAAAACCCCTTccccccccccg                                                                                                                                                              >  1:910278/1‑229 (MQ=255)
                                                                                   atCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCCTGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGCCGCTGGCAACCCAACGAGAAAACGTTGTGATTTACCAGTGGCCCATACCCATATGACCACCCACACACCAACCGCCGTAGCCAATCCCCATACCGGCCCCGCGCGGAAAACCCCCCTCAcccccaccaaaaccacgt                                                                                                                                                     >  1:612040/1‑224 (MQ=255)
                                                                                   atCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCCCGCGGAATACCCGCTTCAGCCATCGCGATATCCCGt                                                                                                                                                     >  2:64311/1‑239 (MQ=255)
                                                                                       aGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCCGGCAACCCAACGAGAAACGGTTCTGATTTACCAGTGGCCCATACCCATATGAACGCACCCACCGCCGCCGCCGTAGCCCCTGCCCATTCCGGCTCCGCGCCGAATAACCCCCTCAACCCCCCCGGTAACCCGTTTc                                                                                                                                                  >  1:717691/1‑238 (MQ=255)
                                                                                        gcgcCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCa                                                                                                                                                 >  2:28360/1‑238 (MQ=255)
                                                                                         cgccgcCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATCCCGGCTCCGCGCGGGATACCCGCTTCAGCCATCGCGATATCTCGTTTCAc                                                                                                                                                >  2:664313/1‑238 (MQ=255)
                                                                                                                 gATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAAATCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACGGCACCAGCAGCAGCCCATGCCCATTCCGGCTCCGCCCCGAATACCCGCTTCCGCCAACGAGATATCTCGTTTCAACCGCaaacaacacaaaaaccgacgc                                                                                                                       >  1:633677/1‑220 (MQ=255)
                                                                                                                    gTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCCTTCCGGCTCCGCGCGGCATACCCCCTTCAGCCATCGCGATATCTCGTTTCACCCGTACACCATCTAACGACTGACGCa                                                                                                                      >  2:1148632/1‑237 (MQ=255)
                                                                                                                    gTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCaaa                                                                                                                    >  2:175000/1‑239 (MQ=255)
                                                                                                                                 ggCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCtt                                                                                                       >  2:461550/1‑239 (MQ=255)
                                                                                                                                   cTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAACCCATCCCCATTCCGGCTCCGCCCGGAATACCCCCCTCCGCCATCGCGAAATCTCGTTTCACCCCTAATCCATTTAACGCCTGCCCCAAATTCTCCACTCtttt                                                                                                     >  1:1043659/1‑237 (MQ=255)
                                                                                                                                   cTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTCACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGCATACCCCCTTCAGCCATCGCGATATCTCGTTTCACCCCTAACTCATCTACCGACTGTCGCACATTCTcccccccttt                                                                                                      >  2:331476/1‑230 (MQ=255)
                                                                                                                                            cccGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGCTTTCTGTCGCTGGCAACCCAACGAGACAACGTTCTGATTTCCCAGTGGCCCCTACCCATATGACCGCCAACACCCCAGCCGCCGTAGCCCATGCCCATTACGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCCCGTTTCACCCGTAACTCATCTAACGCCTGCCGCAAATTCTCCCTCTCTTTGGcacccg                                                                                              >  1:885165/1‑232 (MQ=255)
                                                                                                                                            cccGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACcccg                                                                                            >  2:1070802/1‑239 (MQ=255)
                                                                                                                                                    ttCGTAATTCCGGTTAGGGGCAACCTTTGATTTCTGGCGCTGGACACCCAACGAGAAAACGTTCTGATTTCCCAGTGGCCCATCCCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAAttt                                                                                     <  2:1104878/238‑1 (MQ=255)
                                                                                                                                                      cGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCCCGCTACCTCGTTTCACCCGTAAATCATCTAACGACTACCCCCAATTCCCCAACTTTTTGGCGAccccctcacatttcc                                                                                  >  1:138207/1‑226 (MQ=255)
                                                                                                                                                           ttgcTGTTAGTGGAAACCTCTGTTTTTGGTCGCTGCCACCCCAACGAAAAAACGTTCTGATTTACCAGGGGCCCATCCCCATATGCCCCCCCCCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTa                                                                             <  2:1312075/236‑1 (MQ=255)
                                                                                                                                                           cacaggTTAGTGGCAATCTCTGTTTCCGGTCGCTGGCAACCCAACGATAAACCGTTCTTATTTACCAGTGGCCCATCCCAATATGCCCGCCCCCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCAACCGCGTTAATTTtgctgc                                                                               <  2:622283/232‑1 (MQ=255)
                                                                                                                                                                gTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGTTGCTATCCgg                                                                        >  2:635480/1‑239 (MQ=255)
                                                                                                                                                                   aGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGGAGCCCAGGCCCATTCCGGCTCCGCGCGGAATACCCCCTTCCGCCAACGCGATATCCCGTTTCACCCGTAACTCATCTAACGTCTCACGCCCATGCCCCCTCTCTTTTGCGACCCAGTTAATTTtgttgcgaccgggggg                                                                     >  1:672223/1‑229 (MQ=255)
                                                                                                                                                                       gCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCCCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGtt                                                                  >  2:1158439/1‑238 (MQ=255)
                                                                                                                                                                        cAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGTTcg                                                                >  2:910224/1‑239 (MQ=255)
                                                                                                                                                                           aCTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCCCGGAATACCCGCCTCAGCCATCCCGATATCTCGTTTCACCCGTACCCCACCTAACGACTGCACCAAATTCTCCATCTCTTTGCCACACGCGTTAATTTTCCTCCTATCCgaggggctgtg                                                               >  1:622283/1‑227 (MQ=255)
                                                                                                                                                                                tgattcgttgcGCCGGCCACCCAACGAGAAAACGTTCTGATTCACCAGTGGCCCATCCCCATATGACCGCCACCCCCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGTTCGCGGCTaa                                                         <  2:757472/229‑1 (MQ=255)
                                                                                                                                                                                tGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGTTCTCGGCTaa                                                         >  2:361066/1‑239 (MQ=255)
                                                                                                                                                                                tGATTTCTGCCGCTGACAACCCAACGAAAAACCTTTCTTATTTCCCGGGGCCCAATCCCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGTTCGCGGCTaa                                                         <  2:231515/239‑1 (MQ=255)
                                                                                                                                                                                  ttttGTGTCCCTGGAAACCCAAGGGAAAAACGTTCTGTTTTACCAGTGGCCCATCCCCATAGCACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGTTCGCGGCTAACa                                                       <  2:1043659/238‑1 (MQ=255)
                                                                                                                                                                                  aTTTCGGTCGCGGGCACCCCAAGGAAAAACCTTTCTGATTTCCCAGGGGCCCATCCCCATATGACCGCCCCCCCCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGTTCGCGGCTAACa                                                       <  2:150213/239‑1 (MQ=255)
                                                                                                                                                                                      cTGTCCATGGAACCCCAGCGAAAAACCGTTCTTATTTCCCAGGGCCCCATACCCATATGCCCGCCCCCCCCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGTTCGCGGCTAACAGGGc                                                   <  2:1116084/239‑1 (MQ=255)
                                                                                                                                                                                        gTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCCCCCGTAACTCACCTAACGACTGACGCAAATTTTCCATCTCTTTGGCGCCCGCGTTAATTTTGCTGCATTCCGGAGGGttggggggaaaacgggcat                                                 >  1:297130/1‑222 (MQ=255)
                                                                                                                                                                                        gTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCCCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTTGCGACCGCGTTACTTTTGCGACTATCCGGGGGGTTCCAGGCTAACAGGGCAt                                                 >  1:533641/1‑239 (MQ=255)
                                                                                                                                                                                         ccgcggccAACCAAACGAGAAACCGTTTGGATTACCCAGTGGCCCATCCCCATAGGACCGCCCCCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGTTCGCGGCTAACAGGGCAtt                                                <  2:885165/232‑1 (MQ=255)
                                                                                                                                                                                           ccTGGCACCCCAACGAGAAACCTTTCTGATTTCCCAGTGCCCCATCCCCATAGGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGTTCGCGGCTAACAGGGCATTGt                                              <  2:292411/238‑1 (MQ=255)
                                                                                                                                                                                             gggcacccaaacagaaaaaCGTTCTGATTTCCAGGTGCCCCATACCCATAGGCCCGCCACCCCCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGTTCGCGGCTAACAGGGCATTGTAt                                            <  2:1017887/224‑1 (MQ=255)
                                                                                                                                                                                                      aaCGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATACCTCGTTTCACCCGTAACTCATCTAACACCTGACGCAAATTCCCCACCCCTTTGGGGACCCCGTTATCTTTGCTGCTATCCGGGGGGTTCGCGGCTAATCGGGCATTGTATTCATAAc                                     >  1:1053102/1‑237 (MQ=255)
                                                                                                                                                                                                       aCGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCt                                                                                                                                                                                            >  2:590527/1‑85 (MQ=255)
                                                                                                                                                                                                                  ctcttgttcaCCAGGGGCCCATCCCCATAGGACCGCCACCCGCCCACCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGTTCGCGGCTAACAGGGCATTGTATTCATAACGCTTCGTCACCAGt                       <  2:270182/230‑1 (MQ=255)
                                                                                                                                                                                                                   tCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGTTCGCGGCTAACAGGGCATTGTATTCATAACGCTTCGTCACCAGtt                      >  2:1150180/1‑239 (MQ=255)
                                                                                                                                                                                                                          tACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGTTCGCGGCTAACAGGGCATTGGATTCATAACGCTTCGTCACCAGTTGCTGTTg               >  2:1167438/1‑239 (MQ=255)
                                                                                                                                                                                                                              ggTGCCCCATACCCATATGACCGCCACCCCCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGTTCGCGGCTAACAGGGCATTGTATTCATAACGCTTCGTCACCAGTTGCTGTTGCAGt           <  2:1291871/238‑1 (MQ=255)
                                                                                                                                                                                                                               gTGGCCCATACCCATATGACCGCCACCACCTCAGCAGCCGTAGCCCATGCCCATTCCGGCGCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCCAATTCTCCATCTCTTTTGCGACCGCGTTAATTTTGCTGTTATACGGTGCGTTCGCGGCTACCAGGGCTTTGTATTCATAACGCTTCGTCACCCATTTCTTTTGCAggg          >  1:905002/1‑237 (MQ=255)
                                                                                                                                                                                                                                 ggCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGTTCGCGGCTAACAGGGCATTGTATTCATaa                                      >  2:166434/1‑209 (MQ=255)
                                                                                                                                                                                                                                 ggCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGTTCGCGGCTAACAGGGCATTGGATTCATAACGCTTCGTCACCAGTTGCTGTTGCAGTgc         >  2:23296/1‑238 (MQ=255)
                                                                                                                                                                                                                                 ggCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGTTCGCGGCTAAACGGGCCTTGTTTTCATAACCCTTCGTCACCATTTCCGGTTgagggggg        >  1:1164465/1‑232 (MQ=255)
                                                                                                                                                                                                                                 ggCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCAATCCGGTGGGTTCGCGGCTAACAGGGCATTGTATTCATAACGCTTCGCAACCAGTTGCTGTTGCACTgc         >  2:583732/1‑238 (MQ=255)
                                                                                                                                                                                                                                 ggCCCATACCCATAAGACCGCCACCACAGCAGCCGCCGTAGCCCATGCCCATTCCAGCTCCACGCGGTATACCCGCTTCAGCCAACGCGATATCTCGTTTCACCCGCAACTCATCTAACGCCTGACGCACATTCTCCAACACTTATGCGAACGCCTTAAGTATGTTGCTACCCGACGGCTTCCCGGCTAACTGGGCACTGACTTCATAACGCTTCGTCAcaaataccccttaccacgc         >  1:616498/1‑220 (MQ=255)
                                                                                                                                                                                                                                  gCCCATCCCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGTTCGCGGCTAACAGGGCATTGTATTCATAACGCTTCGTCACCAGTTGCTGTTGCAGTgcg        <  2:1109617/238‑1 (MQ=255)
                                                                                                                                                                                                                                       tacccataTGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGTTCGCGGCTAACAGGGCATTGTATTCATAACGCTTCGTCACCAGTTGCTGTTGCAGTGCGCTGCtt  >  2:312687/1‑239 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             
GCTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAACTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGTTCGCGGCTAACAGGGCATTGTATTCATAACGCTTCGTCACCAGTTGCTGTTGCAGTGCGCTGCTT  >  NC_000913/4201046‑4201514

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