Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 643,952 C→G M19I (ATG→ATC acoB ← acetoin:2,6‑dichlorophenolindophenol oxidoreductase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947643,9520CG97.8% 616.8 / 2.6 184M19I (ATG→ATCacoBacetoin:2,6‑dichlorophenolindophenol oxidoreductase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base G (79/101);  total (80/104)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

Only 100 of 185 total aligned reads displayed.

GGCGCGTCAAGCACGCGGCCGGGGAATTGGTGATACAGGCCTTTGGTCACGCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGCATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATC  >  NC_002947/643809‑644098
                                                                                                                                               |                                                                                                                                                  
ggcgCGTCAAGCACGCGGCCGGGGAATTGGTGATACAGGCCTTTGGTCACGCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTcc                                                                                                                                              <  1:3286542/150‑1 (MQ=255)
  cgcgTCAAGCACGCGGCCGGGGAATTGGTGATACAGGCCTTTGGTCACGCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTg                                                                                                                                            <  2:1045089/150‑1 (MQ=255)
   gcgTCAAGCACGCGGCCGGGGAATTGGTGATACAGGCCTTTGGTCACGCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTgg                                                                                                                                           <  1:3585875/150‑1 (MQ=255)
     gTCAAGCACGCGGCCGGGGAATTGGTGATACAGGCCTTTGGTCACGCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGc                                                                                                                                         <  2:189264/150‑1 (MQ=255)
         aGCACGCGGCCGGGGAATTGGTGATACAGGCCTTTGGTCACGCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAgg                                                                                                                                     <  1:2220537/150‑1 (MQ=255)
          gCACGCGGCCGGGGAATTGGTGATACAGGCCTTTGGTCACGCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAggg                                                                                                                                    >  1:274157/1‑150 (MQ=255)
                 gCCGGGGAATTGGTGATACAGGCCTTTGGTCACGCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACTTCTTCGCCAATGATTAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGtt                                                                                                                             <  1:1144725/150‑1 (MQ=255)
                 gCCGGGGAATTGGTGATACAGGCCTTTGGTCACGCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGtt                                                                                                                             <  2:2743680/150‑1 (MQ=255)
                 gCCGGGGAATTGGTGATACAGGCCTTTGGTCACGCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGCATTTCCTGGGCCAGGGCCTCGtt                                                                                                                             <  2:1013805/150‑1 (MQ=255)
                   cGGGGAATTGGTGATACAGGCCTTTGGTCACGCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGAttt                                                                                                                                                <  1:314066/129‑1 (MQ=255)
                    ggggAATTGGTGATACAGGCCTTTGGTCACGCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGAt                                                                                                                          >  2:657729/1‑150 (MQ=255)
                     gggAATTGGTGATACAGGCCTTTGGTCACGCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATg                                                                                                                         <  1:1522639/150‑1 (MQ=255)
                     gggAATTGGTGATACAGGCCTTTGGTCACGCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGCATTTCCTGGGCCAGGGCCTCGTTGATg                                                                                                                         <  2:3524147/150‑1 (MQ=255)
                        aaTTGGTGATACAGGCCTTTGGTCACGCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGcc                                                                                                                      >  2:779406/1‑150 (MQ=255)
                            ggTGATACAGGCCTTTGGTCACGCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCctgct                                                                                                                  <  2:1176071/150‑1 (MQ=255)
                               gATACAGGCCTTTGGTCACGCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGt                                                                                                               >  1:2856345/1‑150 (MQ=255)
                                     ggCCTTTGGTCACGCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGa                                                                                                         <  2:1592161/150‑1 (MQ=255)
                                      gCCTTTGGTCACGCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGAt                                                                                                        >  1:696719/1‑150 (MQ=255)
                                        cTTTGGTCACGCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGAtct                                                                                                      >  1:2459325/1‑150 (MQ=255)
                                            tgTCACGCCCAACACGCCACCCCCGGCGTACTCGTCGCGTGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGAtctttct                                                                                                  <  2:920054/149‑1 (MQ=255)
                                             gTCACGCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGAtctttctt                                                                                                 <  1:38928/150‑1 (MQ=255)
                                              tCACGCCCAACACGCCACCCCAGGCGTCATCCTCGACAAGTGCAGGGGTAGCGGCGACGAGGGGGTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTg                                                                                                >  1:1817311/1‑150 (MQ=255)
                                               cACGCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGc                                                                                               >  2:630798/1‑150 (MQ=255)
                                                  gCCCAACACGCCCCCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGACAATGATGAATACTGTGCTGGCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGGTGGCCGGCTGGTAGGTGATCTTTTTTGCCAt                                                                                            >  1:536470/1‑150 (MQ=255)
                                                  gCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTgg                                                                                                                >  1:141577/1‑130 (MQ=255)
                                                  gCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCAt                                                                                            >  2:459935/1‑150 (MQ=255)
                                                  gCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCAt                                                                                            >  1:316128/1‑150 (MQ=255)
                                                    ccAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGt                                                                                                                              <  1:3395367/114‑1 (MQ=255)
                                                    ccAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGt                                                                                                                              >  2:3395367/1‑114 (MQ=255)
                                                      aaCACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGtt                                                                                                                             <  2:1608020/113‑1 (MQ=255)
                                                      aaCACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCAtggtg                                                                                        >  1:2146317/1‑150 (MQ=255)
                                                      aaCACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTAGtt                                                                                                                             >  1:1608020/1‑113 (MQ=255)
                                                         acGCCCCCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGtt                                                                                     <  2:2459325/150‑1 (MQ=255)
                                                         acGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGtt                                                                                     <  2:3059652/150‑1 (MQ=255)
                                                         acGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGtt                                                                                     <  1:421428/150‑1 (MQ=255)
                                                           gCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTtct                                                                                   >  2:1127966/1‑150 (MQ=255)
                                                            ccACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTtctc                                                                                  <  1:1040374/150‑1 (MQ=255)
                                                            ccACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTtctc                                                                                  <  1:630798/150‑1 (MQ=255)
                                                               ccccCGGCGCCCTCCTCGAAGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCCCCtg                                                                               <  2:2550405/150‑1 (MQ=255)
                                                                cccAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGAtctttctt                                                                                                 <  1:2060664/131‑1 (MQ=255)
                                                                cccAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGAtctttctt                                                                                                 >  2:2060664/1‑131 (MQ=255)
                                                                cccAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCtgt                                                                              >  2:1663458/1‑150 (MQ=255)
                                                                 ccAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCAtggtg                                                                                        <  2:1248389/139‑1 (MQ=255)
                                                                 ccAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCAtggtg                                                                                        >  1:1248389/1‑139 (MQ=255)
                                                                 ccAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCtgtt                                                                             <  2:316128/150‑1 (MQ=255)
                                                                 ccAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCtgtt                                                                             <  1:779406/150‑1 (MQ=255)
                                                                    ggCGTCATCCTCGCCGGGCGCACCGGCACCGGCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTgttgtt                                                                          <  1:1211866/150‑1 (MQ=255)
                                                                     gCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGAtctt                                                                                                     >  2:544132/1‑122 (MQ=255)
                                                                       gTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGCATTTCCTGGGCCAGGGCCTCGTTGATg                                                                                                                         <  2:2866718/100‑1 (MQ=255)
                                                                        tCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTg                                                                      <  1:486704/150‑1 (MQ=255)
                                                                         cATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATg                                                                                                                         >  2:1944981/1‑98 (MQ=255)
                                                                           tCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTc                                                                         <  2:1457052/144‑1 (MQ=255)
                                                                           tCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTaa                                                                   >  1:1713681/1‑150 (MQ=255)
                                                                              tCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTGATcc                                                                >  2:2366991/1‑150 (MQ=255)
                                                                              tCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATcc                                                                >  2:1347315/1‑150 (MQ=255)
                                                                                gCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCtt                                                                       <  1:514250/141‑1 (MQ=255)
                                                                                  cGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCgg                                                            <  2:1817311/150‑1 (MQ=255)
                                                                                  cGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTAGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCgg                                                            >  1:2875648/1‑150 (MQ=255)
                                                                                   gggCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGc                                                           <  1:1347315/150‑1 (MQ=255)
                                                                                    ggCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCt                                                          >  1:646078/1‑150 (MQ=255)
                                                                                     gcgcACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGt                                                                                                                              <  1:2801608/81‑1 (MQ=255)
                                                                                        caccggcaccgCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAggg                                                     <  2:490725/150‑1 (MQ=255)
                                                                                         accggcaccgCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAggg                                                     <  2:3297489/150‑1 (MQ=255)
                                                                                          ccggcaccgCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATATCATTGGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAgggg                                                    >  1:217065/1‑150 (MQ=255)
                                                                                           cggcaccgCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATg                                                                                                                         >  1:3020509/1‑80 (MQ=255)
                                                                                                      gCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGGGTGTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAg                                        >  1:1971230/1‑150 (MQ=255)
                                                                                                        gACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTcc                                                                                 <  1:1428198/107‑1 (MQ=255)
                                                                                                          cGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGt                                    >  2:3027022/1‑150 (MQ=255)
                                                                                                           gTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCtgtt                                                                             >  2:1418033/1‑108 (MQ=255)
                                                                                                            tCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTGATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCg                                  <  1:2366991/150‑1 (MQ=255)
                                                                                                              ttCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTcc                                                                                 >  1:2212392/1‑101 (MQ=255)
                                                                                                              ttCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTcc                                                                                 <  2:2212392/101‑1 (MQ=255)
                                                                                                              ttCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAgg                                                      >  2:2278472/1‑128 (MQ=255)
                                                                                                              ttCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAgg                                                      <  1:2278472/128‑1 (MQ=255)
                                                                                                               tCGCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTGCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTg                               <  2:414910/150‑1 (MQ=255)
                                                                                                                 gCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGc                                               <  1:3449275/132‑1 (MQ=255)
                                                                                                                 gCCAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGc                                               >  2:3449275/1‑132 (MQ=255)
                                                                                                                  ccAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCtt                                                                       <  2:2121392/107‑1 (MQ=255)
                                                                                                                  ccAATGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCtt                                                                       >  1:2121392/1‑107 (MQ=255)
                                                                                                                    aaTGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCtgtt                                                                             <  2:3172254/99‑1 (MQ=255)
                                                                                                                    aaTGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCtgtt                                                                             >  1:3172254/1‑99 (MQ=255)
                                                                                                                    aaTGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAg                          <  1:459935/150‑1 (MQ=255)
                                                                                                                    aaTGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAg                          >  1:2538138/1‑150 (MQ=255)
                                                                                                                    aaTGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAg                          <  1:2653822/150‑1 (MQ=255)
                                                                                                                    aaTGATGAATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAg                          <  1:2653919/150‑1 (MQ=255)
                                                                                                                       gatgaATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGAtctttc                                                                                                   >  2:2440860/1‑74 (MQ=255)
                                                                                                                       gatgaATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGAtctttc                                                                                                   <  1:2440860/74‑1 (MQ=255)
                                                                                                                          gaATGCTGTGCTGTCGCGGCGGATTTCCCGGGCCCGGGCCTGGTTTTTGTTTTGCTTGTAGCTGAGCTTTCTGGCCATGGTGGTGCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTcggc                    <  2:217065/150‑1 (MQ=255)
                                                                                                                          gaATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCcgc                                                              >  2:2348501/1‑108 (MQ=255)
                                                                                                                          gaATACTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCcgc                                                              <  1:2348501/108‑1 (MQ=255)
                                                                                                                              aCTGTGCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGc                <  2:2856345/150‑1 (MQ=255)
                                                                                                                                tgtgCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGt                                          <  1:1382961/122‑1 (MQ=255)
                                                                                                                                tgtgCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGt                                          >  2:1382961/1‑122 (MQ=255)
                                                                                                                                tgtgCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCtg              >  2:1484994/1‑150 (MQ=255)
                                                                                                                                  tgCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGc                                                                                               <  2:2816917/67‑1 (MQ=255)
                                                                                                                                   gCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATc                                                                 >  1:580353/1‑96 (MQ=255)
                                                                                                                                   gCTGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATc                                                                 <  2:580353/96‑1 (MQ=255)
                                                                                                                                     tGTCGCGGCGGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCt         >  1:1064518/1‑150 (MQ=255)
                                                                                                                                       ccgcgGCGGATTTCCGGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTGCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTg       <  2:536470/149‑1 (MQ=255)
                                                                                                                                         gcggcgGATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCGGTGGCTTggg     >  2:1135968/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GGCGCGTCAAGCACGCGGCCGGGGAATTGGTGATACAGGCCTTTGGTCACGCCCAACACGCCACCCCAGGCGTCATCCTCGCCGGGCGCACCGGCACCGCCGGCGACGTCTTCGCCAATGATGAATACTGTGCTGTCGCGGCGCATTTCCTGGGCCAGGGCCTCGTTGATGGCCTGCTGGTAGCTGATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATC  >  NC_002947/643809‑644098

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