Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,592,029 T→G 15.2% intergenic (‑94/‑109) PP_1395 ← / → PP_1396 AraC family transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,592,0290TG15.2% 410.5 / 59.4 163intergenic (‑94/‑109)PP_1395/PP_1396AraC family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (68/70);  new base G (13/12);  total (81/82)
Fisher's exact test for biased strand distribution p-value = 8.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

Only 100 of 166 total aligned reads displayed.

CAGCTCCAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTCACCCACCTTGCTTTCATCCG  >  NC_002947/1591882‑1592178
                                                                                                                                                   |                                                                                                                                                     
cAGCTCCAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGc                                                                                                                                                     <  2:1990316/150‑1 (MQ=255)
 aGCTCCAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCa                                                                                                                                                    >  2:2069001/1‑150 (MQ=255)
  gCTCCAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAg                                                                                                                                                   <  2:3864107/150‑1 (MQ=255)
   cTCCAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAgg                                                                                                                                                  >  1:2735716/1‑150 (MQ=255)
    tCCAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGa                                                                                                                                                 <  1:4277428/150‑1 (MQ=255)
    tCCAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCAGCGCGGCAGGa                                                                                                                                                 <  2:3598679/150‑1 (MQ=255)
      cAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATc                                                                                                                                               <  2:2387801/150‑1 (MQ=255)
      cAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATc                                                                                                                                               >  2:4744717/1‑150 (MQ=255)
       aGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCa                                                                                                                                              <  1:917772/150‑1 (MQ=255)
       aGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCAGCGCGGCAGGATCa                                                                                                                                              <  2:525050/150‑1 (MQ=255)
        ggCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCAGCGCGGCAGGATCAg                                                                                                                                             <  1:1061487/150‑1 (MQ=255)
         gcgcTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAgg                                                                                                                                            <  1:105229/150‑1 (MQ=255)
         gcgcTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAgg                                                                                                                                            <  1:4744717/150‑1 (MQ=255)
         gcgcTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAgg                                                                                                                                            <  1:2842133/150‑1 (MQ=255)
         gcgcTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAACTGCATGTCGATTGCCTGATCCTGCAGCGCGGCAGGATCAgg                                                                                                                                            <  2:3776652/150‑1 (MQ=255)
          cgcTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGc                                                                                                                                           <  2:696938/150‑1 (MQ=255)
              gCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATc                                                                                                                                       <  1:4261152/150‑1 (MQ=255)
               cAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCAGCGCGGCAGGATCAggagtgcg                                                                                                                                      >  2:1522731/1‑144 (MQ=255)
                  ttCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTCGATTGCCTGATCCTGCAGCGCGGCAGGATCAGGCAATCcaca                                                                                                                                   <  1:3177869/150‑1 (MQ=255)
                    ccAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATg                                                                                                                                 <  2:2858740/150‑1 (MQ=255)
                           cggttcggtGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAAt                                                                                                                          <  1:3450279/150‑1 (MQ=255)
                           cggttcggtGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTCGATTGCCTGATCCTGCAGCGCGGCAGGATCAGGCAATCCACATGCAgcttt                                                                                                                          <  1:4442689/150‑5 (MQ=37)
                            ggttcggtGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCt                                                                                                                                                       >  2:708198/1‑120 (MQ=255)
                              ttcggtGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGc                                                                                                                       >  2:2516236/1‑150 (MQ=255)
                              ttcggtGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAACTGCATGTCGATTGCCTGATCCTGCAGCGCGGCAGGATCAGGCAATCCACATGCAgctttcga                                                                                                                       <  1:4893622/150‑8 (MQ=25)
                                     gACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCAGCGCGGCAGGATCa                                                                                                                                              >  2:4358654/1‑120 (MQ=255)
                                     gACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCAGCGCGGCAGGATCa                                                                                                                                              <  1:4358654/120‑1 (MQ=255)
                                       cGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAAcgcgc                                                                                                              >  1:3017413/1‑150 (MQ=255)
                                       cGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAAcgcgc                                                                                                              >  1:3017430/1‑150 (MQ=255)
                                           gCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAAcgcgccgcg                                                                                                          >  1:1332240/1‑150 (MQ=255)
                                             ggtggtCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCcgcgcg                                                                                                        >  2:4953789/1‑150 (MQ=255)
                                             ggtggtCATGATGGCTTTCATCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCTGGCAATCGACATGCAGTAATGGCCTAACGCGCcgcgcg                                                                                                        >  2:3376190/1‑150 (MQ=255)
                                              gtggtCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCgg                                                                                                       <  1:1780028/150‑1 (MQ=255)
                                                ggtCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGAAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGt                                                                                                     <  1:3185918/150‑1 (MQ=255)
                                                ggtCATGAGGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGCTGGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGt                                                                                                     <  1:2713199/150‑1 (MQ=255)
                                                  tCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGa                                                                                                   <  2:3033243/150‑1 (MQ=255)
                                                  tCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCACGATCAGGCAATCGACATGCAGTAATGGCCAAACGCGCCGCGCGGGTGa                                                                                                   >  1:209202/1‑150 (MQ=255)
                                                    atgatgGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAgc                                                                                                 <  1:2097424/150‑1 (MQ=255)
                                                    atgatgGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAgc                                                                                                 <  1:566835/150‑1 (MQ=255)
                                                    atgatgGCTTTCCACTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAgc                                                                                                 <  1:566620/150‑1 (MQ=255)
                                                            tttCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTCAcc                                                                                         >  1:3810329/1‑150 (MQ=255)
                                                            tttCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCAGCGCGGCAGGATCa                                                                                                                                              >  2:3623030/1‑97 (MQ=255)
                                                            tttCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCAGCGCGGCAGGATCa                                                                                                                                              <  1:3623030/97‑1 (MQ=255)
                                                              tCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCtt                                                                                       <  2:5156961/150‑1 (MQ=255)
                                                              tCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCtt                                                                                       <  2:2428686/150‑1 (MQ=255)
                                                                ctctGCTGACCGACGCAGCCTAGGGTGGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGCTCCTGCCGCGCTGCAGGAGCAGGCAAGCGCCCTGCCGTAAGGGCCTCACGCGCCGCGCGGGTGCGCGCCTAAccgggg                                                                                     >  2:4554088/1‑146 (MQ=255)
                                                                 tctGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAgg                                                                                                                                            >  1:1354870/1‑94 (MQ=255)
                                                                 tctGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGAACAGGCAATCGACATGAAGAAATGGCCTAAAGCGCCGCGCGGGTGAGCGCCTAACCTTGGc                                                                                    >  2:2608507/1‑150 (MQ=255)
                                                                    gctgACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGc                                                                                 >  1:176911/1‑150 (MQ=255)
                                                                     ctgACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACAAGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGcc                                                                                >  1:2519235/1‑150 (MQ=255)
                                                                      tgACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCt                                                                               <  2:5197352/150‑1 (MQ=255)
                                                                      tgACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTCGATTGCCTGATCCTGCAGCGCGGCAGGATCAGGCAATCCACATGCAgcttt                                                                                                                          >  2:3801110/1‑103 (MQ=18)
                                                                      tgACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTCGATTGCCTGATCCTGCAGCGCGGCAGGATCAGGCAATCCACATGCAgcttt                                                                                                                          <  1:3801110/107‑5 (MQ=18)
                                                                        aCCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTgg                                                                             >  2:1657813/1‑150 (MQ=255)
                                                                         ccGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGc                                                                            <  2:4783940/150‑1 (MQ=255)
                                                                          cGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGCTCAGGCCATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCGAACCTGGGCAGCCTGGCt                                                                           >  2:5318553/1‑150 (MQ=255)
                                                                          cGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCCGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCgg                                                                                                       <  1:3809264/122‑1 (MQ=255)
                                                                              gcagcaTAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTtc                                                                       >  2:2784600/1‑150 (MQ=255)
                                                                              gcagcaTAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTtc                                                                       >  2:2784606/1‑150 (MQ=255)
                                                                              gcagcaTAGGTGTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTtc                                                                       >  2:4559483/1‑150 (MQ=255)
                                                                               cagcaTAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTtcc                                                                      <  2:2453080/150‑1 (MQ=255)
                                                                                 gcaTAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTtccgt                                                                    >  1:1187464/1‑150 (MQ=255)
                                                                                  caTAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCGGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGcc                                                                                <  1:712684/137‑1 (MQ=255)
                                                                                  caTAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCCATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTtccgtc                                                                   >  1:1892966/1‑150 (MQ=255)
                                                                                   aTAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGt                                                                                                     >  2:4189312/1‑115 (MQ=255)
                                                                                   aTAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTtccgtcc                                                                  >  2:1015830/1‑150 (MQ=255)
                                                                                     aGGTTTGGGCAAAGCCGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGc                                                                <  1:2801132/150‑1 (MQ=255)
                                                                                         ttGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGc                                                            <  1:4418940/150‑1 (MQ=255)
                                                                                           gggCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCGAACCGTGGCAGCCTGGCTTTTCCGTCCGCCGGCtt                                                          >  1:2277486/1‑150 (MQ=255)
                                                                                             gCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCa                                                        >  2:1803329/1‑150 (MQ=255)
                                                                                             gCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGACTTTTCCGTCCGCCTGCTTCa                                                        >  2:1985389/1‑150 (MQ=255)
                                                                                              cAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCGAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCaa                                                       >  1:643179/1‑150 (MQ=255)
                                                                                              cAAAGACGGGCGCTAGTCGAAAGCTGCATGTCGATTGCCTGATCCTGCCGCGCTGCAGGATCAg                                                                                                                                             <  1:752285/64‑1 (MQ=255)
                                                                                               aaaGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAg                                                      <  2:176911/150‑1 (MQ=255)
                                                                                               aaaGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAg                                                      <  2:1635562/150‑1 (MQ=255)
                                                                                                   aCGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGgaga                                                  <  1:2282914/150‑1 (MQ=255)
                                                                                                     gggCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGAtt                                                <  2:1390841/150‑1 (MQ=255)
                                                                                                        cgcTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATtgtg                                             <  2:843963/150‑1 (MQ=255)
                                                                                                        cgcTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATtgtg                                             >  1:4792413/1‑150 (MQ=255)
                                                                                                           tAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCaa                                          >  2:890052/1‑150 (MQ=255)
                                                                                                                 aaaGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAg                                                      <  1:5162963/132‑1 (MQ=255)
                                                                                                                 aaaGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGAccc                                    <  1:2329246/150‑1 (MQ=255)
                                                                                                                     cTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTaa                                <  2:3360274/150‑1 (MQ=255)
                                                                                                                      tGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTaaa                               >  1:3622905/1‑150 (MQ=255)
                                                                                                                          tgtgGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAaacaac                           <  2:4905103/150‑1 (MQ=255)
                                                                                                                           gtgGATTGCCTGCTCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTACCGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAAcc                          <  2:3177624/150‑1 (MQ=255)
                                                                                                                           gtcgATTGCCTGATCCTGCAGCGCGGCAGGATCAGGCAATCGACATGCAGTAATGGCCGAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGAATGTGCAATGACCCTTAAACAAcc                          >  2:752174/4‑150 (MQ=38)
                                                                                                                            tgGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCg                         <  2:1400606/150‑1 (MQ=255)
                                                                                                                            tgGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCg                         >  1:4044799/1‑150 (MQ=255)
                                                                                                                              gATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCCGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGt                       >  2:5067611/1‑150 (MQ=255)
                                                                                                                              gATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGt                       >  1:2954545/1‑150 (MQ=255)
                                                                                                                                   ccTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcaccc                  >  2:5044845/1‑150 (MQ=255)
                                                                                                                                    cTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcaccca                 >  2:928599/1‑150 (MQ=255)
                                                                                                                                    cTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcaccca                 >  2:1889551/1‑150 (MQ=255)
                                                                                                                                    cTGATCCTGCAGCGCGGCAGGATCAGGCAATCCACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcaccca                 >  2:1269208/1‑150 (MQ=255)
                                                                                                                                     tGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCTATGACCCTTAAACAACCGGTcacccac                >  2:3239283/1‑150 (MQ=255)
                                                                                                                                       aTCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTCACCCACCt              >  2:1160393/1‑150 (MQ=255)
                                                                                                                                          cTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTCACCCACCTTGc           <  2:2331524/150‑1 (MQ=255)
                                                                                                                                            gCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTCACCCACCTTGCtt         >  1:3071706/1‑150 (MQ=255)
                                                                                                                                                   tGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGTAATGACCCTTAAACAACCGGTCACCCACCTTGCTTTCATCCg  >  1:2470555/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CAGCTCCAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTCACCCACCTTGCTTTCATCCG  >  NC_002947/1591882‑1592178

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