Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,416,998 T→C 6.9% intergenic (+141/+27) PP_5644 → / ← PP_3899 membrane protein/hicB family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,416,9980TC 6.9% 491.7 / 16.8 160intergenic (+141/+27)PP_5644/PP_3899membrane protein/hicB family protein
Reads supporting (aligned to +/- strand):  ref base T (86/63);  new base C (5/6);  total (91/69)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

Only 100 of 161 total aligned reads displayed.

GAGTAGATGGCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTCACATTGAA  >  NC_002947/4416852‑4417147
                                                                                                                                                  |                                                                                                                                                     
gAGTAGATGGCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGc                                                                                                                                                    >  2:3697598/1‑150 (MQ=255)
gAGTAGATGGCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGc                                                                                                                                                    >  2:3697554/1‑150 (MQ=255)
 aGTAGCTGGCTCTGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCGTCCCTACGTTAGCCCGCATTGCa                                                                                                                                                   <  2:4004213/150‑1 (MQ=255)
 aGTAGATGGCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCa                                                                                                                                                   <  2:4875740/150‑1 (MQ=255)
       tGGCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCACTGCAATGCgg                                                                                                                                             <  2:1990976/150‑1 (MQ=255)
        ggCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCggg                                                                                                                                            <  1:1113446/150‑1 (MQ=255)
         gCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCACTGCAATGCGGGc                                                                                                                                           <  1:4867975/150‑1 (MQ=255)
           tATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCtt                                                                                                                                         >  1:4236736/1‑150 (MQ=255)
           tATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCtt                                                                                                                                         >  1:5042888/1‑150 (MQ=255)
              gAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTc                                                                                                                                      >  2:5019285/1‑150 (MQ=255)
              gAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCACTGCAATGCGGGCtaacg                                                                                                                                      <  2:2293263/150‑5 (MQ=255)
               aaGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCCTTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCACTGCAATGCGGGCtaacgt                                                                                                                                     <  2:3288619/150‑6 (MQ=255)
                 gCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCttt                                                                                                                                   <  1:27795/150‑1 (MQ=255)
                  cACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTg                                                                                                                                  <  1:3464373/150‑1 (MQ=255)
                   aCTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTgg                                                                                                                                 >  2:3625156/1‑150 (MQ=255)
                       ctctTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGggcgg                                                                                                                             <  1:4278578/150‑1 (MQ=255)
                            gAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTgg                                                                                                                                 >  1:5444014/1‑141 (MQ=255)
                            gAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTgg                                                                                                                                 <  2:5444014/141‑1 (MQ=255)
                             aaGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTAcgc                                                                                                                       <  1:4434253/150‑1 (MQ=255)
                                 cgccATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGtga                                                                                                                   <  1:765212/150‑1 (MQ=255)
                                   ccATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTGTTCTTTGGGCGGCTACGCGtgatg                                                                                                                 >  2:3833986/1‑150 (MQ=255)
                                      tCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGAGATGAcc                                                                                                              >  1:336283/1‑150 (MQ=255)
                                        aaGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCa                                                                                                            >  2:3003839/1‑150 (MQ=255)
                                        aaGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCa                                                                                                            <  1:827462/150‑1 (MQ=255)
                                          gTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAgg                                                                                                          >  2:2544386/1‑150 (MQ=255)
                                             aTCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCCTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTAc                                                                                                       >  1:879812/1‑150 (MQ=255)
                                              tCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAgg                                                                                                          >  2:120197/1‑146 (MQ=255)
                                                ccACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCa                                                                                                    >  2:3627237/1‑150 (MQ=255)
                                                ccACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCa                                                                                                    <  1:1482422/150‑1 (MQ=255)
                                                ccACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCa                                                                                                    <  1:1482420/150‑1 (MQ=255)
                                                    aGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTaa                                                                                                >  1:3135557/1‑150 (MQ=255)
                                                      tCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAgc                                                                                              >  1:1476483/1‑150 (MQ=255)
                                                        ccccTCGGCGGGCTTCACCATCCATTCACGCTTTTTTCAAGCCCTCCCCAGCACAGTGAAGGCTCAACCGCTTCCCTACGTTAGCCCGCATAGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAgcgc                                                                                            >  2:4189180/3‑150 (MQ=255)
                                                        cGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGTACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTGTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAgcgc                                                                                            >  2:5270207/1‑150 (MQ=255)
                                                         gCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCt                                                                                           >  2:3482758/1‑150 (MQ=255)
                                                           cTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCttt                                                                                         <  2:1502244/150‑1 (MQ=255)
                                                             cggcggTCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGc                                                                                       <  2:637363/150‑1 (MQ=255)
                                                             cggcggGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGc                                                                                       <  1:1219414/150‑1 (MQ=255)
                                                                    ctgacccATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAg                                                                                >  2:2313837/6‑150 (MQ=255)
                                                                    cTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAg                                                                                >  1:5327952/1‑150 (MQ=255)
                                                                    cTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGATTTTCTTTGGGCGGCTACGCGTGATAACCCAGGTACTCACTAAGCGATTTGCGTGTCAg                                                                                >  1:3507837/1‑150 (MQ=255)
                                                                     ttCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGa                                                                               <  1:239557/150‑1 (MQ=255)
                                                                           aTGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTAcgc                                                                                                                       >  2:3372031/1‑104 (MQ=255)
                                                                           aTGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTAcgc                                                                                                                       <  1:3372031/104‑1 (MQ=255)
                                                                             gCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCaa                                                                       <  1:5406636/150‑1 (MQ=255)
                                                                             gCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCaa                                                                       <  1:2682280/150‑1 (MQ=255)
                                                                              cATTCACGCTTTTTTCACGCCCTCCCCTGCACTGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAg                                                                      <  1:3227896/150‑1 (MQ=255)
                                                                                ttCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGa                                                                    >  2:2909241/1‑150 (MQ=255)
                                                                                   aCGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAAt                                                                 >  1:5387352/1‑150 (MQ=255)
                                                                                   aCGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAAt                                                                 >  2:4733900/1‑150 (MQ=255)
                                                                                   aCGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAAt                                                                 >  1:4498135/1‑150 (MQ=255)
                                                                                   aCGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAAt                                                                 >  1:3525708/1‑150 (MQ=255)
                                                                                     gCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATAc                                                               >  1:3014848/1‑150 (MQ=255)
                                                                                     gCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATAc                                                               >  1:3200388/1‑150 (MQ=255)
                                                                                     gCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATAc                                                               >  1:2174304/1‑150 (MQ=255)
                                                                                     gCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATAc                                                               >  1:3911714/1‑150 (MQ=255)
                                                                                      cTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGAGGACCCCGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACt                                                              >  2:861495/1‑150 (MQ=255)
                                                                                        tttttCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCaa                                                                       <  1:2570902/139‑1 (MQ=255)
                                                                                        tttttCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCaa                                                                       >  2:2571161/1‑139 (MQ=255)
                                                                                          tttCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGc                                                          <  1:1574344/150‑1 (MQ=255)
                                                                                          tttCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGc                                                          <  1:1774424/150‑1 (MQ=255)
                                                                                            tCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCtt                                                        <  1:861495/150‑1 (MQ=255)
                                                                                               cGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATc                                                                             <  2:763015/126‑1 (MQ=255)
                                                                                               cGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTtggt                                                     <  2:2778783/150‑1 (MQ=255)
                                                                                                gCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTtggtg                                                    <  2:3964415/150‑1 (MQ=255)
                                                                                                   cTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCt                                                 <  1:2837431/150‑1 (MQ=255)
                                                                                                    tCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCa                                                                                 >  1:3258449/1‑117 (MQ=255)
                                                                                                    tCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCa                                                                                 <  2:3258449/117‑1 (MQ=255)
                                                                                                     ccccTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGc                                               <  2:4880289/150‑1 (MQ=255)
                                                                                                     ccccTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGc                                               <  2:3200388/150‑1 (MQ=255)
                                                                                                      cccTGCTCAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGaaa                                                                  <  1:4789348/130‑1 (MQ=255)
                                                                                                      cccTGCGCAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGaaa                                                                  >  2:4789348/1‑130 (MQ=255)
                                                                                                      cccTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGc                                               >  2:1337476/1‑149 (MQ=255)
                                                                                                      cccTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGc                                               <  1:1337476/149‑1 (MQ=255)
                                                                                                           cacaGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCt                                         <  2:1904987/150‑1 (MQ=255)
                                                                                                               gTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGcc                                     <  1:1789617/150‑1 (MQ=255)
                                                                                                                tGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCg                                    >  2:2069199/1‑150 (MQ=255)
                                                                                                                 gAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGc                                   >  2:4239623/1‑150 (MQ=255)
                                                                                                                   aGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTa                                 <  1:151481/150‑1 (MQ=255)
                                                                                                                    ggCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAg                                <  2:4592894/150‑1 (MQ=255)
                                                                                                                     gCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAgg                               >  2:4932076/1‑150 (MQ=255)
                                                                                                                      cTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGt                              <  1:3697554/150‑1 (MQ=255)
                                                                                                                      cTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGt                              <  1:3697598/150‑1 (MQ=255)
                                                                                                                      cTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGt                              >  2:5000005/1‑150 (MQ=255)
                                                                                                                       tCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTg                             >  2:4780792/1‑150 (MQ=255)
                                                                                                                       tCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTg                             >  1:714037/1‑150 (MQ=255)
                                                                                                                       tCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTg                             >  2:2496521/1‑150 (MQ=255)
                                                                                                                          tCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAg                          >  2:1179702/1‑150 (MQ=255)
                                                                                                                          tCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAg                          <  2:1598345/150‑1 (MQ=255)
                                                                                                                           ccGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTCCGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTGTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAgg                         >  2:3468770/1‑150 (MQ=255)
                                                                                                                                  ccTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGcc                  <  2:4498135/150‑1 (MQ=255)
                                                                                                                                   cTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGccc                 >  1:1077206/1‑150 (MQ=255)
                                                                                                                                    tACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCt                >  1:4312648/1‑150 (MQ=255)
                                                                                                                                       aaaaGCCCGCACTGCAATGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTAcc             >  2:2662864/4‑150 (MQ=255)
                                                                                                                                        ttAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTg                                                           >  2:4585960/1‑103 (MQ=255)
                                                                                                                                        ttAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTg                                                           <  1:4585960/103‑1 (MQ=255)
                                                                                                                                         aaGCCCGCACTGCAATGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCt           >  1:1787705/2‑150 (MQ=255)
                                                                                                                                           gCCCGCACTGCAATGCGGGCTTTTCTTTGGGCGGATACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTca         >  1:811207/1‑150 (MQ=255)
                                                                                                                                              cGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTCACAt      >  1:3946806/1‑150 (MQ=255)
                                                                                                                                               gCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTCACAtt     <  2:5327952/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GAGTAGATGGCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTCACATTGAA  >  NC_002947/4416852‑4417147

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: