Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,844,263 T→G 5.2% F356V (TTC→GTC)  PP_4259 → iron‑sulfur cluster‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,844,2630TG 5.2% 367.0 / 8.3 126F356V (TTC→GTC) PP_4259iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (54/65);  new base G (4/3);  total (58/68)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 129 total aligned reads displayed.

CATCATGGACAAGATGAACTACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGATCACCTCTACGTGCTGGACGCCGCCGGC  >  NC_002947/4844115‑4844412
                                                                                                                                                    |                                                                                                                                                      
  tcatGGACAAGATGAACTACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACGTAATATTCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCg                                                                                                                                                     >  1:1332961/1‑150 (MQ=255)
  tcatGGACAAGATGAACTACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCg                                                                                                                                                     <  2:2008371/150‑1 (MQ=255)
      ggACAAGATGAACTACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCAc                                                                                                                                                >  2:2349328/1‑150 (MQ=255)
      ggACAAGATGAACTACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCAc                                                                                                                                                >  1:1393486/1‑150 (MQ=255)
      ggACAAGATGAACTACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCAc                                                                                                                                                <  2:1658620/150‑1 (MQ=255)
        aCAAGATGAACTACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCAccc                                                                                                                                              <  1:558990/150‑1 (MQ=255)
           aGATGAACTACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGtt                                                                                                                                           >  1:1551132/1‑150 (MQ=255)
            gATGAACTACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTc                                                                                                                                          >  1:355101/1‑150 (MQ=255)
             aTGAACTACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGGCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCg                                                                                                                                         <  1:1872624/150‑1 (MQ=255)
               gagcTACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCt                                                                                                                                       >  2:1174105/4‑150 (MQ=255)
               gAACTACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCt                                                                                                                                       <  1:968092/150‑1 (MQ=255)
                   tACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGTTc                                                                                                                                   <  2:355101/150‑1 (MQ=255)
                    aCCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCg                                                                                                                                  >  2:944441/1‑150 (MQ=255)
                      cccAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGAAGGCCGTCGCCAacc                                                                                                                                              <  1:26313/136‑1 (MQ=255)
                      cccAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGAAGGCCGTCGCCAacc                                                                                                                                              >  2:26313/1‑136 (MQ=255)
                       ccAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGtt                                                                                                                               <  2:317804/150‑1 (MQ=255)
                       ccAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGtt                                                                                                                               <  1:944441/150‑1 (MQ=255)
                         aaaGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTc                                                                                                                             >  2:1310164/1‑150 (MQ=255)
                         aaaGGGCTGAACAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTc                                                                                                                             >  2:91370/1‑150 (MQ=255)
                           aGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCATGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGa                                                                                                                           <  2:69931/150‑1 (MQ=255)
                             ggCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACg                                                                                                                         >  1:1583484/1‑150 (MQ=255)
                              gCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGt                                                                                                                        <  1:1174105/150‑1 (MQ=255)
                                tGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCa                                                                                                                      >  2:2146669/1‑150 (MQ=255)
                                 gATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATTCTGCGCCCCCGTCTGGTCGGCTATGCCGTGGTACTGCTGGTGATGACGATCGGGTTGGCGACGGCCTTCGCCACACGTTCGCGGGTCGGATTCGACGTccg                                                                                                                     >  2:752979/1‑148 (MQ=255)
                                 gATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAg                                                                                                                     <  1:657770/150‑1 (MQ=255)
                                    cAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCaa                                                                                                                  <  1:2141227/150‑1 (MQ=255)
                                     aGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCaaa                                                                                                                 >  2:1413061/1‑150 (MQ=255)
                                     aGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCaaa                                                                                                                 <  2:726235/150‑1 (MQ=255)
                                         acacCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGAcc                                                                                                             <  1:812917/150‑1 (MQ=255)
                                         acacCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGAcc                                                                                                             >  2:2312094/1‑150 (MQ=255)
                                             caccGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCgt                                                                                                         <  1:1734671/150‑1 (MQ=255)
                                             caccGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCgt                                                                                                         >  2:1359188/1‑150 (MQ=255)
                                              accGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCgtg                                                                                                        >  2:1612644/1‑150 (MQ=255)
                                                cGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGtgtt                                                                                                      >  1:1724768/1‑150 (MQ=255)
                                                               cGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGc                                                                                       <  1:1310164/150‑1 (MQ=255)
                                                                gggCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGcc                                                                                      <  2:1551132/150‑1 (MQ=255)
                                                                      aagaCCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGc                                                                                >  2:415122/1‑150 (MQ=255)
                                                                      aagaCCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGc                                                                                >  1:1527255/1‑150 (MQ=255)
                                                                        gaCCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCg                                                                              <  2:1332961/150‑1 (MQ=255)
                                                                        caCCCACAAGCTGCGCCCCCGTGTGCTCGGTTCCCCTGCGCTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCg                                                                              <  2:23332/149‑1 (MQ=255)
                                                                          ctcacaTGCTGCGCCCCCGTCTGATCGGCTATGCCGCGGGACGGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGa                                                                            <  1:624698/148‑1 (MQ=255)
                                                                          cccACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGa                                                                            <  2:2426879/150‑1 (MQ=255)
                                                                           ccACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTc                                                                                                                                   <  2:1157156/94‑1 (MQ=255)
                                                                           ccACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTc                                                                                                                                   >  1:1157156/1‑94 (MQ=255)
                                                                           ccACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGAt                                                                           >  1:2240522/1‑150 (MQ=255)
                                                                            cacaTGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATc                                                                          >  2:1723137/1‑150 (MQ=255)
                                                                                    gcgcCCCCGTCTGATCGGCTATGCAGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACgt                                                                  <  2:1216213/150‑1 (MQ=255)
                                                                                        ccccGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTGCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTAc                                                              >  1:1991076/1‑150 (MQ=255)
                                                                                         cccGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACa                                                             <  1:96740/150‑1 (MQ=255)
                                                                                         cccGTCGGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCCCCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACa                                                             <  1:1806122/150‑1 (MQ=255)
                                                                                          gcGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTAc                                                              >  1:684501/2‑148 (MQ=255)
                                                                                          ccGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTAc                                                              <  2:684501/148‑1 (MQ=255)
                                                                                          ccGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAg                                                            >  1:329668/1‑150 (MQ=255)
                                                                                          ccGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCAGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCaaa                                                                                                                 <  2:402816/97‑1 (MQ=255)
                                                                                            gTCTGATCGGCTATATAGGGGCACGGCTGGGGATTTTTATCGGGTTGGCGACGGCATTCGCCAACCGTTCGCTGGACGCTTTCGCCGTCAGCAAAGACCTCGTGTCGTACTGCGAATACGCCTCGGGCCGGATCGAGAACGTGTACAGcc                                                          <  2:514661/150‑1 (MQ=255)
                                                                                             tCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGtgttg                                                                                                     >  2:266049/1‑106 (MQ=255)
                                                                                             tCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGtgttg                                                                                                     <  1:266049/106‑1 (MQ=255)
                                                                                               tGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCa                                                                                                                      >  1:2068285/1‑87 (MQ=255)
                                                                                               tGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCa                                                                                                                      <  2:2068285/87‑1 (MQ=255)
                                                                                                 aTCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAg                                                     >  2:1713008/1‑150 (MQ=255)
                                                                                                   cGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACAGCGTGTAGTAGCAAGAAAACACCCATGGCCGGATCGAGAACGTGTACAGCCTGAAGGt                                                   <  1:415122/150‑1 (MQ=255)
                                                                                                    ggCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTAc                                                              >  1:1677468/1‑138 (MQ=255)
                                                                                                     gCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTaa                                                 <  1:91370/150‑1 (MQ=255)
                                                                                                     gCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTaa                                                 <  2:2193762/150‑1 (MQ=255)
                                                                                                      cTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTaa                                                 >  2:979838/1‑149 (MQ=255)
                                                                                                            ggTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTAGGTTTCGACGTCAGCAAAGACCGCGTAAAAAAAAAAAAAAACGCCCAGGGCCGGATCGAGAACGTGTCCAGCCTGAAGGGAATGAACaa                                          <  2:2368505/149‑1 (MQ=255)
                                                                                                            cGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACaa                                          <  1:412464/150‑1 (MQ=255)
                                                                                                            cGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACaa                                          <  2:1527255/150‑1 (MQ=255)
                                                                                                             gTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAAc                                                                    <  2:1652190/123‑1 (MQ=255)
                                                                                                             gTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAg                                         <  2:1583484/150‑1 (MQ=255)
                                                                                                                gtcctgcggGTGATGATGATCGGGTTGGCGACGGCCGTCGCCAACCGTTCACTGGTCGGTTTGGACGGAAGCAAACTTCGCGAGGTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGAt                                      <  2:566652/147‑1 (MQ=255)
                                                                                                                gTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGAt                                      <  1:1723137/150‑1 (MQ=255)
                                                                                                                gTACTGCTGGTGAGGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGAt                                      <  1:1036136/150‑1 (MQ=255)
                                                                                                                  aCTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGGTTCGGCGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCa                                    <  2:917571/150‑1 (MQ=255)
                                                                                                                    tgctgGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAgc                                  >  1:77985/1‑150 (MQ=255)
                                                                                                                    tgctgGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAACAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACgtg                                                                 >  1:2149160/1‑119 (MQ=255)
                                                                                                                      ctgGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAgcgc                                >  2:2413568/1‑150 (MQ=255)
                                                                                                                      ctgGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAgcgc                                >  2:2413776/1‑150 (MQ=255)
                                                                                                                       tgGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGa                                                                                                               <  2:2023322/70‑1 (MQ=255)
                                                                                                                       tgGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGa                                                                                                               >  1:2023322/1‑70 (MQ=255)
                                                                                                                       tgGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAgg                                                                                  <  1:2425810/99‑1 (MQ=255)
                                                                                                                       tgGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAgg                                                                                  <  1:1152415/99‑1 (MQ=255)
                                                                                                                         gTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGtgttgt                                                                                                    >  1:2342446/1‑79 (MQ=255)
                                                                                                                           gatgatgatCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAgc                                  <  2:402475/143‑1 (MQ=255)
                                                                                                                           gatgatgatCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAgc                                  >  1:402475/1‑143 (MQ=255)
                                                                                                                                  atCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGATCACCTCTAc                    <  1:2413776/150‑1 (MQ=255)
                                                                                                                                       gTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCCGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGTGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGATCACCTCTACGTGCt               <  1:1051194/150‑1 (MQ=255)
                                                                                                                                        ttGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGATCACCTCTACGTGCTg              <  2:1217234/150‑1 (MQ=255)
                                                                                                                                         tGGCGAAGGCCGTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTaa                                                 >  2:972044/1‑114 (MQ=255)
                                                                                                                                         tGGCGAAGGCCGTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTaa                                                 <  1:972044/114‑1 (MQ=255)
                                                                                                                                             gACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAg                                                     >  2:1015404/1‑106 (MQ=255)
                                                                                                                                             gACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAg                                                     <  1:1015404/106‑1 (MQ=255)
                                                                                                                                             gACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGATCACCTCTAc                    >  1:990866/1‑139 (MQ=255)
                                                                                                                                             gACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGATCACCTCTAc                    <  2:990866/139‑1 (MQ=255)
                                                                                                                                               aggCCGTCGCCAACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGATCACCTCTACGTGCTGGAcgcc        >  1:479256/2‑150 (MQ=255)
                                                                                                                                                   cTTCGCCACGCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCGCAGGGCCGGATc                                                                          >  2:1038211/1‑79 (MQ=255)
                                                                                                                                                   cTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCCGCAAGGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGATCACCTCTACGTGCTGGACGCCGCCgg   <  1:1341954/150‑1 (MQ=255)
                                                                                                                                                   cTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATc                                                                          <  1:1038211/79‑1 (MQ=255)
                                                                                                                                                   cTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCga                                                                        >  2:1878630/1‑81 (MQ=255)
                                                                                                                                                   cTTCGCCACCCGTTCGCTGGTCGGTTTAGACTTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCga                                                                        <  1:1878630/81‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CATCATGGACAAGATGAACTACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATGATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGATCACCTCTACGTGCTGGACGCCGCCGGC  >  NC_002947/4844115‑4844412

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