Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,293,719 A→G 17.8% E107G (GAG→GGG)  lpoA → OM lipoprotein stimulator of MrcA transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,293,7190AG17.8% 149.1 / 3.3 70E107G (GAG→GGG) lpoAOM lipoprotein stimulator of MrcA transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (16/39);  new base G (12/2);  total (29/41)
Fisher's exact test for biased strand distribution p-value = 1.72e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCCGATCAGTCCACTGCTTATATGCAGGGCACGGCGCAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTG‑GCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGA  >  NC_000913/3293492‑3293952
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              
tCCCGATCAGTCCACTGCTTATATGCAGGGCACGGCGCAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCCACTACCGCAAGAACTGGACAAGGCTCAGCGTCCCCGGAAAACACGACCGGCGGGAGGGAGTTAAATgg                                                                                                                                                                                                                                     >  1:654635/1‑239 (MQ=255)
ccccGATCAGTCCCCCTCTTAGCTTCGGGGCACGTCGCAGGCGTATTCTGTCTTTTAGCTTGAGCAGATGCCGCAAAGCCTTGATGATGCCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTGTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTgg                                                                                                                                                                                                                                     <  2:586726/238‑1 (MQ=255)
 ccccATCATTCCACTGCTTATAGGCCGGGCACGGCGCAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACGCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGgc                                                                                                                                                                                                                                    <  2:567823/239‑1 (MQ=255)
          tCCACTGCTTATATGCAGGGCACGGCGCAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTTGGTTGTTTAACCAACTACCGCAAGAACTGAACAATGCTCCGCGTCGCGGGGAAAAACTGCTGGCGGGAGGAATTAAACTGACG‑CAGaaaag                                                                                                                                                                                                                          >  1:148822/1‑237 (MQ=255)
          tCCACTGCTTATATGCAGGGCACGGCGCAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTTGGTTGGTTTACCAACACCCCGAAAAACTGAACGATGCGCAGCGGCGCGGGAAAACACCGCCGGCGGGAGGGAATAAAATGGGG‑CAaaaaaa                                                                                                                                                                                                                          >  2:1089015/1‑237 (MQ=255)
          tCCACTGCTTATATGCAGGGCACGGCGCAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGGGAAAACACTCCTGGCGGTAGGGATTAAACTGGGG‑CAGaaaaa                                                                                                                                                                                                                          >  1:331443/1‑237 (MQ=255)
                    atatGCAGGGCACGGCGCAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTTAGGTTTTTAACCAAATACACCAAAAACTGAAGCATGCTCAGCGGTGCGGGGAAAAACTCACGGGGGGAAGGATAAAACCGGGC‑CAGAAAAATTTTTTTGgc                                                                                                                                                                                                                >  2:737238/1‑239 (MQ=255)
                    atatGCAGGGCACGGCGCAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGGGAAAACACTGCTGGGGGGAGGGATTAAACTTGCG‑CCGAAAGAATTTTCTggg                                                                                                                                                                                                                >  1:25443/1‑238 (MQ=255)
                    atatGCAGGGCACGGCACAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTTGAGATTAAACTGGGG‑CCGGAAGAttttttcggc                                                                                                                                                                                                                >  1:690815/1‑233 (MQ=255)
                       tggcggggccGGCGCCGGTTGATTTTGCCTTTTTTTTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGAACAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGgcgcg                                                                                                                                                                                                             <  2:916841/230‑1 (MQ=255)
                         cAGGGCACGGCGCAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACTGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGGGAAAACACTGCTGGCGGTAGGGGTTTAAATTGGG‑CAGAAAGAtttttttgggcgcgc                                                                                                                                                                                                           >  1:568091/1‑228 (MQ=255)
                         cAGGGCACGGCGCAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGGGAAAACACTGCTGGCGGTAAGGGTTAAACTGGCG‑CGGAAAGATTTTGCTGTCGGGCa                                                                                                                                                                                                           >  1:370870/1‑239 (MQ=255)
                           gggCACGGCGCAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGTGAAAACACTGCTGGCGGCTGAGGTTAAAATGGTC‑CAGAAAAATTTTGCTGGCGCGCaaa                                                                                                                                                                                                         >  1:443639/1‑239 (MQ=255)
                           gggCACGGCGCAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCCCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CCGAAAAATTTTTCTGgcggcccaa                                                                                                                                                                                                         >  1:279224/1‑233 (MQ=255)
                              cACGGCGCAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGGTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGGAGGGATTAAACAAGCG‑GAAAAAAAGTTTTCTGGCGCGCAAAACt                                                                                                                                                                                                      >  1:766080/1‑239 (MQ=255)
                                   cgcAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTgg                                                                                                                                                                                                                                     <  1:460034/204‑1 (MQ=255)
                                   cgcAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTg                                                                                                                                                                                                                  <  1:16652/222‑1 (MQ=255)
                                       ggCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCCCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCg                                                                                                                                                                                              <  1:282901/238‑1 (MQ=255)
                                        gCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGGGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAACGTCGCGAGAAAACACTTCTGGCGGGAAGGATTAAAATGGCGCCAGAAA‑ATTTTTCTGGGCGGCCAAACTTGTTggggga                                                                                                                                                                                            >  1:508132/1‑235 (MQ=255)
                                         cTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCg                                                                                                                                                                                              <  1:198084/236‑1 (MQ=255)
                                             ttCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAACGTCGCCAGAAAACACTGCTGGCGGTAGAGATTAAACTGGGG‑ACGAAAAATTTTTCTGGCGGCCAAAACTTTCTGGCGAAAATca                                                                                                                                                                                       >  1:32230/1‑239 (MQ=255)
                                               cTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATcaca                                                                                                                                                                                     <  1:100986/239‑1 (MQ=255)
                                                    ttttttCTTCAGCAGATGCCGCAAAGCTTTGATGATACCAGGATCAACTGGCAATTCCGCGCCGTTCGTGAACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGc                                                                                                                                                                                <  2:279224/239‑1 (MQ=255)
                                                    ttttATCTTCAGCCGATGCCGCAAAGCTCTGAGGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTGGAGTTGGTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACCCCTGc                                                                                                                                                                                <  2:455065/239‑1 (MQ=255)
                                                    ttttATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGGAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGc                                                                                                                                                                                <  1:572872/239‑1 (MQ=255)
                                                    ttttATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACGGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGc                                                                                                                                                                                <  1:737238/239‑1 (MQ=255)
                                                    ttttATCGTCAGCAGATGCAGCAAAGCTCTGATGAAACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGc                                                                                                                                                                                <  1:842190/239‑1 (MQ=255)
                                                            tCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGGAGAGATTAAACTGGCG‑CCGAAAAATTTTTCTGGGCGGCAAAACTTGTTGGGGGAAAACACAccctcccgtttta                                                                                                                                                                        >  1:224223/1‑228 (MQ=255)
                                                                 agATGCAGCAAAGCTTTGATGATACCAGTATCAACTGGCAATTGCTCGCAATGCGTGCACTGGTGAAAGAAGGGAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGGAGAGATTAAACTGGGG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGGGAAAATCACACCTGCCGATTTAGa‑a                                                                                                                                                                     <  1:953423/236‑1 (MQ=255)
                                                                                      attccAGGATCAACGGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGATTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTgg                                                                                                                                            <  1:53835/236‑1 (MQ=255)
                                                                                                  aCTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGGTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTTGGG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCCCACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTAATGGCC‑AGGaaaaa                                                                                                                                 >  1:655619/1‑237 (MQ=255)
                                                                                                  aCTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGGTTAACCAACTACCGCAAGAAATGAACGATGCTCACCGTCGCGAGAAAACACTGCTGCCGGTAGGGATTAAACTGGCG‑CCGAAAGATTTTTCTGGCGCGCAAAACTTGCGGGCGCAAAACACACCCCCCGCTTTTGACACAAAACACAACAA‑C‑CGTTACCGGCA‑GGaaaaaaaa                                                                                                                               >  2:882185/1‑237 (MQ=255)
                                                                                                        aaTTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAgc                                                                                                                                                      <  1:733545/211‑1 (MQ=255)
                                                                                                           ttcTCGCCTTTCGTGCACTGGGGAAAGAGGTAAAAACCGGGCAGGCGTTTGAGTTGTTAAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATgccagc                                                                                                                      <  2:370870/237‑1 (MQ=255)
                                                                                                                      cGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGTGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGTcc                                                                                                           <  1:147627/239‑1 (MQ=255)
                                                                                                                         tccttGGTAAAACAAGGGAAACCCGGCCAGGCGTTTGGGTTTTTAAACCAATTACCGGAAGAATTGAACGATGCCCGGCGTCGCGATAAACCCCTGCTGGCGGTAGAGATTAAACTGGGG‑CAGAAGGATTTTGCGGGCGCGCAAAATGTGGTGGCGAAAATCACCCCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGTCCTTc                                                                                                        <  2:803000/235‑1 (MQ=255)
                                                                                                                         gCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGGAAGAACGGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGTCCTTc                                                                                                        <  1:494202/239‑1 (MQ=255)
                                                                                                                           acgggTAAAAGATGGTAAACCCGGGAAGGGGTTGGTGTTTTTTACCCATCTACCGCAAGACCTGAACGATGCTCAGCTTGGCGAGAAAACCCGGCTGGCGTTAGAGATTAAACTGGCG‑CAGAAAGATTTGGCTGGCGCGCAAAACTTGCTGGCGAAAATCCCACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCa                                                                                                      <  2:655619/236‑1 (MQ=255)
                                                                                                                           aCTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCa                                                                                                      <  1:271078/239‑1 (MQ=255)
                                                                                                                               gtaaaaGAGGGTAAAACCGGGCAGGCGTTTGAGTGTTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCCCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGTcc                                                                                                           <  1:392904/227‑1 (MQ=255)
                                                                                                                               gTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCCAAACTTGCTGGCGCAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGGAAAAATACATGCCAGCCAAGGGCGcccttcccttag                                                                                                  >  1:964922/1‑237 (MQ=255)
                                                                                                                               gTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGAATTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGGTTTTGA‑ACAAAAC‑CCGCAAGCGCGTTACTTGCA‑GGGCAAAAACGGAGCCAGCCAGGGGGGGCCCTCCCTTGa                                                                                                  >  1:368601/1‑239 (MQ=255)
                                                                                                                                        ggTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAATTACCGCAAGAACTGAACGATGCTCAGCTACGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTGGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATg                                                                                                                           <  1:262367/205‑1 (MQ=255)
                                                                                                                                         gTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACCCTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCa                                                                                                      <  1:168328/225‑1 (MQ=255)
                                                                                                                                         gTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTTGCGGTAGGGGTTAAAATTGCC‑CAGAAAGATTTTTCCGGCGCGCAAAACTTGCTGGGGAAAATAAACACTGACCATTTTAA‑ACAAAAC‑CAGCAAACGCGTTAATTGCAGGGCAAAAAAC‑ATACCAACCAGCGGCGGCCCTCCACTTATTTaaaggag                                                                                        >  2:687772/1‑233 (MQ=255)
                                                                                                                                          tAAACCGGGGCCGGCGTTTGGTTTGTTAACCCAACTCCCGAAAGACCTGAACGATGGTCGGGGGCGGGAGAAAACCCTGCGGGCGGTAGAGATTAACCTGGCG‑CAGAAAGATTTTGCTGGCGCGAAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTgcg                                                                                         <  2:690815/237‑1 (MQ=255)
                                                                                                                                                cGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCCCGGCGGCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTaa                                                                                 <  1:708272/239‑1 (MQ=255)
                                                                                                                                                 gggCAGGTGGTTGAGTTGTTTAACCAACTACCGCAAGGACTGAACGATGCTCAGCGTCGCGAGGAAACACTGCTGGCGGTAGCGCTTAAACTCGCG‑CAGAAAGAATTTTCTGGAGCGACCAACTTGCTGGGGAAAAACACAACGCACGATTTTTA‑ACAAAAA‑CAACAAGCGCGTTACTGGAC‑GGCAAAAATACATGCCCCCCACGGGCGCGCTTCCACTTAGTTACAGcccgcgaaaat                                                                                >  2:815333/1‑234 (MQ=255)
                                                                                                                                                 gggCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACCCTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAAt                                                                                <  1:614423/239‑1 (MQ=255)
                                                                                                                                                         ggTTGGGTTTTTTAACCAATTTCCGCAAGAACTGAACGAGGCTCGGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAgg                                                                        <  2:368601/239‑1 (MQ=255)
                                                                                                                                                                         cAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACCCTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCACGCAAAACTTGCTGGAGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTG‑gcgc                                                       <  1:312020/239‑1 (MQ=255)
                                                                                                                                                                             tACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGc                                                                                                                                                          <  1:443058/138‑1 (MQ=255)
                                                                                                                                                                             tACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGGTTTACTTCGCCCGTTAAATGCTCAGGAACCGCTGCTTT‑GCGCcaaa                                                   >  1:412680/1‑239 (MQ=255)
                                                                                                                                                                             tACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTACTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGc                                                                                                                                                          >  2:443058/1‑138 (MQ=255)
                                                                                                                                                                                   aaagaCTTAACGTTGCACAGCGTCGCGGGAAACCCTGGCTGGCGGTAGAAATTAACCTGGCC‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGGGAAAATCACACCTTCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTG‑GCGCGAAAGaaaaa                                             <  1:1085733/235‑1 (MQ=255)
                                                                                                                                                                                   aaaaaCGGAACGATGCCAAGGGTCGCGAGAAAACATTGCTGGCGGTAGGGATTAAACTGGCC‑CAGAAAGATTTTCCTGGCGCGAAAAACTTGCTGGCGAAAATCACCCCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCATGAACCGCTGCTTG‑GCGCGAAAGaaaaa                                             <  1:187281/236‑1 (MQ=255)
                                                                                                                                                                                      aaCTGAACGATGCTCAGCGTCGCGAAAAAACCCTGGTGGCGGTAGAAATTAAACTGGCG‑CAGAAATATTTTGCTGGAGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTG‑GCGCGAAAGAAAAAcag                                          <  1:270666/239‑1 (MQ=255)
                                                                                                                                                                                       aCTGAACGATGCTCATCTTCGCGATAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTGGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTG‑GCGCGAAAGAAAAAcag                                          <  1:211152/238‑1 (MQ=255)
                                                                                                                                                                                         tGAACGATACTCAGGGTCGCGAAAAAACCCTGCTGGCGGTAGAAATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACCCCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTG‑GCGCGAAAGAAAAAcagcag                                       <  1:429779/239‑1 (MQ=255)
                                                                                                                                                                                           aaCGATACTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGACGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGGCAGCACGGGGCGTCCTGCCATTGATTTACTGCGGCCGTTAATTGCTCAGGAACCGCCGCTTG‑G‑GCGAaaaaaaaaaaacaaaaa                                    >  2:429779/1‑230 (MQ=255)
                                                                                                                                                                                                  cccAGGGTCGGGAGAAAACCCTGCTGGCGGTAGAGATTAACCTGGCG‑CAGAAAGATTTTGCTGGCGGGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTG‑GCGCGAAAGAAAAACAGCAGAATATTGa                               <  1:882185/236‑1 (MQ=255)
                                                                                                                                                                                                   tCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGGCCTTCCATTGATTTACTGCGCGCGTTAATTGCCCAGGAACCGCTGCTTG‑GCGGGAAAGAAAAACAGCCGAATATTGAt                              >  1:547341/1‑238 (MQ=255)
                                                                                                                                                                                                          cgcgAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAAGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTG‑GCGCGGAAGAAAAACAGCAGAATATTGATGCCACCTg                      >  1:177022/1‑239 (MQ=255)
                                                                                                                                                                                                          cgcgAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCCCGTTACTGGCA‑GGCAAAAAACGATTCCAGCCAGGGGCGGCCTTCCCTTTATTTACTCGGGCCGGTAATTTCCCCAGAACCCCTCCTTG‑GCGCGCAAAAAAAACAGCAAAAAATTTATGCCAcccg                      >  2:1019687/1‑237 (MQ=255)
                                                                                                                                                                                                            cgAGAAAACACTGCTGGCGGTAGAGATTAAAATGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTTCCGATTTAGA‑ACACAAC‑CAGCAAGCGCGTTACCTGCC‑GGCAAAAATCAACACCAGCCAGGGGCGGACTTCCCTCGATTTACTGCGCGCGTTAATTGCCCAGCACCCGCTGCCTGTCCGCCAAGAAAAAACA‑CAGAAAATTTCTACCACCCGGc                    >  2:874352/1‑239 (MQ=255)
                                                                                                                                                                                                                    ccctgctgGCGGTAGAGATAAAACTGGCG‑CAGAAAGATTTTGCTGGCGCACAAAACTTGCTGGCGAAAATCACCCCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTG‑GCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGctc            <  1:118567/237‑1 (MQ=255)
                                                                                                                                                                                                                    ccctgcggGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAACCTTGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCAGTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTG‑GCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGctc            <  1:455777/237‑1 (MQ=255)
                                                                                                                                                                                                                    cccggctgGCGGTAGGAATTAAACTGGCG‑AAGAAAGATTTTGCTGGCGCCAAAAACTTGCTGGCGAAAACCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTG‑GCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGctc            <  1:330204/235‑1 (MQ=255)
                                                                                                                                                                                                                    cacTGATGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGAAAAACTTGCTGGCGAAAATCCCACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTG‑GCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGctc            <  1:768937/239‑1 (MQ=255)
                                                                                                                                                                                                                              ggTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCCGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTG‑GCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGGGCCTtcctccctct  >  1:336958/1‑235 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              
TCCCGATCAGTCCACTGCTTATATGCAGGGCACGGCGCAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCG‑CAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGA‑ACAAAAC‑CAGCAAGCGCGTTACTGGCA‑GGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTG‑GCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGA  >  NC_000913/3293492‑3293952

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: ‑