Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,843,922 T→G 9.0% I242S (ATC→AGC)  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,843,9220TG 9.0% 311.2 / 20.4 113I242S (ATC→AGC) PP_4259iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (49/53);  new base G (5/5);  total (54/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

Only 100 of 114 total aligned reads displayed.

GGCCAGGCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGCCTGCAGA  >  NC_002947/4843774‑4844071
                                                                                                                                                    |                                                                                                                                                     
ggCCAGGCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATc                                                                                                                                                      <  1:1183964/150‑1 (MQ=255)
gcccAGGCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATc                                                                                                                                                      >  1:84831/3‑150 (MQ=255)
     ggCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTc                                                                                                                                                 >  2:462485/1‑150 (MQ=255)
        cGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTa                                                                                                                                              >  2:1703617/1‑150 (MQ=255)
        cGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTa                                                                                                                                              >  1:253945/1‑150 (MQ=255)
        cGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGATCAGCGTGTCCtt                                                                                                                                              <  1:1100346/150‑2 (MQ=255)
        cGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGATCAGCGTGTCCtt                                                                                                                                              <  1:415365/150‑2 (MQ=255)
           tGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGa                                                                                                                                           <  2:2166264/150‑1 (MQ=255)
            ggCGGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGAc                                                                                                                                          <  2:551877/150‑1 (MQ=255)
             gCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCgt                                                                                                                                                    >  1:2095933/1‑139 (MQ=255)
             gCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCgt                                                                                                                                                    <  2:2095933/139‑1 (MQ=255)
             gCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGAcc                                                                                                                                         <  2:2187107/150‑1 (MQ=255)
                   ccTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGcc                                                                                                                                   >  1:249032/1‑150 (MQ=255)
                     tACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGt                                                                                                                                 <  1:1402986/150‑1 (MQ=255)
                       cTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTgg                                                                                                                               <  1:1105601/150‑1 (MQ=255)
                           tGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGaa                                                                                                                           >  2:840438/1‑150 (MQ=255)
                             ggTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTTGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAAc                                                                                                                         <  2:230953/150‑1 (MQ=255)
                              gTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCgtg                                                                                                                                                   >  1:1110306/1‑123 (MQ=255)
                              gTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCgtg                                                                                                                                                   <  2:1110306/123‑1 (MQ=255)
                              gTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAAcc                                                                                                                        <  1:462485/150‑1 (MQ=255)
                                 gtcttcttCACCCTCGCCACTTCCGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTAGGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCcgc                                                                                                                     <  2:1966209/148‑1 (MQ=255)
                                  gCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCcgcg                                                                                                                    <  1:847567/150‑1 (MQ=255)
                                  gCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCcgcg                                                                                                                    <  1:1703617/150‑1 (MQ=255)
                                    ttcttcACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGc                                                                                                                  >  1:1648838/1‑150 (MQ=255)
                                    ttcttcACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGATCAGCGTGTCCtt                                                                                                                                              <  1:2276315/122‑2 (MQ=255)
                                    ttcttcACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGATCAGCGTGTCCtt                                                                                                                                              >  2:2276315/1‑121 (MQ=255)
                                           cccTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCa                                                                                                           >  2:1407291/1‑150 (MQ=255)
                                            ccTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCaa                                                                                                          >  2:1643279/1‑150 (MQ=255)
                                            ccTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGATCAGCGTGTCCtt                                                                                                                                              >  1:1486697/1‑113 (MQ=37)
                                            ccTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGATCAGCGTGTCCtt                                                                                                                                              <  2:1486697/114‑2 (MQ=37)
                                             cTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCaag                                                                                                         >  1:1579596/1‑150 (MQ=255)
                                                gCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCaagaaa                                                                                                      <  1:1824729/150‑1 (MQ=255)
                                                    cTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATa                                                                                                  <  2:253945/150‑1 (MQ=255)
                                                    cTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATa                                                                                                  <  2:1991678/150‑1 (MQ=255)
                                                     ttACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGATCAGCGTGTCCtt                                                                                                                                              <  2:452838/105‑2 (MQ=37)
                                                     ttACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGATCAGCGTGTCCtt                                                                                                                                              >  1:452838/1‑104 (MQ=37)
                                                        cGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGa                                                                                              <  1:1423269/150‑1 (MQ=255)
                                                         ggCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGTGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGAc                                                                                             <  1:670520/150‑1 (MQ=255)
                                                         ggCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTgg                                                                                                                               <  1:1851247/116‑1 (MQ=255)
                                                         ggCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTgg                                                                                                                               >  2:1851247/1‑116 (MQ=255)
                                                          gCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCaagaaa                                                                                                      >  2:73212/1‑140 (MQ=255)
                                                                 cGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAg                                                                                       <  1:1345039/148‑1 (MQ=255)
                                                                 cGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAg                                                                                       >  2:1345039/1‑148 (MQ=255)
                                                                     tggctACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACCCATGGCACGTGTCCTTTTGCCGATCGCGTTGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAg                                                                                 <  1:1499199/150‑1 (MQ=25)
                                                                     tggctACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTAATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAg                                                                                 >  2:718274/1‑150 (MQ=255)
                                                                      ggctACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAATAAAGATACCGACTACAAGGCCCAgg                                                                                <  1:1643279/150‑1 (MQ=255)
                                                                         tACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGcc                                                                             >  1:1793010/1‑150 (MQ=255)
                                                                         tACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTAATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGcc                                                                             <  1:718274/150‑1 (MQ=255)
                                                                           cgcgAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATATTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGTCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTg                                                                           <  1:2300346/150‑1 (MQ=255)
                                                                            gcgAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTgg                                                                          <  2:1793010/150‑1 (MQ=255)
                                                                            gcgAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTgg                                                                          <  1:2230297/150‑1 (MQ=255)
                                                                              gAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGt                                                                        <  2:1579596/150‑1 (MQ=255)
                                                                              gAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGt                                                                        <  1:840438/150‑1 (MQ=255)
                                                                               aaCAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTgg                                                                          >  2:1284613/1‑147 (MQ=255)
                                                                               aaCAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTgg                                                                          <  1:1284613/147‑1 (MQ=255)
                                                                                    gtgtGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGc                                                                  >  2:1671617/1‑150 (MQ=255)
                                                                                         cGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGa                                                             <  1:111451/150‑1 (MQ=255)
                                                                                          gtgtACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGAc                                                                                             <  2:2324226/117‑1 (MQ=255)
                                                                                          gtgtACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGAc                                                                                             >  1:2324226/1‑117 (MQ=255)
                                                                                          gtgtACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGAc                                                            <  1:1676321/150‑1 (MQ=255)
                                                                                          gtgtACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGAc                                                            >  1:18697/1‑150 (MQ=255)
                                                                                           tgtACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACt                                                           <  1:40532/150‑1 (MQ=255)
                                                                                            gtACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTg                                                          >  1:1972221/1‑150 (MQ=255)
                                                                                              aCATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAACGACCCGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCa                                                        >  1:2225816/1‑150 (MQ=255)
                                                                                                  gtgCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCAt                                                    >  2:1905943/1‑150 (MQ=255)
                                                                                                    gCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCAtgc                                                  <  2:1648838/150‑2 (MQ=255)
                                                                                                      cccTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCAcgcacatg                                                >  2:1033342/1‑143 (MQ=255)
                                                                                                      cccTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATgtgc                                                >  1:1628661/1‑150 (MQ=255)
                                                                                                        cTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATgtgcgt                                              >  2:1151109/1‑150 (MQ=255)
                                                                                                          aTGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATgtgcgtgc                                            <  2:854907/150‑1 (MQ=255)
                                                                                                           tGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCACGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCa                                           >  1:1833815/1‑150 (MQ=255)
                                                                                                           tGCCCGCTTCCAGAGCGTAATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCa                                           <  2:2225816/150‑1 (MQ=255)
                                                                                                            gCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAg                                          <  2:1972221/150‑1 (MQ=255)
                                                                                                              ccGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAg                                                                                 >  1:638330/1‑109 (MQ=255)
                                                                                                               cGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGGGGCGAAACCCGCGGCTCGCGCAAGAAAAAGACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATAGACTGCACCATGTGCGTGCAGGTc                                       >  1:274983/1‑150 (MQ=37)
                                                                                                                  ttCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGc                                    >  1:357045/1‑150 (MQ=255)
                                                                                                                    ccAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGccc                                  <  2:18697/150‑1 (MQ=255)
                                                                                                                        agCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGAcc                              <  2:357045/150‑1 (MQ=255)
                                                                                                                            tGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACTGACAACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCa                          >  1:1957183/1‑150 (MQ=255)
                                                                                                                             gATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCAt                         <  2:1628661/150‑1 (MQ=255)
                                                                                                                               tGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATa                                                                                                  <  1:1543989/75‑1 (MQ=255)
                                                                                                                               tGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATa                                                                                                  >  2:1543989/1‑75 (MQ=255)
                                                                                                                               tGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAg                                          <  2:508140/131‑1 (MQ=255)
                                                                                                                               tGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAg                                          >  1:508140/1‑131 (MQ=255)
                                                                                                                                gTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCaa                                                                                                          >  1:867265/1‑66 (MQ=255)
                                                                                                                                gTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCaa                                                                                                          <  2:867265/66‑1 (MQ=255)
                                                                                                                                  tCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACa                    >  1:1948859/1‑150 (MQ=255)
                                                                                                                                    gACAAGGACACGCTGATCCTATCCTACGACCCACCCCGTGGCGAAACCCAGGGCCCGCGCAACAAATATACCGACTACAAGGCCCAGGGCCTGGGCGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACaaa                  <  1:953389/150‑3 (MQ=255)
                                                                                                                                     cgtaGGACACGATCAGCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATcc                 >  1:663722/4‑150 (MQ=255)
                                                                                                                                       taGGACACGATCAGCGTGTCCTACGACCCACGCCGTGGAGAAACCCGCGGCCGGCGCACGAAAGATACCGAATACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCACGACTGGCATCGACATCCGt               >  2:1076769/2‑150 (MQ=255)
                                                                                                                                       aaGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGt               >  2:1947312/1‑150 (MQ=255)
                                                                                                                                            cacGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTg                                                                           >  2:1214465/1‑85 (MQ=255)
                                                                                                                                            cacGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTg                                                                           <  1:1214465/85‑1 (MQ=255)
                                                                                                                                            cacGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACgg          <  1:1690549/150‑1 (MQ=255)
                                                                                                                                               gCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCTGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGCCt       >  2:851051/1‑150 (MQ=255)
                                                                                                                                                 tGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGAc                                                            <  2:342670/95‑1 (MQ=255)
                                                                                                                                                 tGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGAc                                                            >  1:342670/1‑95 (MQ=255)
                                                                                                                                                 tGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGCCTGc     <  1:1407291/150‑1 (MQ=255)
                                                                                                                                                  gATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGCCTGCa    >  2:801313/1‑150 (MQ=255)
                                                                                                                                                    tCGTGTCCTAGGACCCACTCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATCGAGACTCCTCGTAGAATTACCTGCGTGACTGCATGAACTGCAGCATGTGCGTGCAGGTCTCTGCGACCGGCATCGACATCCGTGACGGCCTGCAGa  <  1:1151109/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCCAGGCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGCCTGCAGA  >  NC_002947/4843774‑4844071

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