Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,378,929 A→G 19.8% L6P (CTC→CCC)  ampC ← penicillin‑binding protein, beta‑lactamase, intrinsically weak

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,378,9290AG19.8% 121.7 / 2.6 77L6P (CTC→CCC) ampCpenicillin‑binding protein, beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base A (38/20);  new base G (3/16);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 1.83e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACGGTGTAAACCACACCACAGCGGCAGAACGATCATCAGGAAC  >  NC_000913/4378697‑4379164
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         
tttGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGTTGCACAAAATCGTTGATTTGTTGGGGGGCAGCAAATGAGGAGCAAGAGGGGGTAAATAAAAAGGAGCaaagaggcg                                                                                                                                                                                                                                                     >  1:170526/1‑234 (MQ=255)
tttGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtcg                                                                                                                                                                                                                                                     >  2:190696/1‑239 (MQ=255)
 ttGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGGGCAGAGcgccgt                                                                                                                                                                                                                                                    >  2:328064/1‑239 (MQ=255)
  tGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGGGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtt                                                                                                                                                                                                                                                   >  2:795099/1‑239 (MQ=255)
       aCCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGGAATTAATAAGGCGCAGAGCGTCGTTTTGaa                                                                                                                                                                                                                                            >  2:1083592/1‑239 (MQ=255)
          gAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGGGCAAGGCGACGTTTTGaccca                                                                                                                                                                                                                                         >  2:506721/1‑235 (MQ=255)
                                      tgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAACTGTGGGGCAAGAGGCGGTAATTAATAAGGGGGCGAGCCGCGTTTTTGACAAAAGGACTGGGTTCAaaaaaaaaaaggcc                                                                                                                                                                                                           >  1:431327/1‑233 (MQ=255)
                                      tgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTTATAATGCGCAGAGGCTCGTTTTTAACATAACGTGAGGGTTCCATaaaaaagagcgc                                                                                                                                                                                                           >  1:592579/1‑233 (MQ=255)
                                      tgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGGCGGTTTGAACAAAAGGGCTGGGTTTCAAAAAAAAAGAccc                                                                                                                                                                                                           >  1:64656/1‑239 (MQ=255)
                                      tgACGGGCTGCTTTTTGCCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGAGGAGCAAGATGCGGCACTTAATAAGGGGCCGAGCGTCGTTTTTAACAAAGGGACTGGGTTCCAACaaaaaaagccc                                                                                                                                                                                                          >  1:579921/1‑232 (MQ=255)
                                              tgctTTTTGGCGCTGTACGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGAAATTAAATAAGCAGCGAGGGGCGATTTTAACATAAGGGCCGGATTTCAAAAAAAAAGGccccgcaaagg                                                                                                                                                                                                   >  1:1197378/1‑231 (MQ=255)
                                               gctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGAAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGGAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGGTTCCATACAAAACAGCCCGCCCAAGGc                                                                                                                                                                                                 >  2:1207213/1‑239 (MQ=255)
                                                ctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTTAACATAGGGCCTGGGTTCCTTACAAAACGGCCCGCCATAAGCg                                                                                                                                                                                               >  2:330624/1‑239 (MQ=255)
                                                          aTGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGGACATAGGGTCTGGTTTCTCTACAAAACGGCCCGCCATATGCGGGGCGGAgtt                                                                                                                                                                                   >  2:53537/1‑239 (MQ=255)
                                                          aTGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGAGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGCACATAGGGTCTGGGTTCCATACAAAACGGCCCGCCATAGACGGGCCCGActt                                                                                                                                                                                   >  2:397371/1‑239 (MQ=255)
                                                                                aGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGGTTTTCaattggggttggtttgattg                                                                                                                                                             >  2:313404/1‑220 (MQ=255)
                                                                                aGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAAGGTTTGGTTTCCATACAAAACCGCCCGCCCACGGCGGGGCCGGAttataatggtgaagcggagagagt                                                                                                                                                              >  1:457749/1‑217 (MQ=255)
                                                                                     taaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGGTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTTGCGATGCGGTTGATTGGAAc                                                                                                                                                         >  2:943219/1‑238 (MQ=255)
                                                                                          tttACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGGATTTAATTGGCggggccgttaattttaacaccacc                                                                                                                                                   >  2:983757/1‑216 (MQ=255)
                                                                                            tACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGGGCAGAGCGGCGGTTTGAACATAGGGGCTGGTTTCAATAAAAAAAGGCCCGCCATAGGGGGGACCGGTTTTACATTGGCGAGGCGTAAATTGAAacaacaacca                                                                                                                                                 >  1:317515/1‑235 (MQ=255)
                                                                                            tACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAGTTACATTGGGGATGGGTTAgtatttctccttttaac                                                                                                                                                 >  2:62094/1‑223 (MQ=255)
                                                                                            tACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCCGGTTTTCAATGGGGGTGCGTTAGATTTTTACGAccccca                                                                                                                                                 >  2:725658/1‑237 (MQ=255)
                                                                                            tACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCAGGGTTTTCATTTGCGATGGGTTAGAATGTAAAGAACCCaa                                                                                                                                                 >  2:268178/1‑239 (MQ=255)
                                                                                            tACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGGTTCCATACAAAACGGCCCGCCATAGGCCGGCCGGGGTTAAATTTGCGATTCGTTAGCATGTAACAacacaca                                                                                                                                                 >  2:842412/1‑236 (MQ=255)
                                                                                            tACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGGCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGGGATGCGGTTGCTTGTAAAGACACCaa                                                                                                                                                 >  2:98704/1‑239 (MQ=255)
                                                                                            tACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGGCTGGTTTCCATACAAAACGGCCCGCCAAAGGCGGGCCGGGTTTTAATTTGTGAGGCGTTAGATTTTAACggaaccca                                                                                                                                                 >  2:1157102/1‑232 (MQ=255)
                                                                                            tACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGGGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTTATAAGGCGCAGGACGGCCTTTTGGACCATAGGCCTGGGTTCCTAAAAAAAGGGCCCCCCACAGGCGGTCCCGGGTTTCATTGGTGGCACGTTAGAATGTAATAacacaca                                                                                                                                                 >  1:768859/1‑236 (MQ=255)
                                                                                                   aTAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTATTAAGGGGCGAAGCGGCGTTTTGCACAACAGGGCTGGTTTTCATAACAAACGGGCCGCCAATAGCCGGCTGGGGTTTACATTGTGGCAGCGTGAGATTTAAACAAAACCAAACAGcc                                                                                                                                           >  1:580828/1‑238 (MQ=255)
                                                                                                   aTAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACCTAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGGTTTTAATTTGGGATTCGTTTGAtctgttcttttttaaaatttc                                                                                                                                           >  2:1199385/1‑218 (MQ=255)
                                                                                                   aTAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGGTTTACATTTGGGAAGTGTTTGATTTTAACGACACCCACCAGCgg                                                                                                                                          >  2:1188182/1‑238 (MQ=255)
                                                                                                   aTAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTTCATTGTCGATGGGTTAGATTTTAACGACACCAATCACCgg                                                                                                                                          >  2:239753/1‑238 (MQ=255)
                                                                                                      actccccccccccgccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGTTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGTTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGAc                                                                                                                                       <  2:170526/228‑1 (MQ=255)
                                                                                                       gcccccccccccgccaTACCGGGGTTTTTTTGTGGCTTTATAGGGGGGGTAATGGGGAGGTGCAAAATATCGTGGAGGTGTGGGGGGGGAGAAAATGGGGAGAAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGGTTGGTTTCCATCAAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACa                                                                                                                                      <  1:943219/229‑1 (MQ=255)
                                                                                                             gccaccgccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGGTTTACATTTGCGATGCGGTAGATTTTAACGAACCCaaaccacggggcaaacgtc                                                                                                                                >  2:817480/1‑222 (MQ=255)
                                                                                                                  cgccaTACCCGGTATTTTTTGTTGCTCTAAAAGGGGGGTAATTGGGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑At                                                                                                                        <  2:1096850/239‑1 (MQ=255)
                                                                                                                             ggATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGt                                                                                                                                                                                                                                    <  2:604654/129‑1 (MQ=255)
                                                                                                                                   ttttgtttttttaaaGGGGGGAAATTGGGGGATGCAAAATTGTTTTGATTTGGTGGGGGGCAGCAAATGGGGGGCAAGGGGGGGTAATTAATAAGGGTAAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGa                                                                                                     <  1:330624/226‑1 (MQ=255)
                                                                                                                                            ctATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGGTTTACATTGGCGATGCGGTAGATTGTAACGCACCCAATAAACGTGGCAACTGTCAGG‑‑ATAGCAGCCAGACCGTCGCAAACCCAtt                                                                                            >  2:444471/1‑239 (MQ=255)
                                                                                                                                             gtgagggggtgatatttggggATGCAAAATTGTGTGGATGTGTGGGGGGGCAGAAAAGGGGGAGAGGGGGGCTGTAATTAATAAGGCGCAGGGCGTGTTTTGAAACAGAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCAttt                                                                                           <  1:444471/219‑1 (MQ=255)
                                                                                                                                                aggCGGGGTTATTGTGCAATGCAAAATATCGTTGTTTTGTTGGGGGGCAGCAAAGGTGGAGCAAGAGGCGGTAATTAATAAGGCGAAGAGCGTCGTTTTGAACATGGGGTCTGGTTTCCATAAAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGcc                                                                                        <  1:817480/237‑1 (MQ=255)
                                                                                                                                                    ggggTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTTCATTGGCGATGCGTTTGATTGTAAACACACCCATTCACGTGACCACCGTCCGGG‑‑ATACAGCCAAACCGTCACAAAACCCAttgcccccc                                                                                     >  2:1075822/1‑232 (MQ=255)
                                                                                                                                                         aaTTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAAc                                                                                                                                                                                                                                           >  2:1303435/1‑94 (MQ=255)
                                                                                                                                                         aaTTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAAc                                                                                                                                                                                                                                           <  1:1303435/94‑1 (MQ=255)
                                                                                                                                                            tgggggaggCAAAATTTCGTTTTTTTGTGGAGGGGAAGCAAAGGGGAAGCAAGTGGCGGTAATTAATAAGGCGAAGGGCGTCGTTTTGAACATAGGTTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTg                                                                           <  1:239753/231‑1 (MQ=255)
                                                                                                                                                            tgggCGATTCCCAATTTGGTTGTTTTGGTGAGGGGTAGAAATTGTGGAGAAAGGGGAGGAAATTAATAGGGCGCAGAGCGTCGTTTTGCACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTg                                                                           <  1:506721/236‑1 (MQ=255)
                                                                                                                                                             aaagagaaagaaaTATCGTTGGTTTGGTGGGGGGCCGCAAATGGGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAGAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTg                                                                           <  2:1298973/227‑1 (MQ=255)
                                                                                                                                                                cGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTTCATTGGCGGTGCGTTAGATTGTAACGACACCAATCAGCCTGACAACTGTCAGG‑‑GTAACAGGCACAACGTAGAAAAACCATTTGCCCGcccggaacgggctt                                                                        >  2:100988/1‑226 (MQ=255)
                                                                                                                                                                   tGCACAATTTCGTTTAGTTGGTGAGGGGCAGGAAAGGGGGAGCAAGGGGCGGTAATTAATAAGGCGCAGAGGGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTc                                                                    <  2:457749/239‑1 (MQ=255)
                                                                                                                                                                   tGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGt                                                                                                                                                                    <  2:645454/149‑1 (MQ=255)
                                                                                                                                                                      acaATTTGGTTGATTGTTTGGGGGGAAGAAAATGTGGAGCAAGGGGCGTTAATTAATAAGGCGGAGGGCGTGGTTTTGAAAATGGGGTTGGGTTTCCATAAAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGa                                                                 <  2:592579/239‑1 (MQ=255)
                                                                                                                                                                          tatCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTaa                                                                                                                                                          <  2:430952/152‑1 (MQ=255)
                                                                                                                                                                             cGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAAGCGCAGAGCGTCGGTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAACCGGGCCGGGTTTAAATTTGCGATACGTTAAATTGTAACGAAACCACACACCGGGAAAAACGTTAGG‑‑GAAGTCAGCCGAACCCACAAAAACCCATTTGCCCCCAGGCCGCGGAtttttacaacagggc                                                          >  1:166/1‑228 (MQ=255)
                                                                                                                                                                               ttGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGGTTTACATTTGCGATGCGTTAGATTGTAACGACCCCAATCAACGTGACAACTGTCAGG‑‑ATAGAAGCCCAACCGTAGGAAACCCCATTTCCCCCACGGACCGGGGTTTTCaaaacggggca                                                         >  2:766426/1‑228 (MQ=255)
                                                                                                                                                                                 tttttGTTGGGGGGTAGAAAATGGGGAGCATGGGGGGGTAATTAATAGGGCGCAGGGCGTCGTTTTGAACAGGGGGTCTGGTTTCCATACAAACCGGGCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACCCCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATcgc                                                      <  1:1157102/237‑1 (MQ=255)
                                                                                                                                                                                 gattgttgggggggCAGAAAATGGGGACCAAGGGGCTGTAATTAATAAGGCGCAGGGCGTGGTTTTGAACATAGGGTCTGGTTTCCATAAAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATcgc                                                      <  2:1217535/229‑1 (MQ=255)
                                                                                                                                                                                  aTTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGGTTTACATTGGCGATGCGGTAGATTGTAACGACACCAATCAGCGTGGCAACTGTCAGG‑‑GTAGCAGCCAGACCGTAGAAAACCCATTTTCCCGCAGGTACGTGGAGTTTCCGATTGTGCATcccg                                                     >  2:371896/1‑239 (MQ=255)
                                                                                                                                                                                   tttttggggggtgattaaaTGTGGGGCGAGGGGGGTTAATTAAGAAGGCGCAGGGCTGCTTTTAGACCATGGGGTGTGTTTTCCATACAAAAGGGCCCGCCCTAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGGTTTTCAGATCGTGCATCGCgt                                                    <  1:313404/223‑1 (MQ=255)
                                                                                                                                                                                   tttGTTGGGGGGCAACAAATGTGGAGCAAGGGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgt                                                    <  2:1197378/239‑1 (MQ=255)
                                                                                                                                                                                   tttGTTGGGGGGAAGCAAATGGGGAGCAAGAGGTGGTAATTAATAAGGCGAAGAGCTACGTTTGGAAAATGGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgt                                                    <  2:837856/239‑1 (MQ=255)
                                                                                                                                                                                   tttGTTGGGGGGAAGAAAATGGGGAGAAAGAGGCGGTAATTAAAAAGGCGCAGGGCGTCGTTTTAAACATAGGGTTTGGTTTCCAAAAAAACCGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgt                                                    <  2:625492/239‑1 (MQ=255)
                                                                                                                                                                                   tttGTTGAGGGGACGCAAATGTGGAGCAAGGGGCGGAAATTAATAAGGCGCAGAGCGTCGTTTTAAACATAGGGTCTGGTTTCCATACAAACCGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgt                                                    <  2:1307955/239‑1 (MQ=255)
                                                                                                                                                                                   ttggttgGGGGGCAGCAAAGGTGGAGCAAGAGGCGGTAATAAATAAGGCGCAGAGCGTGTTTTGGAACATAGGGTCTGGTTTCCATACAAACCGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgt                                                    <  2:579921/236‑1 (MQ=255)
                                                                                                                                                                                   ttggtggAGGGGCAGCAAAGGTGAATCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGGTGGAAACATAGGTTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACAGTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgt                                                    <  2:527104/233‑1 (MQ=255)
                                                                                                                                                                                    ttgttgAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGGGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGGTTTACATTGGGCGTGCGGTAGATTGTAACAGCCCCAATACGCGGGAAAACCGTGAGGGCAAAGCACCAACCCAGAAAAAACCCATTTCCCCCACGGAACGGGGAttttaaaaaggcgccacggcg                                                    >  1:468511/1‑222 (MQ=255)
                                                                                                                                                                                    ggggaaggtgtaaaagggaagggAGAAAGGGGCGGTAATTAAAAAGGCCGAGGGGTTCTTTTTGAACATAGGGTCGGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgt                                                    <  2:186950/217‑1 (MQ=255)
                                                                                                                                                                                            gggCAGCAATTGGGAAGCAGGGGGCTGTAATAAATAAGGCGCAGGGCGGCGTTTTAAAAATGGGTTCTGGTTTCCAAACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATa                                            <  1:190696/238‑1 (MQ=255)
                                                                                                                                                                                            gggCAGAAATGGGGGAGCAGGGGGGGTAAATGAATAAGGCGGAGGGCGTCTTTTGTAGCATGGGGTCGGGTTGCCAACCAAACGGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAATCCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATa                                            <  1:842412/238‑1 (MQ=255)
                                                                                                                                                                                            gggCAGAAAATGGGGAGAAAGGGGGGGTAATTAATAAGGCGCAGGGGGTCTTTTTGAACATAGGGTCTGGTTGCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATa                                            <  2:1108469/238‑1 (MQ=255)
                                                                                                                                                                                               gccgaaaagtggggagagtgggggTAAATAAAAAAGGCCAGAGGGTGCTTTTGAAGATAGGGTTTGGGTTTCCAAAAAAAAGGGCCCCCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACg                                          <  1:1188182/216‑1 (MQ=255)
                                                                                                                                                                                                    aaTGGGGAGCAAGGGGCGGTAATTAAAAGGGCCCAGGGGGTCGTTTTGAACATAGGGTCTGGTTTCAATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACGGTGTaaa                                   <  2:159116/239‑1 (MQ=255)
                                                                                                                                                                                                         ggAGCAAGAGGCGGTAATTAATAAGGCGCAGGGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACGGTGTAAaccac                               <  2:967736/238‑1 (MQ=255)
                                                                                                                                                                                                               aggggCGGTAATTAATAAGGCGGAGGGGGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAACCGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACGGTGTAAACCACACCACAg                        <  2:318838/236‑1 (MQ=255)
                                                                                                                                                                                                                     ggTAATTAATGAGGCGCGGAGTGTCTTTTTGAACATAGGTTCTGTGTTCCACTACAAACGGCCGGCCTTAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACGGTGTAAACCACACCACAGCGGCa                   <  1:328064/238‑1 (MQ=255)
                                                                                                                                                                                                                       taaataatAAGGCGCAGGGCGTCGTTTTGAACATAGGGGTTGGTTTGCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACGgtg                                       <  2:480118/216‑1 (MQ=255)
                                                                                                                                                                                                                                aaGGCGGCGGGCGTCGTTTGAAAAATGGTGTTGGGTTTCCTAACAAAGCGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACGGTGTAAACCACACCACAGCGGCAGAACGatcatc        <  1:1075822/239‑1 (MQ=255)
                                                                                                                                                                                                                                aaGGCCCAGAGCTTCTTTTTAAAAAAGGGGTTTGGTTTCCATAAAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACGGTGTAAACCACACCACAGCGGCAGAACGatcatc        <  2:1106151/239‑1 (MQ=255)
                                                                                                                                                                                                                                      cAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGCTGCGTTAGATTTGAACGACACCAACCCGCGTGACAACCGTCAGG‑‑ATAGCAGCCAGGACGAAGAAAAACCACTTGCCCCCCGGTACCGTGATTTTTCAGACCGGCAACCCGCGGTCCCTACGGTGGAAAccacaccacacaggccaaaaagacacacagaac  >  1:68664/1‑216 (MQ=255)
                                                                                                                                                                                                                                      cAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGGACGTGGATTTTCCGATCGTGCATCGCGTGGTGCATACGGTGTAAACCACACCACAACGGCAGAACGATCATcaggaaa  >  2:174105/1‑238 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         
TTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGG‑‑ATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACGGTGTAAACCACACCACAGCGGCAGAACGATCATCAGGAAC  >  NC_000913/4378697‑4379164

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