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→GGG) zraPZn‑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

GCTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGTTCGCGGCTAACAGGGCATTGTATTCATAACGCTTCGTCACCAGTTGCTGTTGCAGTGCGCTGCTT  >  NC_000913/4201046‑4201514
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             
gcTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCCTGGCTGGAACCGCCCGAAAGTTCGCAATTACTGTTAGTGGCAAA‑CCCTGTTTTCTGTCGCTGGCAACCAAACGGACAAACGGCTTGATTTCCACGTGGcccacccccataa                                                                                                                                                                                                                                        >  1:1312075/1‑238 (MQ=255)
gcTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAATGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGccacaacccatag                                                                                                                                                                                                                                        >  2:892099/1‑238 (MQ=255)
 cTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGcccataccccgtt                                                                                                                                                                                                                                        >  2:203410/1‑234 (MQ=255)
 cTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATg                                                                                                                                                                                                                                       >  2:58146/1‑239 (MQ=255)
  tGGCATGCGATGGAACAGGCCGTGACCACCTGGGAAGGGCAACCCATTAAACCTGATTTAGATTCGGAGGCGATAGCGCATTTCGAGCGCCGCCAGACCCGCTTTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATg                                                                                                                                                                                                                                       <  2:963658/238‑1 (MQ=255)
       gggggtgaaactgcacggcacacccTGCGAAGGGCAAACGATTAACCCTAATTGAGATGCGCAGGCAATAGGGCACAGCGGGGGCCGCCAGGGCAGCTTTTCTGGTGATTTCAGCGTGGCATGGCTGGAAGCGACCGAAGGTTCGCATTTACGGTTAGGGGAAAC‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCGCTACCCATATGAgcgc                                                                                                                                                                                                                                  <  1:524423/215‑1 (MQ=255)
           gggtaaCAGTCCGGCACCCCCGGGGAGGGGAAACCCATCAACCCTGATTTAGATGCGCAGGCGTTAGGCAATTGCGGGCGCCCGCAGGGCCGCTCTTCTGTTGAGGTCAGCGGGGCATGGCGGGAAGCGCCCGAAGGTCCGAAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGccacc                                                                                                                                                                                                                              <  1:772623/234‑1 (MQ=255)
                    aCCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAgcc                                                                                                                                                                                                                     >  2:633352/1‑238 (MQ=255)
                    aCCGTGACCACCTGCGAAGGGCAAACCATTAAACCGGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGACGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCGGTCAACCCAACCAGAAAAGGGTCTGTTTTACAAATGGcccctacccaacagacacccccaaacccaagc                                                                                                                                                                                                                     >  1:1104878/1‑216 (MQ=255)
                      cGTGACCACCTGCGAAGGGCCAACCATTAAACCTGATTTAGATTCGCAGGCGAGAGCGCATATCGAGCGCCGCAAGAGCCGCTCTTCTGTTGATGTCAGCGGGGCATGGCGGGAACCGCCGGAAGGCTCGCACTTACTGTTAGTTGCACA‑CGCTGCTTTCTGTCCCAGGCACGCAAACAAGGAACGCTTGTGGTTTACCAGTGGGCCCCCCCCATATGaacacaaccccagcaccacgg                                                                                                                                                                                                                  >  1:1058289/1‑219 (MQ=255)
                      cGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCCCCTCAgccgcc                                                                                                                                                                                                                  >  2:882541/1‑239 (MQ=255)
                      cGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTCATTTACCAGTGGCCAATACCCCTTTTACCGCccccactcacaccccc                                                                                                                                                                                                                  >  2:1141487/1‑229 (MQ=255)
                       ggcacccatgtagaatgcaatagcATTAAACGTTATTTAGGTGCGGAGGTGATCGCGATAATGGAGCGCCTCCAGACCCGCTCTGCTTTTGATGTCAGTGGGGCATGGCTGGAAGCGCCGAAAGGTCCGAAATTACT‑TTTGTTGCAAATCTCTGATTTCTGTCCCTGGCAACCCAACGAGAAAACGTTCTGTTTTACCATTGCCCCATCCCCATTTGACCGCCACCTCCGCAGccgccg                                                                                                                                                                                                                 <  1:1196086/216‑1 (MQ=255)
                        tGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCACA‑CTCTGATTTCTGTCGCTGGCAACCCACCGAGAAAACGTTCTGATTTACCACTTGCCACTACCCATATGACCCCCaccacccaccccccccc                                                                                                                                                                                                                >  1:292411/1‑228 (MQ=255)
                        tGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGt                                                                                                                                                                                                                >  2:229338/1‑239 (MQ=255)
                        tGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCACTGGCCCATACCCATATGACCGCCACAACCGCAACCGCCGt                                                                                                                                                                                                                >  2:967991/1‑239 (MQ=255)
                         gACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGGTTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGccgccgga                                                                                                                                                                                                               >  2:666237/1‑237 (MQ=255)
                         gACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTCTTGATTTTCCAGTGGCCCCTCCCCATATGACGCCCACCACCACGGCAACCGTa                                                                                                                                                                                                               >  2:194377/1‑239 (MQ=255)
                         gACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAAAGTTCTGATTTACCAGTGGCCCATACCCATAGGACCGCCACCACCGGACCCGCCGTa                                                                                                                                                                                                               >  2:987011/1‑239 (MQ=255)
                                  cGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGCAAGCGACCGAAGGTTCGCAATTACTGTTACTGGCCAA‑CTCCGAGTCCTGCCCCTGGCAACCCAACGGGAAAACGTTTTGATTTACCAGAGGCCCCAACCCCAAACAAGCCCACCACCGCAgcccaccaagccccatc                                                                                                                                                                                                       >  1:231515/1‑224 (MQ=255)
                                   gAAGGGAAAACCATTAACCCTGATTTAGTTTCGCGGGCGATAGCGCATATCGAGCGCCGCAAGGGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCCCCCGAAGGTTCGAAATTACTGTTAGTGGAAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGGGGCCCATCCCCATATGACCGCCACCACCGCAGCCGCCGTAgcccatgcc                                                                                                                                                                                                      <  2:1008630/238‑1 (MQ=255)
                                         cAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGCCCGCCACCACCGCAGCCGCCGTAACCCATGCCCATTcc                                                                                                                                                                                                >  2:757656/1‑238 (MQ=255)
                                                 tAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CCCTGATTTCTGTCGCTGGCAACCCAACGAGCAAACGATCTGATTTACCAGTGGCCCATCCCCATATGACCCCCCCCACCGCCTCCGCCGCAACCCACTCCCATTCCCGCCCcccg                                                                                                                                                                                       >  1:1017887/1‑239 (MQ=255)
                                                        gATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGACGCGCCCGAAGGTTCGCACTTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTCCCAGTGGCCCATACCCATATGACCCCCACACCCCCACCCCCCGAACCCCCTCCCCCTTCCgcgtcccccccgaaa                                                                                                                                                                                 >  1:1194786/1‑224 (MQ=255)
                                                        gATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAAt                                                                                                                                                                                 >  2:640485/1‑238 (MQ=255)
                                                         aTTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATa                                                                                                                                                                                >  2:561122/1‑238 (MQ=255)
                                                               aTTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCCCCGCAGCCGCCGTACCCCACGCCCATTCCgcgcccccccggaaaaacgcct                                                                                                                                                                          >  1:150213/1‑217 (MQ=255)
                                                                   gCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATTCCCATTCCGGCTCCGCCCGGAATACCCGCCTCAg                                                                                                                                                                      >  2:254367/1‑238 (MQ=255)
                                                                          aTAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAACCGTTCTGATTTACCACTGGCCCCTACCCATATGACCGCCACCACCGCAGCCGCCGGAGCCCATACCCCTTCCGCCTCCCCACGGCAAAACCCCTTccccccccccg                                                                                                                                                              >  1:910278/1‑229 (MQ=255)
                                                                                   atCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCCTGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGCCGCTGGCAACCCAACGAGAAAACGTTGTGATTTACCAGTGGCCCATACCCATATGACCACCCACACACCAACCGCCGTAGCCAATCCCCATACCGGCCCCGCGCGGAAAACCCCCCTCAcccccaccaaaaccacgt                                                                                                                                                     >  1:612040/1‑224 (MQ=255)
                                                                                   atCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCCCGCGGAATACCCGCTTCAGCCATCGCGATATCCCGt                                                                                                                                                     >  2:64311/1‑239 (MQ=255)
                                                                                       aGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCCGGCAACCCAACGAGAAACGGTTCTGATTTACCAGTGGCCCATACCCATATGAACGCACCCACCGCCGCCGCCGTAGCCCCTGCCCATTCCGGCTCCGCGCCGAATAACCCCCTCAACCCCCCCGGTAACCCGTTTc                                                                                                                                                  >  1:717691/1‑238 (MQ=255)
                                                                                        gcgcCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCa                                                                                                                                                 >  2:28360/1‑238 (MQ=255)
                                                                                         cgccgcCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATCCCGGCTCCGCGCGGGATACCCGCTTCAGCCATCGCGATATCTCGTTTCAc                                                                                                                                                >  2:664313/1‑238 (MQ=255)
                                                                                                                 gATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑ATCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACGGCACCAGCAGCAGCCCATGCCCATTCCGGCTCCGCCCCGAATACCCGCTTCCGCCAACGAGATATCTCGTTTCAACCGCaaacaacacaaaaaccgacgc                                                                                                                       >  1:633677/1‑220 (MQ=255)
                                                                                                                    gTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCCTTCCGGCTCCGCGCGGCATACCCCCTTCAGCCATCGCGATATCTCGTTTCACCCGTACACCATCTAACGACTGACGCa                                                                                                                      >  2:1148632/1‑237 (MQ=255)
                                                                                                                    gTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCaaa                                                                                                                    >  2:175000/1‑239 (MQ=255)
                                                                                                                                 ggCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCtt                                                                                                       >  2:461550/1‑239 (MQ=255)
                                                                                                                                   cTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAACCCATCCCCATTCCGGCTCCGCCCGGAATACCCCCCTCCGCCATCGCGAAATCTCGTTTCACCCCTAATCCATTTAACGCCTGCCCCAAATTCTCCA‑CTCtttt                                                                                                     >  1:1043659/1‑237 (MQ=255)
                                                                                                                                   cTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTCACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGCATACCCCCTTCAGCCATCGCGATATCTCGTTTCACCCCTAACTCATCTACCGACTGTCGCACATTCTcccccccttt                                                                                                      >  2:331476/1‑230 (MQ=255)
                                                                                                                                            cccGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGCTTTCTGTCGCTGGCAACCCAACGAGACAACGTTCTGATTTCCCAGTGGCCCCTACCCATATGACCGCCAACACCCCAGCCGCCGTAGCCCATGCCCATTACGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCCCGTTTCACCCGTAACTCATCTAACGCCTGCCGCAAATTCTCCCTCTCTTTGGcacccg                                                                                              >  1:885165/1‑232 (MQ=255)
                                                                                                                                            cccGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACcccg                                                                                            >  2:1070802/1‑239 (MQ=255)
                                                                                                                                                    ttCGTAATTCCGGTTAGGGGCAAC‑CTTTGATTTCTGGCGCTGGACACCCAACGAGAAAACGTTCTGATTTCCCAGTGGCCCATCCCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAAttt                                                                                     <  2:1104878/238‑1 (MQ=255)
                                                                                                                                                      cGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCCCGCTACCTCGTTTCACCCGTAAATCATCTAACGACTACCCCCAATTCCCCAACTTTTTGGCGAccccctcacatttcc                                                                                  >  1:138207/1‑226 (MQ=255)
                                                                                                                                                           ttgcTGTTAGTGGAAAC‑CTCTGTTTTTGGTCGCTGCCACCCCAACGAAAAAACGTTCTGATTTACCAGGGGCCCATCCCCATATGCCCCCCCCCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTa                                                                             <  2:1312075/236‑1 (MQ=255)
                                                                                                                                                           cacaggTTAGTGGCAAT‑CTCTGTTTCCGGTCGCTGGCAACCCAACGATAAACCGTTCTTATTTACCAGTGGCCCATCCCAATATGCCCGCCCCCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCAACCGCGTTAATTTtgctgc                                                                               <  2:622283/232‑1 (MQ=255)
                                                                                                                                                                gTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGTTGCTATCCgg                                                                        >  2:635480/1‑239 (MQ=255)
                                                                                                                                                                   aGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGGAGCCCAGGCCCATTCCGGCTCCGCGCGGAATACCCCCTTCCGCCAACGCGATATCCCGTTTCACCCGTAACTCATCTAACGTCTCACGCCCATGCCCCCTCTCTTTTGCGACCCAGTTAATTTtgttgcgaccgggggg                                                                     >  1:672223/1‑229 (MQ=255)
                                                                                                                                                                       gCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCCCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGtt                                                                  >  2:1158439/1‑238 (MQ=255)
                                                                                                                                                                        cAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGTTcg                                                                >  2:910224/1‑239 (MQ=255)
                                                                                                                                                                           a‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCCCGGAATACCCGCCTCAGCCATCCCGATATCTCGTTTCACCCGTACCCCACCTAACGACTGCACCAAATTCTCCATCTCTTTGCCACACGCGTTAATTTTCCTCCTATCCgaggggctgtg                                                               >  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)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             
GCTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGCTCTTCTGTTGATGTCAGCGTGGCATGGCTGGAAGCGCCCGAAGGTTCGCAATTACTGTTAGTGGCAAA‑CTCTGATTTCTGTCGCTGGCAACCCAACGAGAAAACGTTCTGATTTACCAGTGGCCCATACCCATATGACCGCCACCACCGCAGCCGCCGTAGCCCATGCCCATTCCGGCTCCGCGCGGAATACCCGCTTCAGCCATCGCGATATCTCGTTTCACCCGTAACTCATCTAACGACTGACGCAAATTCTCCATCTCTTTGGCGACCGCGTTAATTTTGCTGCTATCCGGTGGGTTCGCGGCTAACAGGGCATTGTATTCATAACGCTTCGTCACCAGTTGCTGTTGCAGTGCGCTGCTT  >  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: ‑