Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,086,024 C→A 5.2% L345I (CTT→ATT)  phoH → ATP‑binding protein; putative PhoH family P‑loop ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,086,0240CA 5.2% 375.1 / 3.7 122L345I (CTT→ATT) phoHATP‑binding protein; putative PhoH family P‑loop ATPase
Reads supporting (aligned to +/- strand):  ref base C (55/58);  new base A (6/2);  total (61/61)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 123 total aligned reads displayed.

ACGTGACGGTTATCGTCAACGGTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCCTGGCTGGTGAAGTGTAC  >  NC_000913/1085875‑1086171
                                                                                                                                                     |                                                                                                                                                    
cgttGACGATAACCGTCAACGGTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCAc                                                                                                                                                      >  1:1569657/4‑150 (MQ=255)
   tGACGGTTATCGTCAACGGTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACttt                                                                                                                                                   >  2:685650/1‑150 (MQ=255)
   tGACGGTTATCGTCAACGGTGATATCACGCAATGCGATTTGCCTCGCGGCGGGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACttt                                                                                                                                                   <  2:651681/150‑1 (MQ=255)
    gACGGTTATCGTCAACGGTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTg                                                                                                                                                  <  1:1537825/150‑1 (MQ=255)
     aCGGTTATCGTCAACGGTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGc                                                                                                                                                 <  1:1456741/150‑1 (MQ=255)
      cGGTTATCGTCAACGGTGATATCACGCAATGCGATTTGCCTCGCTGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATTGTTGGGATCGTCCGGTTCGGTAAAGAGTACTGCGTACGTTCGGCACTTTGcc                                                                                                                                                >  1:309128/1‑150 (MQ=255)
      cGGTTATCGTCAACGGTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGcc                                                                                                                                                >  1:873489/1‑150 (MQ=255)
          tATCGTCAACGGTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACg                                                                                                                                           <  2:1652345/150‑1 (MQ=255)
            tCGTCAACGGTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTa                                                                                                                                         <  1:640727/150‑1 (MQ=255)
             cGTCAACGGTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGCCGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑CACGTAc                                                                                                                                        >  2:1233880/1‑150 (MQ=255)
             cGTCAACGGTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTAc                                                                                                                                        <  1:1252351/150‑1 (MQ=255)
              gTCAACGGTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACg                                                                                                                                       <  1:1285819/150‑1 (MQ=255)
              gTCAACGGTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACg                                                                                                                                       <  1:949420/150‑1 (MQ=255)
               tCAACGGTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTT‑GGCAAAGTGCCGAACGTACGc                                                                                                                                      <  1:393005/150‑1 (MQ=255)
               gcAACGGTGATATCACGCACCGCGATTAGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGCACGTTTCGAAGAAGATGAAATGGGTGGGAGCGTCCGGTGCGGTAAAGAGGACTGCGTACGTT‑GGCAAAGTGCCGAACGTACGc                                                                                                                                      <  1:1533155/149‑1 (MQ=255)
                    ggTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCAt                                                                                                                                 >  1:297365/1‑150 (MQ=255)
                     gTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATg                                                                                                                                >  2:122098/1‑150 (MQ=255)
                     gTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTGGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATg                                                                                                                                <  2:1512068/150‑1 (MQ=255)
                       gATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGcc                                                                                                                              <  1:655504/150‑1 (MQ=255)
                           tCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGTAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACa                                                                                                                          <  1:1326650/150‑1 (MQ=255)
                            cACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAg                                                                                                                         >  2:1508723/1‑150 (MQ=255)
                            cACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAg                                                                                                                         >  1:410680/1‑150 (MQ=255)
                             aCGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTATCGAAGAACATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGt                                                                                                                        >  1:1489141/1‑150 (MQ=255)
                              cGCAATGCGATTTGCCTCGCGGCGTGTGCGCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGtt                                                                                                                       <  1:333349/150‑1 (MQ=255)
                                 actGCGATTTGCCGCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGAGCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAg                                                                                                                    <  2:1772090/148‑1 (MQ=255)
                                 aaTGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCATCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAg                                                                                                                    <  1:189449/150‑1 (MQ=255)
                                     cGATTTGCCTCGCGGCGTGTGCTCCGGATTATGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAgtgtg                                                                                                                <  2:1662810/150‑1 (MQ=255)
                                     cGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAgtgtg                                                                                                                >  1:1076422/1‑150 (MQ=255)
                                     cGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGCACGCTGCATGCCTACAGTTAAgtgtg                                                                                                                <  1:1807975/150‑1 (MQ=255)
                                           gCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCg                                                                                                          >  1:124276/1‑150 (MQ=255)
                                             cTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGt                                                                                                        <  2:124276/150‑1 (MQ=255)
                                             cTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGt                                                                                                        >  1:1302203/1‑150 (MQ=255)
                                             cTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGt                                                                                                        >  2:673185/1‑150 (MQ=255)
                                                   gcgTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCagag                                                                                                  <  1:507030/150‑1 (MQ=255)
                                                      tgtgCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGccc                                                                                               >  1:10698/1‑150 (MQ=255)
                                                         gCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCggg                                                                                            >  1:964605/1‑150 (MQ=255)
                                                            ccGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGa                                                                                         >  1:1406960/1‑150 (MQ=255)
                                                            ccGGATTAAGGGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTAAGCTGCATTCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGa                                                                                         <  1:1532354/150‑1 (MQ=255)
                                                             cGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGaa                                                                                        >  2:1385583/1‑150 (MQ=255)
                                                              ggATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAAc                                                                                       >  2:211236/1‑150 (MQ=255)
                                                                  tAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGAGCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCggg                                                                                   <  1:1310067/150‑1 (MQ=255)
                                                                   aaGTGCCGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGGCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGc                                                                                  <  1:1722937/150‑1 (MQ=255)
                                                                   aaGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGc                                                                                  >  1:204419/1‑150 (MQ=255)
                                                                      tGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCttt                                                                               >  2:1443370/1‑150 (MQ=255)
                                                                      tGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTAGGTAAAGAGGCCTGCGTACGTTCGGCAATTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCttt                                                                               >  2:303831/1‑150 (MQ=255)
                                                                        aCGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGt                                                                             >  1:1156208/1‑150 (MQ=255)
                                                                           cgcTGGAACGTTTAGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAATAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTTAATAGCCCTGGCGAACCGGGCTTTGtttt                                                                          >  2:1305206/1‑150 (MQ=255)
                                                                           cgcTGGAACGATTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCGTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGtttt                                                                          >  1:260577/1‑150 (MQ=255)
                                                                             cTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTgg                                                                        <  2:467109/150‑1 (MQ=255)
                                                                               ggAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGgt                                                                      >  2:1316288/1‑150 (MQ=255)
                                                                                 aaCGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTCCGTACGTTCCGCAATTTGCC‑AACGTACGCTCCATGACTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCTAACCGGGCTTTCTTTTGGgtgt                                                                    >  1:108748/1‑150 (MQ=255)
                                                                                   cGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGttt                                                                  >  1:1348149/1‑150 (MQ=255)
                                                                                      ttCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATg                                                               >  1:1473954/1‑150 (MQ=255)
                                                                                       tCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGc                                                              >  1:1472607/1‑150 (MQ=255)
                                                                                           agaagaTGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCgg                                                          >  1:1626684/1‑150 (MQ=255)
                                                                                             aagaTGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGAc                                                        >  1:1254832/1‑150 (MQ=255)
                                                                                              agaTGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACt                                                       >  2:550034/1‑150 (MQ=255)
                                                                                              agaTGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACt                                                       >  1:1750869/1‑150 (MQ=255)
                                                                                               gaTGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTa                                                      >  2:425212/1‑150 (MQ=255)
                                                                                                  gacctggtggGGAGCGTCCGGTGCGGGAAAGAGGTCGGCGTACGTGCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAgcg                                                   <  2:312125/146‑1 (MQ=255)
                                                                                                  gAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAgcg                                                   >  2:1676213/1‑150 (MQ=255)
                                                                                                  caaaagcttGGGATCGTCCGGTGCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTAATCTGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAgcg                                                   <  2:307564/143‑1 (MQ=255)
                                                                                                        gTTGGGAGCGTCCGGGTCGGTAAAGAGGACTGCGTTCGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGGGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCttttt                                             <  1:769741/150‑1 (MQ=255)
                                                                                                         ttGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTc                                            >  1:1311706/1‑150 (MQ=255)
                                                                                                         ttGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTc                                            >  1:1404128/1‑150 (MQ=255)
                                                                                                          tGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCt                                           <  2:1422351/150‑1 (MQ=255)
                                                                                                           gggATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTg                                          <  1:673185/150‑1 (MQ=255)
                                                                                                            ggATCGTCCGGTTGGGTAAAGAGGACTGCGTACGATCGGCAATTAGCC‑AACGTACGCTGCATGCCTACAGTAAAGTGTGTTATCGGTGAAGAGCCCGGGCGAACCGGCCTTTGTTTTGGGTGTGTATGCCCGGACTAGCGCATTTTCTGa                                         >  2:1579172/1‑150 (MQ=255)
                                                                                                             gATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGaa                                        <  1:1513241/150‑1 (MQ=255)
                                                                                                              aTCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGaaa                                       <  2:873489/150‑1 (MQ=255)
                                                                                                              aTCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGaaa                                       <  1:1677606/150‑1 (MQ=255)
                                                                                                                cGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAaaca                                     <  2:1262603/150‑1 (MQ=255)
                                                                                                                 gTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAaacac                                    >  2:1320775/1‑149 (MQ=255)
                                                                                                                  tCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAaacaac                                   >  1:208992/1‑150 (MQ=255)
                                                                                                                  tCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAaacaac                                   <  1:1676213/150‑1 (MQ=255)
                                                                                                                        tCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCAtttt                             <  1:603196/150‑1 (MQ=255)
                                                                                                                          gttAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTAAGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGAAGAGCCCGGGCGAACAGGGCTTAGTTATGGGTGTTTAAGCCCGCACTAGCGCTTTTTCTGAAACAACCATATTTa                           >  2:517703/3‑150 (MQ=255)
                                                                                                                             actgtggACTGCGTACGTTCGGCAGTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTAttt                        <  2:1489141/145‑1 (MQ=255)
                                                                                                                              aaGAGGACTGCGTACGTTCGGGACTTTGCC‑ATCGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGTGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTg                       <  2:716656/150‑1 (MQ=255)
                                                                                                                              aaGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTg                       >  2:1072545/1‑150 (MQ=255)
                                                                                                                                gagGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTG‑TTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGccc                    >  2:1799459/1‑150 (MQ=255)
                                                                                                                                   tgcagCGTACGTT‑GGCAAAGTGCCGAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCct                  >  1:643087/5‑150 (MQ=255)
                                                                                                                                   tgcagCGTACGTT‑GGCAAAGTGCCGAACGTACGCAGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGCACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGccccg                  >  1:437014/5‑149 (MQ=38)
                                                                                                                                   gACTGTGTAGGTTCGGCTCTTTGCC‑AACGTACGCTGCATTCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGTCGAACCGGGGTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCct                  <  2:309128/150‑1 (MQ=255)
                                                                                                                                    aCTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCctg                 <  2:1521366/150‑1 (MQ=255)
                                                                                                                                    aCTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCctg                 <  1:1790261/150‑1 (MQ=255)
                                                                                                                                    aCTGCGTACGTTCGGCACTTTGCC‑AACGTACGCGGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCctg                 <  1:158675/150‑1 (MQ=255)
                                                                                                                                    aCTGCGTACGTTCGGCACTTTGCC‑AACGTACGCGGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCctg                 <  1:122098/150‑1 (MQ=255)
                                                                                                                                     cTGCGTCCGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTATTCTGAAACAACCATTTTTATTTGCCCctgg                <  1:70893/150‑1 (MQ=255)
                                                                                                                                      tGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCctggc               <  2:1519480/150‑1 (MQ=255)
                                                                                                                                      tGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCctggc               >  1:1761047/1‑150 (MQ=255)
                                                                                                                                        ggTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCctggctg             <  2:123331/149‑1 (MQ=255)
                                                                                                                                         gTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCctggctgg            >  2:207177/1‑150 (MQ=255)
                                                                                                                                           aCGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCCTGGCTGGTg          <  1:842170/150‑1 (MQ=255)
                                                                                                                                           aCGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCCTGGCTGGTg          >  2:1458161/1‑150 (MQ=255)
                                                                                                                                             gTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCCTGGCTGGTGaa        <  1:1356576/150‑1 (MQ=255)
                                                                                                                                                 ggCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCCTGGCTGGTGAAgtgt    <  2:1569657/150‑1 (MQ=255)
                                                                                                                                                 ggCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCACTGGCTGGTGAAgtgt    >  1:1083646/1‑150 (MQ=255)
                                                                                                                                                   cACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCCTGGCTGGTGAAGTGTAc  <  2:10698/150‑1 (MQ=255)
                                                                                                                                                   cACTTTGCC‑AACGGACGCTGCAGGCCTACAGTTAAGTGTGTTAACGGGGCAGAGACCGGGCTAACCGGGCGTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCCTGGCTGGTGAAGTGTAc  <  1:322564/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACGTGACGGTTATCGTCAACGGTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCC‑AACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCCTGGCTGGTGAAGTGTAC  >  NC_000913/1085875‑1086171

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑