Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,024,075 A→C 5.1% C44G (TGC→GGC)  dadX ← alanine racemase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,024,0750AC 5.1% 1306.9 / 40.0 379C44G (TGC→GGC) dadXalanine racemase
Reads supporting (aligned to +/- strand):  ref base A (192/167);  new base C (12/8);  total (204/175)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 411 total aligned reads displayed.

GTGCTCGGCGATCAGCGCCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGT  >  NC_002947/6023926‑6024224
                                                                                                                                                     |                                                                                                                                                     
gTGCTCGGCGATCAGCGCCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGgcgca                                                                                                                                                       <  2:1375080/150‑1 (MQ=255)
  gCTCGGCGATCAGCGCCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGgcgcagc                                                                                                                                                     >  2:3635333/1‑150 (MQ=255)
    tCGGCGATCAGCGCCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGgcgcagcgc                                                                                                                                                   >  2:4840429/1‑150 (MQ=255)
      ggCGATCAGCGCCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCAc                                                                                                                                                 >  1:1758753/1‑150 (MQ=255)
         gATCAGCGCCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGgc                                                                                                                                              <  2:2065853/150‑1 (MQ=255)
              gcgcCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGt                                                                                                                                         >  1:1359086/1‑150 (MQ=255)
              gcgcCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGt                                                                                                                                         >  1:6755969/1‑150 (MQ=255)
              gcgcCAGGTCGCTGGCTTCAAAAAAACCTGCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCACGGCGCCGTGCGCTGCGccct                                                                                                                                         <  1:5194066/150‑3 (MQ=255)
               cgcCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTg                                                                                                                                        <  2:1285430/150‑1 (MQ=255)
               cgcCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGGACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCACGGCGCCGTGCGCTGCGCCCTg                                                                                                                                        <  1:3893793/150‑1 (MQ=255)
                gcCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTgg                                                                                                                                       >  1:2813554/1‑150 (MQ=255)
                  cAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGcc                                                                                                                                     >  1:184255/1‑150 (MQ=255)
                    gTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGt                                                                                                                                   <  2:3370261/150‑1 (MQ=255)
                     ttCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGACGGCCACGGCGCCGTGCGCTGCGCCCTGGccctc                                                                                                                                  <  1:2157715/150‑4 (MQ=255)
                        gCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGc                                                                                                                               >  1:708831/1‑150 (MQ=255)
                          tGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCAt                                                                                                                             >  2:1100051/1‑150 (MQ=255)
                            gCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCg                                                                                                                           <  1:6428625/150‑1 (MQ=255)
                              ttCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGc                                                                                                                         >  1:5016880/1‑150 (MQ=255)
                              ttCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGc                                                                                                                         >  2:2778564/1‑150 (MQ=255)
                              ttCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGc                                                                                                                         <  1:6472393/150‑1 (MQ=255)
                              ttCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGc                                                                                                                         >  2:3911260/1‑150 (MQ=255)
                               tCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGcc                                                                                                                        >  2:1583560/1‑150 (MQ=255)
                               tCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGcc                                                                                                                        >  1:2335520/1‑150 (MQ=255)
                                cAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCt                                                                                                                       <  2:30746/150‑1 (MQ=255)
                                cAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCt                                                                                                                       <  2:5211529/150‑1 (MQ=255)
                                 aaaaaaaCCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCtt                                                                                                                      >  2:5483668/1‑150 (MQ=255)
                                 aaaaaaaCCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCtt                                                                                                                      >  1:4045101/1‑150 (MQ=255)
                                 aaaaaaaCCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCtt                                                                                                                      >  1:1574551/1‑150 (MQ=255)
                                 aaaaaaaCCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCtt                                                                                                                      <  1:6310349/150‑1 (MQ=255)
                                  aaaaaaCCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTg                                                                                                                     <  2:5912571/150‑1 (MQ=255)
                                   aaaaaCCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGa                                                                                                                    >  2:6319351/1‑150 (MQ=255)
                                   aaaaaCCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGa                                                                                                                    >  2:4625097/1‑150 (MQ=255)
                                   aaaaaCCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGa                                                                                                                    >  2:4342223/1‑150 (MQ=255)
                                        ccTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACg                                                                                                               >  2:1930482/1‑150 (MQ=255)
                                         cTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACgg                                                                                                              >  2:1007156/1‑150 (MQ=255)
                                          ttCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGc                                                                                                             >  1:2414662/1‑150 (MQ=255)
                                          ttCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGc                                                                                                             >  1:3556435/1‑150 (MQ=255)
                                          ttCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGc                                                                                                             >  1:6688613/1‑150 (MQ=255)
                                          ttCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGc                                                                                                             >  2:3237396/1‑150 (MQ=255)
                                           tCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCg                                                                                                            >  1:4787719/1‑150 (MQ=255)
                                              agcagcagTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAgg                                                                                                         <  1:4929441/150‑1 (MQ=255)
                                              agcagcagTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAgg                                                                                                         <  1:4471250/150‑1 (MQ=255)
                                              agcagcagTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGACGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAgg                                                                                                         >  2:1895233/1‑150 (MQ=255)
                                                   cagTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCa                                                                                                                 <  2:4173647/137‑1 (MQ=255)
                                                   cagTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCttt                                                                                                    <  1:4614674/150‑1 (MQ=255)
                                                   cagGACCGGGGCTTTGAGACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCttt                                                                                                    <  2:211423/150‑1 (MQ=255)
                                                    agTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTg                                                                                                   <  2:6640311/150‑1 (MQ=255)
                                                      tACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGgcgcag                                                                                                                                                      >  1:4851635/1‑97 (MQ=255)
                                                         cGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCAcc                                                                                              >  1:600970/1‑150 (MQ=255)
                                                          ggggCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCAccc                                                                                             <  2:4567494/150‑1 (MQ=255)
                                                          ggggCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCAccc                                                                                             <  1:5955660/150‑1 (MQ=255)
                                                              cTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCa                                                                                         <  1:1000338/150‑1 (MQ=255)
                                                              cTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCa                                                                                         <  2:5117395/150‑1 (MQ=255)
                                                              cTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCa                                                                                         <  2:1825878/150‑1 (MQ=255)
                                                              cTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCa                                                                                         <  1:1195616/150‑1 (MQ=255)
                                                               tttGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAg                                                                                        >  1:308588/1‑150 (MQ=255)
                                                               tttGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAg                                                                                        >  1:3199563/1‑150 (MQ=255)
                                                               tttGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCGCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTCGGCATCGGCCGTGATCACGGCGAGGGCTTTGGCACCCGTccg                                                                                        >  2:6541014/1‑148 (MQ=255)
                                                                    tACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCAc                                                                                   >  1:862585/1‑150 (MQ=255)
                                                                    tACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCAc                                                                                   >  2:1202621/1‑150 (MQ=255)
                                                                    tACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCAc                                                                                   >  2:3690152/1‑150 (MQ=255)
                                                                    tACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTGTGGCACCCGTCAGTTCAc                                                                                   >  1:3432657/1‑150 (MQ=255)
                                                                     aCCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACg                                                                                  <  2:574405/150‑1 (MQ=255)
                                                                     aCCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACg                                                                                  <  1:3675449/150‑1 (MQ=255)
                                                                      ccGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACgg                                                                                 <  2:6030641/150‑1 (MQ=255)
                                                                       cggcggCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACggg                                                                                <  2:2594845/150‑1 (MQ=255)
                                                                       cggcggCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACggg                                                                                <  2:4234689/150‑1 (MQ=255)
                                                                               cgcAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACgg                                                                        <  2:5299292/150‑1 (MQ=255)
                                                                               cgcAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACgg                                                                        <  2:4956557/150‑1 (MQ=255)
                                                                                gcAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGt                                                                       >  1:2243328/1‑150 (MQ=255)
                                                                                 cAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTa                                                                      <  1:3335081/150‑1 (MQ=255)
                                                                                    ttCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGtt                                                                   >  1:946175/1‑150 (MQ=255)
                                                                                       cAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTtgtg                                                                >  2:2406954/1‑150 (MQ=255)
                                                                                       cAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGGAGTtgtg                                                                >  1:1496398/1‑150 (MQ=255)
                                                                                         gcgcCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGgc                                                              >  1:2498761/1‑150 (MQ=255)
                                                                                           gcCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGgcgc                                                            >  2:6358294/1‑150 (MQ=255)
                                                                                              tctTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCCCGGCGCCGTGGCCGTCGGCAGCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGGAGTGGTGGCGCAgg                                                         >  2:2864203/1‑150 (MQ=255)
                                                                                                ttCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTCGTTGCGGCGCAGGGc                                                       >  2:672254/1‑150 (MQ=255)
                                                                                                ttCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGc                                                       <  1:5843946/150‑1 (MQ=255)
                                                                                                  cGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGGAGTTGTGGCGCAGGGCCt                                                     >  1:2039350/1‑150 (MQ=255)
                                                                                                     tGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGa                                                  >  1:5516706/1‑150 (MQ=255)
                                                                                                     tGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGa                                                  <  1:4184727/150‑1 (MQ=255)
                                                                                                        aGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGt                                               <  1:6599137/150‑1 (MQ=255)
                                                                                                             aCGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGGAGTTGTGGCGCAGGGCCTGGAGGTCGATc                                          >  1:5737785/1‑150 (MQ=255)
                                                                                                                      ccGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCCGGAGGt                                               <  2:619086/136‑1 (MQ=255)
                                                                                                                            gCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGa                           <  2:3095740/150‑1 (MQ=255)
                                                                                                                              ttCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACg                         >  1:1866830/1‑150 (MQ=255)
                                                                                                                               tCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGc                        <  2:1509874/150‑1 (MQ=255)
                                                                                                                                  gCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATa                     >  1:2616421/1‑150 (MQ=255)
                                                                                                                                  gCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATa                     <  2:6631971/150‑1 (MQ=255)
                                                                                                                                   ccTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAg                    >  1:1330499/1‑150 (MQ=255)
                                                                                                                                      ctacggCCACGGCGCCGTGCGCTGCGCCCTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCga                 >  1:6374147/5‑149 (MQ=255)
                                                                                                                                       gAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGc                >  1:3004387/1‑150 (MQ=255)
                                                                                                                                       gAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGc                <  1:6079170/150‑1 (MQ=255)
                                                                                                                                        acggCCACGGCGCCGTGCGCTGCGCCCTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCt               >  2:6483728/3‑150 (MQ=255)
                                                                                                                                        aGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCt               <  1:1047027/150‑1 (MQ=255)
                                                                                                                                            ccAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCt           >  2:3966633/1‑150 (MQ=255)
                                                                                                                                             cacggCGCCGTGCGCTGCGCCCTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCtt          >  2:2264973/4‑150 (MQ=255)
                                                                                                                                              aGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTa         <  1:5252053/150‑1 (MQ=255)
                                                                                                                                               gggCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTAt        <  2:708831/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGCTCGGCGATCAGCGCCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGT  >  NC_002947/6023926‑6024224

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