Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 635,845 A→C 7.0% E22D (GAA→GAC) ‡ mpl → UDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947635,8450AC 7.0% 564.2 / 24.1 182E22D (GAA→GAC) ‡mplUDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base A (75/94);  new base C (7/6);  total (82/100)
Fisher's exact test for biased strand distribution p-value = 5.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

Only 100 of 186 total aligned reads displayed.

ATTACCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTC  >  NC_002947/635696‑635994
                                                                                                                                                     |                                                                                                                                                     
aTTACCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGaa                                                                                                                                                       <  2:4231675/150‑1 (MQ=255)
aTTACCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGaa                                                                                                                                                       >  2:364853/1‑150 (MQ=255)
  tACCAGGGGAGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACt                                                                                                                                                     <  1:3222049/150‑1 (MQ=255)
   aCCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACtt                                                                                                                                                    <  2:2361802/150‑1 (MQ=255)
    ccAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTg                                                                                                                                                   <  2:4416566/150‑1 (MQ=255)
        gggCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGTTCCTTGGCCa                                                                                                                                               <  1:2950344/150‑1 (MQ=255)
        gggCGACTGAGCGAAACAGCGGGTGCGATCTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGTTCCTTGGCCa                                                                                                                                               <  2:4519538/150‑1 (MQ=255)
           cGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGTTCCTTGGCCagca                                                                                                                                            <  2:5067174/150‑4 (MQ=255)
           cGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACcg                                                                                                                                            <  1:626186/150‑1 (MQ=255)
            gACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGAGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACcgc                                                                                                                                           >  2:5063137/1‑150 (MQ=255)
               tGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTc                                                                                                                                        >  2:2468098/1‑150 (MQ=255)
               tGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATACACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTc                                                                                                                                        >  2:865669/1‑150 (MQ=255)
                     aaaCAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGc                                                                                                                                  <  2:3055270/150‑1 (MQ=255)
                        cAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGctct                                                                                                                               <  1:3946048/150‑1 (MQ=255)
                         aGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTg                                                                                                                              >  1:4917364/1‑150 (MQ=255)
                           cGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGAc                                                                                                                            <  2:3596885/150‑1 (MQ=255)
                            gggTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACg                                                                                                                           <  2:71179/150‑1 (MQ=255)
                             ggTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGc                                                                                                                          >  2:4200940/1‑150 (MQ=255)
                               tGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATAGTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCa                                                                                                                        <  2:4643174/150‑1 (MQ=255)
                                  gATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACg                                                                                                                     >  2:1536882/1‑150 (MQ=255)
                                        tCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATc                                                                                                               <  1:100194/150‑1 (MQ=255)
                                          gcgcCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATccc                                                                                                             <  1:2251560/150‑1 (MQ=255)
                                          gcgcCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATccc                                                                                                             <  1:4469038/150‑1 (MQ=255)
                                          gcgcCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACAGCGTCACCGGCTCTGACGCCAACGTCTATccc                                                                                                             >  1:1582895/1‑150 (MQ=255)
                                              ccTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGa                                                                                                         >  1:4353023/1‑150 (MQ=255)
                                               cTGTGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTTACGCCAACGTCTATCCCCCGAt                                                                                                        <  2:276202/150‑1 (MQ=255)
                                               cTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGAt                                                                                                        <  1:2051234/150‑1 (MQ=255)
                                               cTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGAt                                                                                                        <  2:4687632/150‑1 (MQ=255)
                                                 gggCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATAGTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGa                                                                                                      <  2:1733892/150‑1 (MQ=255)
                                                   gCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGc                                                                                                    <  2:2537495/150‑1 (MQ=255)
                                                     tACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCAc                                                                                                  <  1:2107325/150‑1 (MQ=255)
                                                      aCCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCAcc                                                                                                 <  1:760672/150‑1 (MQ=255)
                                                         cTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAg                                                                                              <  2:4353023/150‑1 (MQ=255)
                                                           gCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCt                                                                                            >  2:4530083/1‑150 (MQ=255)
                                                             cTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCg                                                                                          >  2:3185441/1‑150 (MQ=255)
                                                                gATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAg                                                                                       <  1:4679794/150‑1 (MQ=255)
                                                                gATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAg                                                                                       >  1:3565703/1‑150 (MQ=255)
                                                                gATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAg                                                                                       <  1:3749167/150‑1 (MQ=255)
                                                                      ttGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGTTCCTTGGCCa                                                                                                                                               >  1:2434706/1‑88 (MQ=255)
                                                                      ttGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAgg                                                                                 <  2:268024/150‑1 (MQ=255)
                                                                      ttGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAgg                                                                                 <  1:1316912/150‑1 (MQ=255)
                                                                      ttGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAATTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAgg                                                                                 <  1:5146962/150‑1 (MQ=255)
                                                                      ttGCGAGGTCGACCATGCATATTCACAGTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAgg                                                                                 <  2:1948418/150‑1 (MQ=255)
                                                                         cGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTa                                                                              <  2:4058089/150‑1 (MQ=255)
                                                                          gAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAATGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTAt                                                                             >  2:573492/1‑150 (MQ=255)
                                                                            ggTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATcccc                                                                                                            >  1:1172865/1‑117 (MQ=255)
                                                                              tCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGaa                                                                         >  1:2337742/1‑150 (MQ=255)
                                                                                gACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACt                                                                       >  2:2439493/1‑150 (MQ=255)
                                                                                    aTGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATc                                                                            <  2:2798575/141‑1 (MQ=255)
                                                                                    aTGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATc                                                                            >  1:2798575/1‑141 (MQ=255)
                                                                                    aTGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGAcc                                                                   >  2:5118344/1‑150 (MQ=255)
                                                                                    aTGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGAcc                                                                   >  2:5261779/1‑150 (MQ=255)
                                                                                     tGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGAccc                                                                  <  2:437462/150‑1 (MQ=255)
                                                                                      gCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCa                                                                 >  2:3808777/1‑150 (MQ=255)
                                                                                         tatTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAgg                                                              <  1:917159/150‑1 (MQ=255)
                                                                                            tCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAATGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTa                                                           <  1:573492/150‑1 (MQ=255)
                                                                                            tCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGaa                                                                         >  2:1987041/1‑136 (MQ=255)
                                                                                            tCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTa                                                           >  2:3631797/1‑150 (MQ=255)
                                                                                            tCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTa                                                           >  1:1792954/1‑150 (MQ=255)
                                                                                             cacaTTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTAt                                                          <  2:4936929/150‑1 (MQ=255)
                                                                                              acaTTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATg                                                         >  1:2393395/1‑150 (MQ=255)
                                                                                                    tcGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCgg                                                   >  1:2275934/1‑150 (MQ=255)
                                                                                                    tcGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCgg                                                   >  1:1265281/1‑150 (MQ=255)
                                                                                                       gTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGccc                                                >  2:1184535/1‑150 (MQ=255)
                                                                                                          tttGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGc                                             <  1:76145/150‑1 (MQ=255)
                                                                                                           ttGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCt                                            <  1:3092145/150‑1 (MQ=255)
                                                                                                           ttGCGGCACTTTCATGGGGTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCt                                            <  2:5194873/150‑1 (MQ=255)
                                                                                                            tGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTg                                           <  1:1496720/150‑1 (MQ=255)
                                                                                                            tGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTg                                           <  1:1842796/150‑1 (MQ=255)
                                                                                                            tGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTg                                           <  1:1842839/150‑1 (MQ=255)
                                                                                                             gCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTgg                                          <  1:3951719/150‑1 (MQ=255)
                                                                                                              cGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGa                                         <  1:4197916/150‑1 (MQ=255)
                                                                                                                gCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATc                                       <  1:998833/150‑1 (MQ=255)
                                                                                                                  aCTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCg                                     <  2:2813031/150‑1 (MQ=255)
                                                                                                                    tttCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTg                                           <  2:1433423/142‑1 (MQ=255)
                                                                                                                    tttCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGgc                                   <  2:3552577/150‑1 (MQ=255)
                                                                                                                    tttCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACGTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGAGGCGCACCCAGCTCGATGCCCAGGGTAGCGAACTGACCCAAGGCTCTGAGCCGGCCCAGCTGGCTCCggg                                   >  2:5393087/1‑149 (MQ=255)
                                                                                                                      tCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGgcgc                                 <  1:222712/150‑1 (MQ=255)
                                                                                                                         tGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAg                              >  1:3971683/1‑150 (MQ=255)
                                                                                                                             ttCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGt                                                                                                                    <  1:3604009/60‑1 (MQ=255)
                                                                                                                              tCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCtg                         <  2:4124245/150‑1 (MQ=255)
                                                                                                                                gCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCtggt                       <  1:2853366/150‑1 (MQ=255)
                                                                                                                                gCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCtggg                       >  2:3305301/1‑149 (MQ=255)
                                                                                                                                 cTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCtggtg                      <  2:1696332/150‑1 (MQ=255)
                                                                                                                                  tGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAGGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCtggtgg                     <  1:1051566/150‑1 (MQ=255)
                                                                                                                                  tGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCtggtgg                     <  2:2053264/150‑1 (MQ=255)
                                                                                                                                   ggCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTggtggt                    >  2:4787012/1‑150 (MQ=255)
                                                                                                                                   ggCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTggtggt                    >  1:4858378/1‑150 (MQ=255)
                                                                                                                                    gCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTc                   >  2:2854752/1‑150 (MQ=255)
                                                                                                                                       gcggtGGCCAAGTTCCTTGGCCAGCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCtggtgg                     <  2:2138029/141‑1 (MQ=255)
                                                                                                                                       gcggtGGCCAAGTTCCTTGGCCAGCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATc                >  1:4723180/5‑150 (MQ=255)
                                                                                                                                       gcggtGGCCAAGTTCCATGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGGATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATc                >  2:5292963/5‑150 (MQ=255)
                                                                                                                                         ggtGGCCAAGTTCCTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCgg              >  1:5308011/3‑150 (MQ=255)
                                                                                                                                         gCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGAACAGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCgg              >  2:4897182/1‑150 (MQ=255)
                                                                                                                                          gtGGCCAAGTTCCTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGc             >  2:3393576/2‑150 (MQ=255)
                                                                                                                                          gtGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGc             >  1:4858474/2‑150 (MQ=255)
                                                                                                                                           tGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCa            >  2:1100809/1‑150 (MQ=255)
                                                                                                                                             gCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAAc          >  2:3494483/1‑150 (MQ=255)
                                                                                                                                                 aGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATACCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCa      <  2:2693584/150‑1 (MQ=255)
                                                                                                                                                     aCTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTc  >  1:1404801/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTACCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGTC  >  NC_002947/635696‑635994

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: