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.

ACGTGACGGTTATCGTCAACGGTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCCTGGCTGGTGAAGTGTAC  >  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)
          tATCGTCAACGGTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACg                                                                                                                                           <  2:1652345/150‑1 (MQ=255)
            tCGTCAACGGTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTa                                                                                                                                         <  1:640727/150‑1 (MQ=255)
             cGTCAACGGTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGCCGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCCACGTAc                                                                                                                                        >  2:1233880/1‑150 (MQ=255)
             cGTCAACGGTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTAc                                                                                                                                        <  1:1252351/150‑1 (MQ=255)
              gTCAACGGTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACg                                                                                                                                       <  1:1285819/150‑1 (MQ=255)
              gTCAACGGTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACg                                                                                                                                       <  1:949420/150‑1 (MQ=255)
               tCAACGGTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTGGCAAAGTGCCGAACGTACGc                                                                                                                                      <  1:393005/150‑1 (MQ=255)
               gcAACGGTGATATCACGCACCGCGATTAGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGCACGTTTCGAAGAAGATGAAATGGGTGGGAGCGTCCGGTGCGGTAAAGAGGACTGCGTACGTTGGCAAAGTGCCGAACGTACGc                                                                                                                                      <  1:1533155/149‑1 (MQ=255)
                    ggTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCAt                                                                                                                                 >  1:297365/1‑150 (MQ=255)
                     gTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATg                                                                                                                                >  2:122098/1‑150 (MQ=255)
                     gTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTGGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATg                                                                                                                                <  2:1512068/150‑1 (MQ=255)
                       gATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGcc                                                                                                                              <  1:655504/150‑1 (MQ=255)
                           tCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGTAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACa                                                                                                                          <  1:1326650/150‑1 (MQ=255)
                            cACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAg                                                                                                                         >  2:1508723/1‑150 (MQ=255)
                            cACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAg                                                                                                                         >  1:410680/1‑150 (MQ=255)
                             aCGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTATCGAAGAACATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGt                                                                                                                        >  1:1489141/1‑150 (MQ=255)
                              cGCAATGCGATTTGCCTCGCGGCGTGTGCGCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGtt                                                                                                                       <  1:333349/150‑1 (MQ=255)
                                 actGCGATTTGCCGCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGAGCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAg                                                                                                                    <  2:1772090/148‑1 (MQ=255)
                                 aaTGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCATCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAg                                                                                                                    <  1:189449/150‑1 (MQ=255)
                                     cGATTTGCCTCGCGGCGTGTGCTCCGGATTATGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAgtgtg                                                                                                                <  2:1662810/150‑1 (MQ=255)
                                     cGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAgtgtg                                                                                                                >  1:1076422/1‑150 (MQ=255)
                                     cGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGCACGCTGCATGCCTACAGTTAAgtgtg                                                                                                                <  1:1807975/150‑1 (MQ=255)
                                           gCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCg                                                                                                          >  1:124276/1‑150 (MQ=255)
                                             cTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGt                                                                                                        <  2:124276/150‑1 (MQ=255)
                                             cTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGt                                                                                                        >  1:1302203/1‑150 (MQ=255)
                                             cTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGt                                                                                                        >  2:673185/1‑150 (MQ=255)
                                                   gcgTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCagag                                                                                                  <  1:507030/150‑1 (MQ=255)
                                                      tgtgCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGccc                                                                                               >  1:10698/1‑150 (MQ=255)
                                                         gCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCggg                                                                                            >  1:964605/1‑150 (MQ=255)
                                                            ccGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGa                                                                                         >  1:1406960/1‑150 (MQ=255)
                                                            ccGGATTAAGGGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTAAGCTGCATTCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGa                                                                                         <  1:1532354/150‑1 (MQ=255)
                                                             cGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGaa                                                                                        >  2:1385583/1‑150 (MQ=255)
                                                              ggATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAAc                                                                                       >  2:211236/1‑150 (MQ=255)
                                                                  tAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGAGCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCggg                                                                                   <  1:1310067/150‑1 (MQ=255)
                                                                   aaGTGCCGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGGCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGc                                                                                  <  1:1722937/150‑1 (MQ=255)
                                                                   aaGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGc                                                                                  >  1:204419/1‑150 (MQ=255)
                                                                      tGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCttt                                                                               >  2:1443370/1‑150 (MQ=255)
                                                                      tGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTAGGTAAAGAGGCCTGCGTACGTTCGGCAATTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCttt                                                                               >  2:303831/1‑150 (MQ=255)
                                                                        aCGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGt                                                                             >  1:1156208/1‑150 (MQ=255)
                                                                           cgcTGGAACGTTTAGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAATAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTTAATAGCCCTGGCGAACCGGGCTTTGtttt                                                                          >  2:1305206/1‑150 (MQ=255)
                                                                           cgcTGGAACGATTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCGTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGtttt                                                                          >  1:260577/1‑150 (MQ=255)
                                                                             cTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTgg                                                                        <  2:467109/150‑1 (MQ=255)
                                                                               ggAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGgt                                                                      >  2:1316288/1‑150 (MQ=255)
                                                                                 aaCGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTCCGTACGTTCCGCAATTTGCCAACGTACGCTCCATGACTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCTAACCGGGCTTTCTTTTGGgtgt                                                                    >  1:108748/1‑150 (MQ=255)
                                                                                   cGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGttt                                                                  >  1:1348149/1‑150 (MQ=255)
                                                                                      ttCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATg                                                               >  1:1473954/1‑150 (MQ=255)
                                                                                       tCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGc                                                              >  1:1472607/1‑150 (MQ=255)
                                                                                           agaagaTGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCgg                                                          >  1:1626684/1‑150 (MQ=255)
                                                                                             aagaTGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGAc                                                        >  1:1254832/1‑150 (MQ=255)
                                                                                              agaTGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACt                                                       >  2:550034/1‑150 (MQ=255)
                                                                                              agaTGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACt                                                       >  1:1750869/1‑150 (MQ=255)
                                                                                               gaTGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTa                                                      >  2:425212/1‑150 (MQ=255)
                                                                                                  gacctggtggGGAGCGTCCGGTGCGGGAAAGAGGTCGGCGTACGTGCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAgcg                                                   <  2:312125/146‑1 (MQ=255)
                                                                                                  gAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAgcg                                                   >  2:1676213/1‑150 (MQ=255)
                                                                                                  caaaagcttGGGATCGTCCGGTGCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTAATCTGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAgcg                                                   <  2:307564/143‑1 (MQ=255)
                                                                                                        gTTGGGAGCGTCCGGGTCGGTAAAGAGGACTGCGTTCGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGGGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCttttt                                             <  1:769741/150‑1 (MQ=255)
                                                                                                         ttGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTc                                            >  1:1311706/1‑150 (MQ=255)
                                                                                                         ttGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTc                                            >  1:1404128/1‑150 (MQ=255)
                                                                                                          tGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCt                                           <  2:1422351/150‑1 (MQ=255)
                                                                                                           gggATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTg                                          <  1:673185/150‑1 (MQ=255)
                                                                                                            ggATCGTCCGGTTGGGTAAAGAGGACTGCGTACGATCGGCAATTAGCCAACGTACGCTGCATGCCTACAGTAAAGTGTGTTATCGGTGAAGAGCCCGGGCGAACCGGCCTTTGTTTTGGGTGTGTATGCCCGGACTAGCGCATTTTCTGa                                         >  2:1579172/1‑150 (MQ=255)
                                                                                                             gATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGaa                                        <  1:1513241/150‑1 (MQ=255)
                                                                                                              aTCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGaaa                                       <  2:873489/150‑1 (MQ=255)
                                                                                                              aTCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGaaa                                       <  1:1677606/150‑1 (MQ=255)
                                                                                                                cGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAaaca                                     <  2:1262603/150‑1 (MQ=255)
                                                                                                                 gTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAaacac                                    >  2:1320775/1‑149 (MQ=255)
                                                                                                                  tCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAaacaac                                   >  1:208992/1‑150 (MQ=255)
                                                                                                                  tCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAaacaac                                   <  1:1676213/150‑1 (MQ=255)
                                                                                                                        tCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCAtttt                             <  1:603196/150‑1 (MQ=255)
                                                                                                                          gttAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTAAGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGAAGAGCCCGGGCGAACAGGGCTTAGTTATGGGTGTTTAAGCCCGCACTAGCGCTTTTTCTGAAACAACCATATTTa                           >  2:517703/3‑150 (MQ=255)
                                                                                                                             actgtggACTGCGTACGTTCGGCAGTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTAttt                        <  2:1489141/145‑1 (MQ=255)
                                                                                                                              aaGAGGACTGCGTACGTTCGGGACTTTGCCATCGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGTGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTg                       <  2:716656/150‑1 (MQ=255)
                                                                                                                              aaGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTg                       >  2:1072545/1‑150 (MQ=255)
                                                                                                                                gagGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGccc                    >  2:1799459/1‑150 (MQ=255)
                                                                                                                                   tgcagCGTACGTTGGCAAAGTGCCGAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCct                  >  1:643087/5‑150 (MQ=255)
                                                                                                                                   tgcagCGTACGTTGGCAAAGTGCCGAACGTACGCAGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGCACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGccccg                  >  1:437014/5‑149 (MQ=38)
                                                                                                                                   gACTGTGTAGGTTCGGCTCTTTGCCAACGTACGCTGCATTCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGTCGAACCGGGGTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCct                  <  2:309128/150‑1 (MQ=255)
                                                                                                                                    aCTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCctg                 <  2:1521366/150‑1 (MQ=255)
                                                                                                                                    aCTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCctg                 <  1:1790261/150‑1 (MQ=255)
                                                                                                                                    aCTGCGTACGTTCGGCACTTTGCCAACGTACGCGGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCctg                 <  1:158675/150‑1 (MQ=255)
                                                                                                                                    aCTGCGTACGTTCGGCACTTTGCCAACGTACGCGGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCctg                 <  1:122098/150‑1 (MQ=255)
                                                                                                                                     cTGCGTCCGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTATTCTGAAACAACCATTTTTATTTGCCCctgg                <  1:70893/150‑1 (MQ=255)
                                                                                                                                      tGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCctggc               <  2:1519480/150‑1 (MQ=255)
                                                                                                                                      tGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCctggc               >  1:1761047/1‑150 (MQ=255)
                                                                                                                                        ggTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCctggctg             <  2:123331/149‑1 (MQ=255)
                                                                                                                                         gTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCctggctgg            >  2:207177/1‑150 (MQ=255)
                                                                                                                                           aCGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCCTGGCTGGTg          <  1:842170/150‑1 (MQ=255)
                                                                                                                                           aCGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCCTGGCTGGTg          >  2:1458161/1‑150 (MQ=255)
                                                                                                                                             gTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCCTGGCTGGTGaa        <  1:1356576/150‑1 (MQ=255)
                                                                                                                                                 ggCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCCTGGCTGGTGAAgtgt    <  2:1569657/150‑1 (MQ=255)
                                                                                                                                                 ggCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCACTGGCTGGTGAAgtgt    >  1:1083646/1‑150 (MQ=255)
                                                                                                                                                   cACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCCTGGCTGGTGAAGTGTAc  <  2:10698/150‑1 (MQ=255)
                                                                                                                                                   cACTTTGCCAACGGACGCTGCAGGCCTACAGTTAAGTGTGTTAACGGGGCAGAGACCGGGCTAACCGGGCGTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCCTGGCTGGTGAAGTGTAc  <  1:322564/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACGTGACGGTTATCGTCAACGGTGATATCACGCAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCCTGGCTGGTGAAGTGTAC  >  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: